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Summary 



The results presented in this thesis have in part already been published in Refs. [1, 2, 
3, 4, 5] listed overleaf (page v). We consider physics beyond the Standard Model which 
implies the existence a of long-lived electromagnetically charged massive particle species 
(CHAMP) which we denote by X^. We discuss in detail the unique sensitivity the early 
Universe exhibits on the mere presence and on the decay of such a particle. A CHAMP 
can be realized in supersymmetric (SUSY) extensions of the Standard Model. We carry 
out a detailed study of gravitino (G) dark matter scenarios in which the lighter scalar 
tau (stau, Ti) is the lightest Standard Model superpartner so that ti = X. We also 
provide a thorough investigation of the thermal freeze-out process of ti . 

The thesis is divided into three parts: 

Part I: In this part we consider a generic but weak-scale CHAMP. In Chapter 1 we 
set the stage for the coming investigations by shortly reviewing the framework of Big 
Bang Nucleosynthesis (BBN), by working out the typical CHAMP freeze-out abundance, 
and by reviewing the stringent constraints arising from such a decaying component 
during/after BBN. We also take a critical look at the BBN constraints arising from 
the hadronic decay modes of an arbitrary exotic. In particular, we develop on a refined 
treatment of the Coulomb stopping mechanism of charged hadrons. 

In Chapter 2 we discuss the physics which emerges when the light elements fused in 
BBN are captured by X~ at the time of primordial nucleosynthesis. Since the associated, 
most striking effects were only discovered recently, we provide a detailed exposition of 
the topic. In particular, we explicitly show how to obtain the rates for bound state 
formation which carry a finite charge radius correction of the nucleus. In the remainder 
of this chapter, which is based on [4], we focus on the catalytic production of ^Li and 
^Be. There, we also discuss the issue of a potential late-time catalysis due to proton- 
CHAMP bound states. Upon solution of the full set of Boltzmann equations we obtain 
stringent constraints on the primordial presence of long-lived X~ from overproduction 
of ^Li. Moreover, setting an upper limit on the abundance of primordial ^Be allows us 



iii 



iv 



to constrain this scenario also from catalytic Be production. 

Part II: The second part is devoted to scenarios in which G is the lightest supersym- 
metric particle (LSP) and n is the next-to-lightest SUSY particle (NLSP). In Chapter 3 
we focus on the gravitino LSP as a dark matter candidate. We recollect the results on 
thermal gravitino production, consider explicitly the post-inflationary reheating process, 
and obtain an update on the upper bound on the reheating temperature of the Universe 
from thermal production. 

In Chapter 4 we then focus on gravitino dark matter scenarios in which ri is the 
NLSP. This chapter resembles many of the results of the research papers [1, 2, 3, 4]. We 
constrain the gravitino-stau scenario by incorporating the BBN bounds from ri-decays 
previously obtained in the literature. In addition, the concrete realization of the long- 
lived CHAMP scenario allows us to employ our results on the catalytic production of 
^Be and ^Li. In the framework of the constrained minimal supersymmetric Standard 
Model (CMSSM) a ti NLSP can be naturally accommodated. There, we show that the 
novel catalytic effects severely constrain the reheating temperature of the Universe and 
potentially imply very heavy superparticle mass spectra which will be hard to probe at 
the upcoming Large Hadron Collider (LHC) experiments. We also consider explicitly 
the possibility of a non-standard cosmological evolution and check for the viability of 
thermal leptogenesis. 

Part III: Chapter 5 constitutes the final part of this thesis and is based on [5]. There, 
we take an in-depth look into the chemical decoupling process of the long-lived ri from 
the primordial plasma. The quantity of interest is the thermal freeze-out abundance of 
the stau. We identify its dependence on the crucial SUSY parameters and also show 
that it sensitively depends on the details of the Higgs sector. Stau annihilation into final 
state Higgses as well as resonant annihilation via the heavy CP even Higgs boson can 
substantially deplete the decoupling yield. Remarkably, we find these features are already 
realized in the CMSSM. In those regions of the parameter space even the most restrictive 
bounds from the thermal catalysis of BBN reactions can potentially be respected. We 
discuss the implications for the gravitino-stau scenario. 
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Chapter 1 



BBN and particle decays 



We start this work with a brief introduction into the framework of Big Bang Nucleosyn- 
thesis (BBN) reviewing abundance predictions of some of the primordial light elements 
and discussing their current observational status (Sects. 1.1 and 1.2). In Sec. 1.3 we then 
carry out a simplified treatment of the chemical decoupling of a long-lived CHAMP {X^ ) . 
This frames the thermal X-abundance region and in Sec. 1.3 we shall see that with it 
are associated strong limits on the energy release from X decays during/after BBN. 

In Sec. 1.5 we investigate in some detail the stopping mechanism of injected particles 
short after the main stage of primordial nucleosynthesis. There, we will recover existing 
results of the literature but also develop on a refined treatment of Coulomb stopping of 
injected charged hadrons. 

1.1 Primordial nucleosynthesis after WMAP 

The cumulative evidence from observations of the Hubble expansion as well as of the 
cosmic microwave background (CMB) radiation has put the hot Big Bang model on 
firm footing. In addition, one of the pillars on which modern day cosmology rests on 
is the framework of BBN. Relying solely on Standard Model physics and a Friedmann- 
Robertson- Walker Universe, an overall agreement between the BBN predictions and 
the observationally inferred primordial abundances of the light elements D, ^He, ^He, 
and ''Li is found. This is truly striking given that those elements span nine orders 
of magnitude in number and that light element observations are performed in vastly 
different astrophysical sites. It is this concordance which provides direct evidence that 
the Universe must once have had a temperature T > 1 MeV. 

Standard BBN (SBBN) has only one free parameter, the baryon-to-photon ratio 
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Vh = ny,/n^- It measures the nucleoli content of the primordial plasma and controls 
the rates of the processes which eventually lead to the fusion of the light elements. 
With the measurements of the Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite experiment [6, 7, 8] unprecedented precision data on the multipoles of the CMB 
angular power spectrum became available. Based on a ACDM model, i.e., a flat Uni- 
verse filled with baryons, cold dark matter, neutrinos, and a cosmological constant, this 
has allowed one to pin the baryon density down to [9] Q^h'^ = 0.02273 ± 0.00062 with 
h = Hq/{WO kmMpc^^s^^) parameterizing the Hubble constant Hq. The value implies 
a baryon-to-photon ratio^ of [9] 



so that we have knowledge of the baryon content of the Universe at the time of photon 
decoupling to ~ 3% accuracy (at 68% C.L.). Using (1.1) and/or other non-BBN de- 
terminations of ?7b as input for primordial nucleosynthesis makes BBN a parameter-free 
theory. When we talk about the SBBN light element predictions in the following we 
shall mean the outlined minimal framework of primordial nucleosynthesis with rjy, fixed 
by the CMB measurements. 

1.2 BBN as a probe for New Physics 

The comparison of SBBN predictions with the observationally inferred primordial light 
element abundances makes the theory of primordial nucleosynthesis a powerful tool to 
test and to constrain models of New Physics. 

A true success of the standard cosmological model is the emerging concordance be- 
tween the SBBN predicted deuterium abundance and the measurements of D/H (in 
number) in hydrogen-rich clouds absorbing the light of background quasars at high red- 
shifts. Those astrophysical sites are believed to be most appropriate to yield an estimate 
on the primordial fraction D/H|p. Including the latest measurement [11] of this ratio, 
the weighted mean of seven determinations reads [12] 



Conversely, with an uncertainty which is comparable to that of weak and nuclear rates 
used in BBN codes, the SBBN deuterium abundance can be predicted in the r/b-range 

'^This follows from the WMAP 5-year data set. For comparison, the 3-year result implied 77b(CMB) = 
6.116io'24g [7]. The conversion from Qhh^ to T^b requires knowledge of the average mass per baryon [10]. 



7?b(CMB) = (6.225 ±0.170) x 10 



(1.1) 




(1.2) 
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of interest as [13] 

/6 X 10"^°\ ^'^ 

D/HI = 2.67(1 ±0.03) X lO^M = (2.52 ± 0.11) x 10"^ . (1.3) 

^ \ J 

In the last expression we have used the CMB inferred baryon-to-photon ratio (1.1) and 
added errors in quadrature. As can be seen, both values agree within their ~ la range. 
Indeed, despite the difficult observations, deuterium is the baryometer of choice. Because 
of its weak binding energy, D is only destroyed in astrophysical environments so that 
its post-BBN evolution is monotonic. Moreover, the SBBN prediction shows a strong 
sensitivity on the baryon-to-photon ratio, D/H|p oc f]^^'^ . Any physical process which 
is triggered by extending the SBBN framework must not spoil the agreement between 
prediction and observation. 

Though the agreement in deuterium is impressive it may still only be a coincidence. 
Let us consider ^He which is the most tightly bound element among the SBBN products. 
The primordial ^He mass fraction is defined as^ 



4mHe/"-H 

1 + 4mHe/nH 



yp ^ ^ 0.25 (1.4) 



and with ~ 25% this makes ^He the second most abundant element after hydrogen. The 
estimate in the last relation already follows from the observation that most neutrons 
available are finally bound in ^He and that the neutron-to-proton ratio in number at 
the onset of BBN is n/p ~ 1/7. Observationally, "^He is inferred from helium and 
hydrogen recombination lines measured by now in more than 80 extragalactic HII regions 
of low-metallicity. Following [13], the estimate for primordial mass fraction reads Yp = 
0.240 lb 0.006 where the large adopted error reflects the fact that systematic uncertainties 
may well dominate; cf. [13] and references therein. Though the value is somewhat low 
there is currently no clear discrepancy with its SBBN prediction, the latest one reading 
Yp = 0.2486 ± 0.0002 [14]. It should be noted, however, that ^He is a poor baryometer 
varying only logarithmically with r/b- Contrariwise, being very sensitive to n/p and thus 
to the Hubble rate, ^He acts as a powerful discriminator between models predicting 
additional relativistic degrees of freedom at the onset of BBN. 

Among the most generic ways how physics beyond the Standard Model can affect the 
output of BBN are, e.g., a change in timing of the reactions caused by new contribu- 
tions to the Hubble expansion rate, non-thermal nuclear reactions from late decays and 
annihilation of heavy particles, and the thermal catalysis of nuclear reactions caused by 



^The convention to call Yp the mass fraction is slightly misleading since mria = 3.97mp with nip (ma) 
denoting the mass of the proton (alpha-particle). However, what is observed is Yp and it represents the 
*He abundance in mass within 1% accuracy. 
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electromagnetic or strongly interacting relics. In this regard, the stable lithium isotopes 
have attracted much attention because they turn out to be very sensitive on the latter 
two effects. 

Standard BBN has a long-standing lithium problem. Let us first discuss the heavier 
and more stable isotope ^Li. At r;b(CMB) it is produced mainly in form of "^Be via 
^He + ^He — > ''Be + 7 which then beta decays via electron capture into ''Li after BBN. 
The cross section for the ''Be fusion also dominates the error on the SBBN prediction. 
With a recent update of the reaction cross section [14] the authors tighten the SBBN 
prediction to 

Vh|p = (5.2<^^) X 10-1°. (1.5) 

Lithium is observed in absorption spectra in the atmospheres of metal-poor stars in the 
galactic halo as well as in stars of galactic globular clusters. A link between the measured 
^Li with a primordial origin was first promoted in [15]. What has become known as the 
"Spite-plateau" was an observed constant lithium abundance of ^(Li) = 2.05 it 0.15 
in halo dwarfs of low metallicity^ —2.4 < [Fe/H] < —1.1 and which corresponds to 
^Li/H = (1.12+°;^^) X IQ-^o using ^(Li) = log^o (Li/H) + 12. Ever since much work 
has been done and other groups found similar values so that there seems to be a clear 
discrepancy with the SBBN prediction (1.5) being a factor of a few too high. Indeed, 
the indication of a correlation of ''Li with [Fe/H] [17, 18] tilts the plateau so when 
extrapolating to smallest metallicities values as low as ^Li/H = 6.3 x 10-^^ [19] have been 
inferred. Moreover, such an increasing discrepancy is not alleviated by the most recent 
observation that the ^Li abundance in extremely metal-poor stars with [Fe/H] < — 3 is 
on average 0.2 dex lower than in those (plateau) stars of higher metallicity [20].^ In this 
work we will not touch the ^Li problem. Instead, we concentrate much of our attention 
to the second stable lithium isotope ®Li. 

The measurements of ^Li in the atmospheres of old stars of low metallicity are ex- 
tremely difficult with only one firm detection in the 1990s [21, 22, 23, 24] whereas other 
measurements of ^Li have changed into upper limits; cf. [19] and references therein. 
More recently ^Li has been observed in 9 more halo dwarfs with —3 < [Fe/H] < — 1 
showing a similar isotopic ratio of ^Li/''Li of ~ 5% [18]. This is tantalizing because it 
suggests the existence of a ^Li plateau mirroring the one for ''Li. At first glance, this 
points to a primordial origin of ^Li at the level of ^Li/H|p ~ few x 10"^^. However, the 
story is complicated by the fact that lithium is produced in galactic cosmic rays and 
may as well have undergone stellar depletion. Whereas in standard stellar models ^Li 

^[Fe/H] = log^o (Fe/H) - log^j (Fe/H)^ with the solar abundance logjo (Fe/H)^ ~ -4.55 [16]. 
''dex denotes the decimal exponent. For example, from A(Li) = 2.0 to A(Li) = 1.8 is 0.2 dex. 
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depletion is negligible [25, 26], ^Li is more fragile and particularly destruction in proton 
burning is more efficient. Indeed, non-standard models leading to lithium destruction, 
e.g., from inward diffusion or from rotationally induced mixing have been considered, 
trying to reconcile ^Li observations with its SBBN prediction; cf. [18] and references 
therein. However, the absence of significant scatter in the stars of the Spite-plateau de- 
mands a uniform depletion thus putting strong constraints on any of such mechanisms. 
When considering upper bounds on the primordial ®Li abundance many papers adopt 
values in the range 

^Li/H|p < 10^11 ^10"^°. (1.6) 

Comparing this with the SBBN output ^Li/H|p ~ 10~^^ (see Sec. 2.4.1) this isotope 
shows a gaping discrepancy between prediction and observation; we refer the reader 
to Sec. 4.1 for a further discussion. 

The lithium problem(s) has (have) particularly inspired non-standard BBN scenarios 
seeking their solution. Most notably in this regard are the possibility of the late-decay 
of a massive particle species [see Sec. 1.4] and the catalysis of nuclear reactions; see 
Chap. 2. In this thesis we will exclusively consider physics beyond the Standard Model 
with a weak-scale long-lived CHAMP which we call X^. We shall see that X-decays as 
well the catalysis of nuclear reactions due to the presence of X~ during BBN will pose 
strong constraints on the CHAMP abundance/lifetime parameter space. 



1.3 Typical CHAMP abundances 

Let us assume that the temperature of the primordial plasma was T ^ rnxf^O with 
mx ^ O (100 GeV) denoting the mass of X. Then, X has once been tracking an equi- 
librium abundance. With dropping temperature, X cannot maintain thermal equilibrium 
so that the species freezes-out.^ This happens approximately at the time when the rate 
of X-annihilation drops below the Hubble expansion rate H{T). 

The key to the freeze-out abundance of X lies in considering the Boltzmann equation 
for the total X number density nx = nx+ + 1 

^ + ?,Hnx = -{a^nnv) [nj, - {n'^f] . (1.7) 
The Hubble rate is given by 



^ ' ^' 90 Mp ^ 



"For a low reheating temperature scenario where X may not achieve thermal equilibrium see [27]. 



8 Chapter 1. BBN and particle decays 



with (7cfr radiation degrees of freedom and Mp denoting the (reduced) Planck mass 
Mp ~ 2.4 X 10^*^ GeV. The quantity {a 

ann^) IS found upon a thermal average of the 
total annihilation cross section o"ann times the "relative velocity" v; for details on the 
exact definition of ciann) {<^a,nnv)-, and v we refer the reader to Sects. 5.2 and 1.5.1. The 
equilibrium number density is denoted by n^. 

For a relic species it is customary to scale out the dilution of the number density 
due to the expansion of the Universe. We define the yield variable Yx by dividing nx 
by the entropy density s(T) = 2Tr'^ hes{T)T^ / 45 where hes is an effective degrees of 
freedom parameter [28] . In absence of X-destroying or -producing events and as long as 
no entropy is released, Yx = nx/s remains constant. From (1.7) one then finds 



where [28] 



'dF " V 45 ^P\^a°n^^) [Yx - (^x ) J (1-9) 



gr = ^il + -T ^ • 1.10 



The exact solution of (1.9) can be rather involved and for the case where X is the 
lighter stau, X = ti, this is presented in great detail in Part III of this thesis. Never- 
theless, in order to get a feeling for the expected abundances of an electromagnetically 
charged relic we can employ a simplified treatment of decoupling which is based on the 
non-relativistic limit for X.^ In this limit, the equilibrium number density is given by 



„- = ,,(!^)"%-"--v/r, (1.11) 



and the thermally averaged cross section may be written as [29, 30] 

(aan„T^)n.r. ^ ^ {'^y^'J^^dvv\a^^^v)e-^^-'/'^ (1.12) 

To find the (approximate) decoupling temperature Tf we equate (cTannt^) (-^f) ~ 
H{T[). With the notation Xf = mx/Tf this yields the standard result 

Let us assume an X annihilation cross section expanded in powers of v, (JannV ~ a+bv'^ 
so that with (1.12) (fiann^) develops the form ((Tann?^)n.r. ~ a + 6bT/mx- Choosing 



^We disregard here effects on Yx such as coannihilation, annihilation on the threshold, or resonant 
annihilations [29]; see, however. Part 111. 
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a = o? Ivr?-^ on dimensional grounds''' and considering s-wave annihilation, (6 ^ a) we 
find numerically j;f ~ 26 24 for rax = 0.1 1 TeV and i^x = 2. This gives the 
abundance at the time of chemical decoupling, Yx{Tf) ~ n^(Tf)/s(rf). For T < Tf we 
can neglect in (1.9) so that when accounting for residual X annihilation from Tf to 
To one finds for the inverse of the freeze-out yield 

with To = 2.725 K [32] denoting the present day photon temperature. Whenever we 
write Yx in the following, Yx(Tq) is understood, i.e., the yield of the species X it would 
have had today had it not decayed. 

Taking into the account the temperature dependence of by interpolating the tab- 
ulated values in [33] and integrating (1.14) yields the following estimate on the X abun- 
dance 

with an approximate linear scaling in mx- Note that an increase in {aannv) contributes 
linearly to Yx{Tq)^^, provided Xf = const. Therefore, we have indicated that Yx in (1.15) 
is a value more towards the upper end, corresponding to a guaranteed annihilation cross 
section of electromagnetic strength. A stronger coupling will allow X to stay longer in 
equilibrium, thus receiving an additional Boltzmann suppression. 

Conversely, we can constrain Yx from below by assuming the maximum cross section 
of mutual X^ annihilation which is given by the unitarity limit [34], = vrA^ (s-wave); 
A denotes the de Broglie wavelength of the relative motion. Using A = 1 /{m^cdv) together 
with the reduced mass m^ed = ^xf^ one finds with (1.12) 



(0"ann,ut')n.r. = —\ (1-16) 



4^7? [rnx 

A 

Employing this cross section yields an estimate on the smallest possible freeze out abun- 
dance for a weak scale electromagnetically charged relic. Using mx = 100 GeV gives 
Xf ~ 40 and 

Yx > W-^\ (1.17) 

Let us see how this lower limit on the decoupling yield compares with experimental 
bounds on charged cosmological relics from (negative) searches of anomalous heavy iso- 
topes of ordinary nuclei. For example, in [35] severe limits on the concentration of X~^ 



For example, the cross section for annihilation into two photons reads a^-yV ~ 2na'^/mji [31]. When 
considering the total abundance this gives ffann = CT77/2 and hence {(Js,nnv) = ira^/m^ +0 (rnx/T). 
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in form of heavy hydrogen as well as X~ in low Z-nuclei have been obtained for a weak 
scale relic in the mass range 100 GeV < mx < 10 TeV. It was found that present day 
abundances in excess of 

Yx+x- < 10-'' (1.18) 

are firmly excluded.® Note that individual limits on either or exist for masses 
ranging from a few GeV to multi-TeV and they can be stronger than (1.18) by more than 
ten orders of magnitude; cf. [32]. Comparing (1.17) with (1.18) forces us to conclude 
that, under our assumption of a standard cosmological evolution, X^ cannot be stable. 
Considering a charged thermal relic of finite lifetime with abundances in the range (1.17) 
to (1.15) sets the stage for our further investigations. 

In the next chapter we shall also quantify the X abundance normalized to baryon 
number rib instead of entropy. We define 

X^^ = lA^iZo) ^ 1.13 X 10^0 Yx, (1.19) 

where riy = 2 C(3) T^/vr^ is the photon number density with ({x) denoting the Riemann 
Zeta function. We have chosen this notation in order to clearly distinguish the two 
different normalizations and it will be clear from the context whether X denotes the 
particle itself or its abundance. The previous estimates (1.15) and (1.17) then translate 
into 

'""^^^'""(ra)- f^-^") 

It is also instructive to compare the X abundance with that of "^He. This will be of 
some importance in the discussion of catalytic BBN effects where bound states of ^He 
with X~ play a key role. Since to a very good approximation it follows from (1.4) that 

^He = n4He/?^b = Yp/'i ~ 0.06 (1.21) 

we see from (1.20) that we can expect that the number density of X is typically smaller 
than that of ^He unless X is rather heavy; here, 600 GeV but in concrete particle physics 
models with X annihilating via a number of channels, a heavier X is required. Indeed, 
when focusing on the particle content of the minimal supersymmetric Standard Model 
(MSSM) (plus a gravitino LSP) with ti = X, the X-abundance is determined by the 
standard chemical decoupling and X < '^He holds unless mx > 0(4 TeV); see Sec. 4.1. 



*We have obtained the constraint from the right end of the ^'*C-line in Fig. 7 of [35] and changed the 
normalization of nx from baryon number to entropy. For a recent compilation of other such limits along 
with a thorough investigation of the decoupling yield of a generic electromagnetically- or color-charged 
particle species confer [36]. 
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1.4 Particle decays during BBN 

Using primordial nucleosynthesis as a consistency check for the existence of long-lived 
particles X has a long-standing history.^ When exotics decay during or after BBN elec- 
tromagnetic and/or hadronic energy is injected into the plasma. Depending on timing, 
energy deposition, and abundance of the decaying species, the light element output can 
be affected significantly. The comparison with the observational bounds yields con- 
straints on the parameter space of X. From early works, e.g., [37, 38], to elaborate 
studies [39, 40, 41] this is still an active field of research of continuing refinement and 
increasing sophistication; for a most recent work see, e.g, [42]. Since we will also make 
use of BBN constraints on electromagnetic and hadronic energy release, we provide here 
a cursory overview pointing out important features. For a more detailed exposition of 
the topic we refer the reader to [39, 40, 41] and references therein. 

Electromagnetic constraints The emerging constraints can be classified with respect 
to the decay mode of the exotic particle. When X decays radiatively into primary 
high-energy photon(s) and/or electrons (charged leptons) an electromagnetic cas- 
cade is induced. The important processes are pair creation (7 + 7bg — > e~ +e"'"), 
photon-photon scattering (7 + 7bg 7 + 7), Compton scattering (7 + ebg ^ 7 + e), 
inverse Compton scattering (e+7bg e+7), and pair creation on nuclei (7+Nbg — > 
e"*" + e~ + N). The subscript 'bg' denotes the particles which are in equilibrium 
with the plasma. Since n^/ne|bg ^ lO^'' the scattering on background photons 
is very frequent. This leads to an efficient thermalization of the cascade so that 
destruction of light elements does not happen frequently. However, once energetic 
photons are degraded below < m1/22T [43] they loose their ability to pair 
create on 7bg. The soft photons of the associated 'break-out' spectrum are then 
capable to efficiently destroy those light elements whose binding energy E^'^ lies 
below the threshold of electron pair creation. For D {E^^ = 2.22 MeV [44]) this 
happens at T < 10 keV whereas "^He (^{J^^ = 28.3 MeV [44]) is destroyed when 
the photon temperature drops below T < 1 keV. This corresponds to respective 
cosmic times of t > 10^ s and t > 10^ s when thermal nucleosynthesis reactions 
have long frozen out. 

Constraints on the electromagnetic energy release for mx = 1 TeV are shown in 
the left Fig. 1.0 which is taken from [41]. The lines represent upper limits on the 
quantity E^isYx, i.e., on the "visible energy" E^is released per decay times their 

^In this section X stands for an arbitrary, not necessarily electromagnetically charged, long-hved 
species. 
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Figure 1.1: Constraints on tlie electromagnetic {left figure) and hadronic {right 
figure) energy release from X-decays for mx = 1 TeV; the figures are taken from 
[41]. The lines represent upper limits on the quantity E^isYx, i.e., on the "visible 
energy" i^vis released per decay times their abundance prior to decay, Yx, and 
are plotted as a function of tx', see main text for discussion. 



abundance prior to decay, Yx , and are plotted as a function of tx ■ Two dominant 
constraints are visible. For tx ^ 10^ s the most restrictive constraint, labeled 
D/H, arises from the destruction of deuterium below its observationally inferred 
primordial level. For larger tx, however, ^He gets dissociated and D (along with 
^He and T) is also created. Since the ^He-target is very abundant, D is indeed 
overproduced for tx ^ 10^ s. Moreover, the combination ^He/D is then always 
produced yielding the most stringent constraint on electromagnetic energy release 
for TX > 106 s. 

Hadronic Constraints A second class of constraints on decaying X during/after BBN 
comes from hadronic energy release into the plasma. For example, even if X 
dominantly decays into photons, a non-vanishing hadronic branching ratio is ex- 
pected from the conversion of a (virtual) photon into a quark-antiquark pair or 
from charged meson production on background photons (if kinematically allowed). 
The partons emitted in the decay are quickly hadronized and the highly energetic 
fragmentation products such as protons (p), neutrons (n), as well as their antipar- 
ticles are released. Also long-lived mesons, namely, charged pions (vr^) and kaons 
{K^, K^J, with lifetimes O (lO~^ s) have a chance to interact with background 
nuclei before decaying. Once an energetic hadron scatters on a background nu- 
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cleus, essentially p or He, a hadronic shower is induced. In particular, He may 
be destroyed with secondaries further participating in interactions with the plasma 
constituents. 

The energetic charged hadrons N are downgraded in energy via electromagnetic in- 
teractions, most importantly, by Coulomb scattering (N+e^g — > N+e^), Compton 
scattering (N+7bg N+7), and Bethe-Heitler scattering (N+7bg — > N+e^+e"). 
Injected neutrons loose their energy mainly by their magnetic-moment interaction 
with e^g. It is clear that emergent constraints on the hadronic energy release will 
sensitively depend on the competition between the rate for hadronic scattering and 
the rates for (electromagnetic) thermalization and/or decay (of unstable particles). 
Moreover, even after hadrons are stopped they may still induce neutron-to-proton 
interconversion processes [37]. 

In the right Fig. 1.0 which is taken from [41] the constraints on i^vis^x due to 
hadro-dissociation as well as n-p interconversion are shown for a particle with 
mx = 1 TeV and hadronic branching ratio i?h = 1- The effects from photo- 
dissociation are not included. Note that this is an unrealistic situation since the 
hadron stopping process itself as well as meson decays induce electromagnetic show- 
ers. For Tx < 100 s, i.e., for T > 100 keV, the emitted hadrons essentially deposit 
all their kinetic energy electromagnetically before interacting with the background 
nuclei. However, interconversion processes which always lead to an increase of 
n/p enhances the ^He output. The associated constraints from "^He overproduc- 
tion for two different observationally adopted limits on the primordial mass frac- 
tion are shown by the dotted lines labeled Yp. For larger lifetimes, tx ^ 100 s 
(r < 100 keV), mesons typically decay before interacting hadronically. However, 
the stopping power for protons and neutrons rapidly decreases with dropping tem- 
perature so that ^He is destined for being destroyed. This yields the stringent D/H 
constraint on hadronic energy release for tx ^ 100 s. Moreover, a small fraction 
of the energetic spallation products T and ^He can scatter again on ambient ^He 
producing ^Li [45] (and ''Li). This non-thermally induced fusion reaction gives the 
hadronic constraint labeled ^Li/H in Fig. 1.0. Since ^Li is efficiently destroyed in 
(thermal) proton burning for temperatures T > 10 keV the constraint becomes the 
dominant one only for tx ^ lO'^ s [46]. 



So far, the discussion has been completely generic with our ignorance on the nature 
of X parameterized by i?vis- Constraining the particle's parameter space requires the 
specification of the couplings of X to Standard Model particles as well as its mass mx- 
This allows for the determination of the decay modes of X along with the computation 



14 Chapter 1. BBN and particle decays 



of the associated average electromagnetic and hadronic energy E'vis released per decay. 
Moreover, the freeze-out abundance Yx can be calculated so that plots like in Fig. 1.0 
can be employed to constrain the model. In part II of the thesis we incorporate the most 
stringent of the constraints for the case of a decaying stau in the gravitino dark matter 
scenario. There, we also provide more details as soon the problem of inclusion of such 
constraints becomes acute. 

1.5 A critical look at hadronic constraints for T < 100 keV 

In the previous section we have noted that for T < 100 keV {tx ^ 100 s) the BBN 
constraints on hadronic energy release sensitively depend on the competition between 
the hadronic and electromagnetic scattering rate (and potentially the lifetime). If in- 
jected nucleons as well as their secondaries — either produced in spallation reactions or 
"up-scattered" in elastic scatterings — are predominantly thermalized by interactions on 
background nucleons or nuclei, those constraints become stringent. Underestimating 
the stopping power due to electromagnetic interactions would lead to overly restrictive 
bounds on the hadronic energy release. 

For example, in the last section we have seen that the non-thermal production of ^Li 
due to the energetic spallation debris T and ^He of destroyed ^He yields the dominant 
hadronic constraint for tx ^ 10^ s. The reactions involved are T + '^Helbg — > ^Li + n 
and ^He + ^He|bg ^ '^Li + p [45]. For T < 30keV, i.e., for t > 10^ s, inverse Compton 
scattering on background photons cannot prevent low- energy hadronic interactions [37] 
above the lithium formation threshold [0(10 MeV)]. The dominant electromagnetic 
degradation mechanism is then Coulomb scattering. However, the rapidly diminishing 
number of background electrons (positrons) with dropping temperature also renders the 
energy loss by Coulomb scattering increasingly inefficient. Furthermore, in [40] it is 
claimed that the non-thermal ^Li output is boosted by a factor of ten because of a 
peculiarity in the Coulomb process: Once the velocity f3 of the energetic mass-three 
nuclei drops below the thermal electron velocity (/3e), the stopping power seems to be 
strongly suppressed. This observation was first made in [37]. 

In light of these comments a critical look on the Coulomb stopping process is war- 
ranted. We shall pay particular attention to the velocity dependence of the cross sections, 
i.e., on /3 and {(3e). In Sec. 1.5.1 (and partly also in Sec. 1.5.2) results from the literature 
are reconciled. In Sec. 1.5.3 we focus on the stopping of charged hadrons and incorporate 
the proper screening-prescription of the Coulomb interaction. In Sec. 1.5.4 we discuss 
the obtained results. 
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1.5.1 Energy transfer in binary collisions 

A thorough investigation of the electromagnetic stopping of hadrons in the context of 
primordial nucleosynthesis has first been presented in [37]. Indeed, the treatments in [47, 
41] (see Fig. 1.0) employ the results of that work. Since the stopping power sensitively 
depends on the velocities of the incident hadron and the target particles, we first reconcile 
the general result on the energy transfer obtained in [37]. Though we encounter minor 
disagreements they turn out to be without relevance. 

Our starting point is the rate of energy loss due to binary scatterings (Al) of Ref. [37] 
^ = E (£)3 / ^^^^'P.- /.(p.) [1 T /./(p./)] ^^^-^^-^1 ' (1-22) 

where AEj denotes the energy transfer between the incident hadron and a (background) 
particle species j with three-momentum pj and gj internal degrees of freedom. The 
transfer is weighted by the center-of-mass (CM) cross section duj/dD. and averaged over 
inital and final state distribution functions fj and respectively. A subtle point is 
the appearance of the M0ller velocity [28] 

v^,, = = [(/3 - (3,f - (/3 X /3,)]'/' (1.23) 

which is the relativistic generalization of the conventional relative velocity frei = 1/3 — 
The respective velocities of the hadron and the target are given by /3 = p/E and j3j = 
Pj/Ej and F = [{p ■ Pj)"^ — m^M-^]^/^ denotes the Flux-factor. Only in the CM frame or 
in the rest frame of one of the incoming particles v^0i coincides with v^^i- We stress that 
p = {E,p)'^ and pj = {Ej,pj)'^ denote the respective four-momenta of the energetic 
nucleus and of the ambient target particle in the rest frame of the thermal bath. In that 
frame, and when j is in thermal equilibrium, the distribution functions fj take on their 
familiar form: = [exp (Ej^^j^-^/T) it 1]^"*^. The upper signs in (1.22) and in the last 

expression refer to fermions, the lower to bosons. Finally, Ejf is the energy of the target 
after scattering and rrij (M) is the mass of the target (indicent nucleus). 

The energy transfer AEj = Ejf — Ej can be obtained by a series of Lorentz trans- 
formations: Since the scattering is elastic, in the CM frame we have E'^™ = E^^- Thus, 
we can obtain Ejf by a Lorentz transformation A = A3A2A1 of pj into the CM frame 
followed by an inverse transformation of back. A is broken up as follows:^" We 
choose f3 to lie along the z-axis and to have an angle a with Pj £ yz-plane 

p = {E, 0, 0, pE)'^ and pj = {Ej,0, (3jEj sin a, (ijEj cos af (1.24) 
^''The explicit forms of Ai, A2, and A3 are given in the Appendix l.A at the end of this chapter. 
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so that a = TT corresponds to a "head-on-head" colhsion. Boosting into the rest frame of 
the incident nucleus gives 

f-ff^Ejil- PPj cos a)\ 

.J". =(^i,P;-f (1-25) 

\l(sEj{Pj cos a - P) J 

where 7^3 = (1 — Under Lorentz transformations the M0ller velocity changes 

as [281 



Pj 



Aipj 



1 



(1.26) 



which can be used to obtain the velocity of the ambient target in the rest frame of the 
incident particle. Confirming the expression given in [37] it reads 



P2 ^ 1/3' I 



1 — P(3j cos a 



(1.27) 



By the same token, the expression for the angle ip between (3^ and the ^;-axis reads 



cosip = (3'j ■ e'^/f32 



P — Pj cos a 



(1.28) 



which differs by a sign from [37] . Instead of explicitly carrying out the rotation pj = A2P'j 
which makes (3" parallel to the z'-axis (vm0i is lengthy) we use our knowledge on the 
form oip": p'- = {Ej,0,0, P2Ej)'^ since Ej = E'J. Boosting into the CM frame using A3 
one finds 

fjcmE'jil - P2pcm)\ 




\ 7cmEjiP2 - Pern) J 

where Pcm is obtained from (p^™)^ = —{p'^'^)z', 7c 



Asp',' 



P 



I 7cmM \ 




\-lcmPcmM J 

(1 - PlJ-^'"^. We find 



(1.29) 



Pc 



E',P2 



P2Ej^p{l - PPj COS a) 



(1.30) 



M + E'j M + -f/3Ej{l - PPj cos a) 
In the CM frame, the scattered target three-momentum p^™ has a scattering angle 9 with 



p^™ and both momenta span a plane with azimuthal angle (p. Thus, p^™ 



is given by 



( lcmE'j{l- p2Pcm) \ 
IcmE'j {P2 - Pcm) sin 9 COS c 

IcmE'j {P2 - Pcm) sin 9 sin (, 

\ lcmE'j{P2 - Pcm) COS 9 J 



cm _cm"iT 



(1.31) 



1.5. A critical look at hadronic constraints for T < 100 keV 17 



and with pjj = A ^p^™ = A^ A2 A3 we can transform back into the rest frame of 
the thermal bath. This yields for the energy of the scattered background particle 

Ejf = £'j7^7cm(l - PPj cos a) 

X {7cm [(1 -/3cm/32)(l +/3/3cmCOS^/;) 

- (/3cm - P2){f3 COS -0 + /3cm) COS 6] 

+ /3(/3cm -/32)sin(9sin(/)sinV} (1.32) 

from which lS.Ej = Ejj — Ej is obtained. We encounter some sign-differences and a 
different angular dependence on cj) in the last line with (A2) of [37]. This may be due to 
a different definition of the coordinate system and turns out to yield the same results. 
Since the target medium is unpolarized, daj/dQ, is independent of (j) and the integration 
over the azimuthal angles in (1.22) can be performed upon which the last line of (1.32) 
drops out. Neglecting the Fermi blocking/Bose enhancement factors in (1.22) we find 
(adopting the notation of [37]), 

^ = lZ§^£dPjrnp] {l-f3])-'f,{P„T)I,iP„(3,T), (1.33) 

/I 
d9 da sine sinaAj-^ , (1.34) 

Aj = Wm0i|7|7cm(l - PPj cos a) 

X [(l-/3cm/32)(l+/3/3cmCOS0) 

-(/?cm-/32)(/?COS0 + /3cm)cOS0]-l} . (1.35) 

1.5.2 Hadron-electron scattering 

Let us now focus on "Coulomb scattering" between an incident hadron and background 
electrons (positrons) and compute daj/dQ for j = e^. Note that also neutral hadrons 
scatter on via their magnetic moment interaction. 

For spin-1/2 hadrons such as nucleons or T and ^He nuclei the hadron-photon vertex 
can be written as [48] F^, = 2M{Ge - Gm)Pt,/P'^ + Gm^Tf, with P = p^"^ + and 
pj^ denoting the (outgoing) four-momentum of the nucleus. The respective electric and 
magnetic form factors Ge and Gm depend on the (squared) four-momentum transfer 
1^ — (P/™ ~ j,cm^2 normalized such that Ge(0) = Z is the charge in units of e 

and that Gm(0) = fi is the magnetic moment in units of e/2M of the hadron. The 

^^The Sachs form factors Ge and Gm axe convenient because no interference terms oc G^Gm appear 
in the cross section (1.36); they are related to the Dirac and Pauli form factors Fi and F2 via Ge ~ 
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(unpolarized) differential cross section for electron-hadron scattering is readily obtained. 
In accordance with [48] (typo in [37]) it reads using the Mandelstam variables s, t, and u 

dai/2 -Kop- 1 



dt [s - (M + me)2][s - (M - me)2] t2(i _ t/(4M2)) 
X \Gl{{s-uf ^{^M'^ -t)t\ 



t 



4M2 



Gl,[{s -uf - {AM^ -t){4ml+t)]^ . (1.36) 



Owing to a different vertex structure for spin-0 hadrons such as pions or ^He, 
FPfj^, where F is the electromagnetic form factor {F(0) = Z), one readily obtains [48] 

dao 7ra^F^[{s-uf + {4M^-t)t] 



dt [s- (M + me)2][s- (M-me)2]t2- 



(1.37) 



We can expand (1.36) and (1.37) in terms of x = -ffsEj/M (typo in [37]). The 
expansion is most likely to fail for scattering of (light) ultra-relativistic nuclei at high 
temperatures of the thermal bath. To see the validity of the expansion consider the 
typical energy of an electron (Ee) = (Ej) by using Maxwell-Boltzmann statistics: 

, „ > „ K^ imp IT) 

{E, =3T + m ' 1.38 

Here, is the modified Bessel function of the first/second kind. For example, with 

T = 0.1 MeV, a kinetic energy Tkin = (7 - 1)M = 100 GeV (1 GeV) of the nucleus, and 
M = nip one finds (x) ~ 0.08 (0.002). Thus for the cases of interest the expansion in x 
is fine. Since the cross-sections are independent of the azimuthal angle^^ 

^ = P±^ with p, = Ip^'"! ~ M/32(l - /3/3eCosQ)x (1.39) 
d\l TT dt 

where Pe = Pj- Neglecting m-g it follows 

s - ~ 2M2x(1 - p(3e cos a) 
s - n ~ 4M^x(l - Pf3e cos a) 

t ~ -2M2/?|x2(l-cos6l)(l-/?/3eCosa)2 



so that we find for the CM cross section for charged spin- 1/2 and spin-0 nuclei 

ddch _ a^Z^ 1 -/?|(1 -cosg)/2 

Hn ~ M2x2 /?2(l -/J/3^cosa)2(l -cos0)2' 



(1.40) 



Fi + F2g^/4M^ and Gm = Fi + F2 [49]. With the definitions for Fi and F2 and for > rrie, in the 
laboratory frame, the Rosenbluth formula [50] follows from (1.36). 

^■^An overall sign has been dropped since it can be fixed by the integration borders. 
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Note that we have made the approximation F,Ge — Z assuming small Let us see if 
this is justified. The maximum energy transfer is realized in a back-to-back collision in 
the CM frame for which q"^ = t^i^ = — 4p^. Considering the case that (3 ^ Pj, i.e., the 
case when the electron is a stationary target, it follows from (1.39) that ~ 7^/3me. 
For example, for a proton with E = 100 GeV this gives 

Q2 = _g2 ^ 0.01 GeV^ for which 
GeiQ"^) is practically unchanged from unity [51]. Moreover, note that setting F,Ge = Z 
usually leads to an overestimation of the cross section with F and Ge decreasing for 
increasing [51, 52]. Consequently, the stopping power is overestimated leading to 
more conservative BBN constraints. For neutral hadrons we find^^ 

da^, ^ o?G\ l + /?|(l-cosg)/2 

dVL 2M2 1-COS0 ■ ^ ' ' 

We disagree in (1.40) and (1-41) with [37] by a factor of /3| in the denominator. The 
disagreement arises as follows: Eq. (A5a) of [37] is actually Eq. (139.5) of [48]. In order 
to arrive at the latter equation ~ (E'^™)^ has been used (pg ~ e\ in the notation of 
[48]). However, it is more accurate to use = (7cm/?cmAf)^ = /?! (-Eg™)^ + O [x^^ ; recah 
that /32 is the velocity of the electron as seen from the rest frame of the nucleus. We 
note in passing that up to corrections O (a;^) one has 7cm — Ij /3cm — x(32{^ — (3(3e cosq), 
E^^ ~ M, and E^"^ ~ Mx{l - P^ecosa). 

We can use the above expansion in x to simplify A in (1.35) for the limiting cases of 
ultra-relativistic and non-relativistic hadrons traversing the background plasma. Consid- 
ering P — 1, i.e., an ultra-relativistic incident particle, and therefore 7^3 ^ 1, P2 cosip ~ 
—1, and Um0i — (1 — PPe cos a) we get to leading order O [x^) 

Af' ~ 7^(1 -/?e cos a)2(l- cos 0) . (1.42) 

Conversely, for /3 <C 1 and therefore 7^ ~ 1 + /?^, /?2 — (/?e — /3cosa)/(l — PPeCosa), 
and v^0i ~ /?£ — /3 cos a we find 

A^'^^ {p^ - PPeCosa){pe - Pcosa){l - cos9) . (1.43) 

Both expressions agree with the ones obtained in [37] with differing signs in (1.35) being 
compensated. 

1.5.3 Cutoff considerations for charged particles 

After having obtained the cross section for Coulomb and magnetic moment scattering for 
hadrons on electrons we make the following observation for charged particles: Though 

^''For example, Gm{0) = —1.91 [32] for the neutron, being entirely anomalous. 
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the energy transfer in a cohision is smallest in the forward direction — > [as can be 
seen by the factor (1 — cos6) in Eqs. (1-42) and (1.43)], the divergence (1 — cos6)~'^ in 
the cross section for charged hadrons (1-40) — arising from the long-range nature of the 
Coulomb interaction — is too strong to be canceled. In this sense, the energy loss due 
to scatterings in the forward direction gives the most efficient contribution. Cutting off 
the angular integration in (1.34) at ^min leads to the well known logarithmic dependence 

on fjnin- 

Of course, ^min has to be motivated. 

In a plasma, i.e., in a gas of charged particles, correlation effects lead to the screening 
of the long-range Coulomb interaction. In Ref. [41] the authors determine the cutoff by 
comparing the energy transfer to the electron with the plasma frequency 

^pi = 1-44 

where n^, denotes the total electronic density 

2 X {mlT/'K'^)K2{m/T) for T > me/26 
7/8 rj-Bfi^ for T < me/26 



TT-e = rip- + n„ 



(1.45) 



In the first line we have neglected the electron chemical potential and in the second 
line we have imposed charge neutrality of the Universe. The upper relation in (1.45) is 
derived by using Maxwell-Boltzmann statistics. For T ^ me, i.e., for ultra-relativistic 
electrons/positrons, this implies an error by a factor of 3C(3)/4 ~ 0.9. Note also that 
electrons freeze out in the temperature region of interest, me/26 ~ 20 keV. 

The plasma frequency does, however, not provide the correct scale [53]. The screening 
of the electric field is a longitudinal phenomenon whereas the notion of the plasma fre- 
quency as an effective photon mass is associated with transverse plasma excitations. 
Electrons as well as the (ionized) light elements contribute to the screening with a 
scale [54] 

;2 ; 2 , , 2 47rane Ana ^-^ 2 n Act\ 

j 

where /cd denotes the Debye scale with Debye length Ad = k^^; nj denotes the number 
density of nuclei with charge number Zj. Note that kj) and Wpi can be very different 
with ajpi/kj) = y^T/me. However, in the temperature region of main interest ujpi and 
kj) are within a factor of a few. Moreover, since the screening scale will enter only 
logarithmically we neglect the contribution of the ions (in particular protons) in the 
following and set ks — kn- 
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We shall distinguish two cases: For (3 > (/3e) the electrons can be viewed as a sta- 
tionary target. Thus, we follow the screening prescription obtained in [53] and replace 
(1.36) via 



da 



fis,t,u,M) 



da 
'dt 



fis,t,u,M) 



for P>{Pe) 



(1.47) 



q2(q2 + A;2) 

(for elastic scattering in the CM frame t = — q^) whereas for (3 < (/3e) electrons have time 
to rearrange so that the scattering resembles one on a Yukawa-like charge distribution 
with screening length Ad- Then the correct prescription reads 



da 
'dt 



f{s,t,u,M) 



q- 



da 
'dt 



f{s,t,u,M ) 
(q2 + fc2 



2 \2 



for (3<m- 



(1.48) 



Given the above considerations we replace the scattering cross section (1.40) for 
charged particles employing the screening prescriptions (1.47) and (1.48) for the respec- 
tive cases /3 > {(3e) and {3 < (/3e)- We find 



da ch 



dn 



2 72 



l-/3|(l-cose)/2 



dO. 



^>^^^^ j;2jV^2^2(l-/3/3eCOSQ)2(l-COSe)2 + (l-COSe)K2 ' 

1 -/3|(1 -cos0)/2 



/3<(/3e> 



with 



^^^^ [/3|(1 -/3/9eCosa)(l -cos^) + k2/(i _/3/3^cosa)]^ 
^2 - 4/(2M2x2) = kll(2^lE]). 



(1.49) 
(1.50) 



K 



In the region where acts as a regulator, i.e. 
observation that 



y^rel/nrel 



for e 



(1.51) 

0, we make the immediate 







1 


d^ 






da ch 




1 — COS 9 


dn 







For scattering in the forward direction the screening prescriptions will yield a numerical 
difference only for (1 — cos 6) < . For our cases of interest is usually a very small 
quantity, e.g., -logK^ ~ 0(10^ 15) for T = 30 keV or 0(3^8) for T = 300 keV. 
Thereby, only in a very small integration regime over 6 both cross sections will be 
significantly different — though the integrand is largest in this area. 

In order to decide which cross section is applicable for a given value of the hadron 
velocity /3, it has to be compared with the average electron/positron velocity 

2T{me + T) 



m'lK2{me/T) 



-me/T 



(1.52) 
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which is obtained by using Maxwell-Boltzmann statistics; for T <^ rrif, the formula 
reduces to the standard result (/?e)n.r. = ^/ST /irme. Note that /3 is related to the kinetic 
energy of the incident hadron via /? = [1 — M^/(Tkin + M)^]^/^. Thus, for example, a 
proton with Tki^ = 1 GeV (50 MeV) has (5 = 0.88 (0.31) so that /? drops below {[5e) for 
r = 400 keV (20 keV). 



1.5.4 Discussion on Coulomb stopping 

In this section we discuss in some detail the results on the energy loss for charged 
particles due to Coulomb interactions with the background electrons (positrons). We 
will also compare with a treatment found in the literature. 

To see the net effect of the different screening prescriptions on the stopping power we 
perform a full numerical integration of (1.34) and (1.33) using the Vegas algorithm [55]. 
For the integration over the electron (positron) velocity knowledge of the distribution 
function /g is required. Though electrons are frozen out for T < me/26 ~ 20 keV they 
remain tightly coupled to the photon bath via Thomson scattering. This ensures that 
electrons maintain kinetic equilibrium so that we can make the approximation 



fe{Ee,T) 



[expiEe/T] + l]-i for T > me/26 

(1-53) 

neff^/n"^ for T < me/26 



where T denotes the photon temperature. For T < me/26 we use ff^ = exp {—E(./T) ~ 
exp[— Pe/(2meT) — me/T], i.e., we resort to Maxwell-Boltzmann statistics in the non- 
relativistic limit. From the definition Ue = ge f f^'^Pe/(2vr)^/e we reproduce the second 
line of (1.45) by using in the form of (1-11) with = 2; for the case T > me/26 we 
use ge = 4. 

In Fig. 1.1 we show the stopping power dE/dt for an injected proton computed by 
numerical integration from (1.33) using the different screening prescriptions. Solid lines 
correspond to (1.49) and dotted lines (hardly visible) are associated with (1.50). In the 
left panel we show —dE/dt in units of MeV/s as a function of the proton velocity /3 
at temperatures T = 100 keV and 10 keV as labeled. In addition, the points (3 = {j3e) 
indicate which screening prescription should be used. In the right panel we show —dE/dt 



^*It is shown in [56] that fe = R~^Te^^^ Nq exp {—\>^/2m^) satisfies the Boltzmann equation in the 
non-relativistic hmit with an elastic coUision term due to Thomson scattering; A'o oc n^R^ — const 
by comparison with (1.45). Defining the temperature of a non-relativistic particle species j with ar- 
bitrary distribution / as {3/2)Tjnj = gj J d^pj /{2Tv)^p'j/{2mj)f{pj) [56] it is found that the electron 
temperature Te tracks T weU until recombination, T ~ 0.3 eV. There, (T - T^)/T = O (lO"'') [57]. 
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Figure 1.2: We show the stopping power dE/dt for an injected proton computed by 
numerical integration from (1.33) using the different screening prescriptions. Solid lines 
correspond to (1.49) and dotted lines (hardly visible) are associated with (1.50). Left: 
Shown is —dE/dt as a function of the proton velocity (3 at temperatures T = 100 keV 
and 10 keV as labeled. In addition, the points (3 — (/3e) indicate which screening 
prescription should be used. Right: —dE/dt as a function of T for a relativistic proton 
(Tkin = 1 GeV, P = 0.9) and a non-relativistic proton (Tkin = 50 MeV, P = 0.3) is 
plotted. In addition, the dashed line shows the energy loss due to the density effect in 
the non-relativistic case. 

as a function of T for a relativistic proton (Tkin = 1 GeV, /3 = 0.9) and a non-relativistic 
proton (Tkin = 50 MeV, (3 = 0.3). 

From Fig. 1.1 we can make a number of observations. An immediate one is that the 
stopping power is essentially insensitive to the employed screening prescription. Con- 
cretely, we find that both prescriptions yield a difference in dE/dt by no more than 
20% for the considered temperature range. From the right panel we see that once the 
velocity (3 of the incident particle drops below the average electron velocity, the stopping 
power starts to decrease rapidly. This confirms the observation made in [37] and it is 
also intuitive since it becomes increasingly difficult to transfer momentum to the — on 
average — faster electrons. Indeed, for (3 < f3e the charged hadron can even gain energy 
in a collision which is indicated by a sign-flip of le [Eq. (1.34)] in the collinear region 
where a — > ('head-on-back collision'). From the left panel we realize that the stopping 
power rapidly decreases with dropping temperature. This is because for T < mg the 
number density of electrons and positrons is Boltzmann suppressed. For T < 20 keV 
the decrease is weaker because the remaining electrons fail to track their exponentially 
decreasing equilibrium abundance. More precisely, dE/dt scales like for such low 
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temperatures because rig ~ a ^ and a ~ T ^ during radiation domination. 

So far, we have only considered (screened) binary collisions of a fast charged parti- 
cle traversing a QED plasma. Such a treatment gives an accurate description for those 
scattering events of the particle with largest energy transfers, i.e., with smallest im- 
pact parameters b. Considering b > Xj), the nucleus scatters simultaneously on many 
electrons. (Note that in our case the Debye length is much larger than the typical 

— 1/3 

inter-particle distance, Xd ^ Ue .) A sweeping external charge, i.e., a perturbation 
Pext = Ze6{r — (3t), induces a macroscopic electric field E in the medium which acts 
back on the particle. The resulting energy loss per unit path length can be found by 
computing the work done on the particle. It equals the force exerted onto the charged 
hadron in direction opposite to its motion 



dE 
dx 



_ 1 dE 

kob>i ^ dt 



= Zee^-B{r = I3t), (1.54) 

fcDb>l 



where is a unit vector in (3 direction. The electric field can be found by consid- 
ering the macroscopic Maxwell equations with dielectric permittivity e. In the non- 
relativistic limit and using^^ ^{^) — 1 — Wpj/w^ for a Maxwellian plasma, the stopping 
power reads [59] 



dE 
'dt 



u;'^i\n(- ^ (1.55) 



As can be seen by the dashed line in the left Fig. 1.1 the contribution due to the 
'density effect' is subleading. Note that the formula was derived under the premise 
that the velocity of the massive particle is large compared to the thermal speed of the 
electrons. Indeed, the argument of the logarithm in (1.55) is greater than unity only for 
P ^ 0.6(/3e)n.r. SO that the line is cut-off for T > 60 keV. We remark that a full relativistic 
treatment of the energy loss of a massive particle due to the dielectric response of the 
medium for arbitrary velocities is complex but will not affect significantly the above 
made conclusions; we refer the reader to Landau's treatment in [58].^^ 

We can compare the full numerical integration of (1.33) with the treatment of Coulomb 
stopping presented in [41]. The authors employ the results of [37] which also we have 
taken as a starting point. Full numerical integration of (1.33) is not feasible when scan- 
ning the (rx, Yx) parameter space so that analytical approximations based on (1.42) and 



^^Here, uj and k (fc = jkj) are the frequency and wave vector of the Fourier transformed fields and the 
expression is the limiting case for uj/{k{f3e)) ^ 1; see [58]. 

^^The case of a hot QED plasma with me <^ eT has been treated within the framework of thermal 
field theory in [60] . 
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(1.43) have been used in [41]. Since we have observed that the employed screening pre- 
scription affects dE/dt only marginally for T < 100 keV and that k-^ and ojp\ — entering 
the stopping power logarithmically — are not too different, it is not surprising that we 
find overall agreement with [41] on the energy degradation rate within a factor of a few. 

We remark that obtaining constraints on the hadronic energy release of decaying X 
involves a fair amount of modelling and computation. After calculation of the hadronic 
branching ratio in the decay of X, each step involves uncertainties and approximations: 
Employing a hadronization algorithm, computing the initial energy spectra of secon- 
daries, following the energy degradation and cascade formation due to electromagnetic 
and hadronic processes, and finally obtaining the yields of non-thermally produced light 
elements. We have seen that already the seemlingly elementary process of Coulomb stop- 
ping can become involved — especially when it is necessary to apply it to a large range 
of incident particle energies and plasma temperatures. In the light of these comments 
we close this chapter by noting that we have not found a radically different picture than 
that of previous considerations which would strongly influence on the strength of the 
hadronic constraints presented in Fig. 1.0. 
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The explicit matrices for the Lorentz transformations performed in 1.5.1 are given below. 
The matrix A2 describes an (active) rotation in counter-clockwise direction around the 
x'-axis when looking towards the origin. The inverse transformations Aj~^, A^^, and A3 ^ 



are obtained by the replacement (3 
A3, respectively. 



Ai 



-il), and j3c 



A. 



-/3cm in Ai, A2, and 
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(1.56) 



(1.57) 



26 Chapter 1. BBN and particle decays 



Chapter 2 

Bound states and catalysis of 
BBN 



In this Chapter we now discuss (some of) the rich physics which emerges when the hght 
elements are captured by during/after the time of primordial nucleosynthesis. We 
start in Sec. 2.1 by reviewing the basic properties of such bound states. Section 2.2 is 
devoted to the calculation of the wave functions associated with the CHAMP-nucleus 
system. This allows us in Sec. 2.3 to obtain recombination cross sections carrying a 
finite nuclear charge radius correction. The detailed exposition in Sects. 2.2 and 2.3 is 
of some value since (apart from exceptions) those rates are not publicly available in the 
literature. 

In Sec. 2.4 we then consider the catalysis of BBN reactions. After a general review 
we employ the results from the literature on the catalyzed production of ^Li and ^Be 
and show explicitly how to incorporate them into a Boltzmann network calculation. In 
Sec. 2.5 we discuss the potential impact of neutral proton-CHAMP bound states on the 
synthesized elements. We close this chapter with Sec. 2.6 in which we first infer an 
upper limit on primordial ^Be and then present the results of our CBBN calculation 
which heavily constrains the X~-abundance/lifetime parameter space. 

2.1 Basic bound state properties 

The presence of negatively charged massive particles X~ during/after primordial nucleo- 
synthesis leads to the formation of bound states (NX^) with the ionized nuclei N of 
the light elements. In this section we shall describe the basic properties of such bound 
states. 
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In order to obtain a first estimate on the physical properties one can immediately 
apply the standard formute for the quantum mechanical motion in a Coulomb field. The 
characteristic size of the (NX~) system is given by the Bohr radius of the system, 

1 29 fm , , 

Ob = ' , (2.1) 

where A and Z are the atomic mass and charge number of the nucleus, respectively. In 
the second relation we have used that the reduced mass 

"T-red = \ (2.2) 

m-N + 

is given to good accuracy by the mass of nucleus, m-red — 'n-N ("ix ^ "^n)) and that 
roughly ~ A where = 938 MeV [32] is the proton mass. The binding energies 
of a point-like nucleus orbiting X~ are given by the well-known formula 



= ::r° - (-25 kev) (2.3) 



—- -25 keV — 2- 



where n denotes the principal quantum number. 

For the case of a proton bound state ■* ~ 29 fm so that the system is a factor 
of mp/rrif, ~ 1800 smaller than a hydrogen atom. Nevertheless, the p-X~ distance is 
still large when compared to the rms charge radius of the proton, (r^)^^ — 0.88 fm [61]. 
The situation changes for heavier nuclei. For example, considering (^LiX^), one finds 
that o!^^^'^ ^ — 1-6 fm whereas the measured ^Li rms charge radius reads {f''^)^-^^ — 
2.54 fm [61]. Thus, we expect corrections to the naive Bohr-like binding energies (2.3) 
once the finite size of the nucleus is taken into account. 

In order to obtain more realistic values for the ground state energy, we need to 
make an assumption on the charge distribution of the nucleus. A compilation thereof is 
presented in [62]. We employ the Gaussian p = eZ'(^/7r)^/-^e~^^^ with radial coordinate 
r from which the potential 

0(r) = erf [y/lr] (2.4) 



47rr 

is obtained upon solution of Poisson's equation^ V'^cj) = p. Requiring (r^)p = (r^)N 
relates the parameter ^ to the rms charge radius; (r^)p = 3/2^. The error function is 
defined by erf (x) = 27r~^/^ e~^^dt. 

The above choice of the charge distribution is particularly convenient because the 
electric potential (2.4) is given in analytical form (2.4). This makes an application of the 



^We use Heaviside-Lorentz units with e — ^/Aiva. 
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Table 2.1: Basic quantities for some selected light elements and their bound states 
with X~ for mx — > oo. If necessary, nuclear masses ttin are derived^ from [44]. Rms 
charge radii {r'^)]^^ are taken from [61], ah denotes the Bohr radius (2.1), E^°q^ is 
the naive Coulomb ground state energy (2.3), and £^™'' provides realistic values of the 
binding energy from the variational principle (2.6). 



bound state 


rriN [MeV] 




ah [fm] 


Kf [keV] 


E^^' [keV] 




938 


0.88 


28.8 


-25 


-25 


{DX-) 


1876 


2.14 


UA 


-50 


-49 


(^HeX-) 


3727 


1.67 


3.6 


-397 


-347 




5601 


2.54 


1.6 


-1342 


-797 




8393 


2.52 


0.8 


-3575 


-1469 



^ Nuclear masses are obtained from atomic masses by subtracting Znic and correcting for the 
total binding energy of all electrons where we follow the prescription given in [44]. 



Rayleigh-Ritz variational method straightforward. Using the (unnormalized) trial wave 
function 

jPir- a, b) = e-""'/"^ (1 + br/ab) (2.5) 

with variational parameters a and b an upper bound on the true ground state energy Ey, 
can be obtained by minimizing the right hand side of 

For the Hamiltonian H of the (NX~) system we use H = — (2mN)~^V^ + ecj), i.e., we 
take mx oo. 

The results of minimization of (2.6) for selected light elements along with some other 
basic quantities are summarized in Table 2.1. Note that the Bohr radii of bound states 
with elements heavier than ^He lie within the nuclear radii. Thereby, the true binding 
energy for those systems is significantly reduced in magnitude as can be seen by com- 
paring E^°Q^ with E^'^'^ . We remark that the binding energy is an important quantity 
since it directly influences on the bound state fraction of the light nuclei. 



2.2 Wave functions of the relative motion 



For the calculation of photo-dissociation and recombination cross sections which include 
the finite charge radius correction, we are in need of the actual wave functions of the 
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N-X~ system. In the following we shall therefore obtain the wave functions for the 
(NX~) bound states as well as for the N~X~ continuum. It will also allow us to see 
how well our variationally obtained upper bounds E™ fit the actual value of the true 
ground-state energy E]^. 

The Schrodinger equation for the radial part R{r) = u{r)/r of the wave function 
tp{r,6,(j)) = Yi m{G,(p)R{f) of the relative motion is given by 



+ 



2m,ed(^ - V) 



(2.7) 



As usual, Y] 



denotes the spherical harmonic with orbital and magnetic quantum 



numbers / and m. For the potential V we make the following choices 



V 



-Za/r 
e (j){r) 

-Za/(2i?o) (3 - r^/Rl) 



point 
gauss 

h.sph (r < i?o) 



(2i 



where "point" stands for the Coulomb potential of a point-like nucleus, "gauss" for a 
Gaussian charge distribution with defined in (2.4), and "h.sph" for a potential of a 
homogeneously charged sphere of squared radius R^ = 5(r^)Ar/3 [62]. For r > Rq, Vh.sph 
is to be continued by ^oint- 



2.2.1 Discrete spectrum 

We solve (2.7) for E <Q and the various choices of V [Eq. (2.8)] numerically. For fixed 
n = nr -|- Z -|- 1 we exploit the fact that the radial function R[r) and thus u{r) vanishes rir 
times; rir is the radial quantum number. The solution of (2.7) is fixed by imposing the 
standard boundary conditions u{5) = (^'+^ and u'{5) = (/ + 1)5^ with 5 <^ {'''c)]^^ and 
normalizing to unity, f \u\'^dr = 1. 

At the top of Fig. 2.0 we show the numerical solutions of the normalized radial wave 
function Rni = Rio for the (^HeX~) ground state (IS in the usual spectral notation) 
for the different choices (2.8) of the potential. An attenuation of the wave functions 
with finite charge radius relative to the Coulomb case can be seen at small r. At large 
radii, the wave functions are Coulomb-like.^ It can further be seen that i?gauss — -Rh.sph 
for all r, i.e., the radial wave function for (^IIeX~) is rather insensitive to the concrete 
choice of the charge distribution. In the middle and at the bottom of Fig. 2.0 we plot 
r^R?, i.e., the probability density of the ^He-X~ distance, for n < 3 and Gaussian 
charge distribution. Except for small radii, the curves essentially resemble distributions 

^Of course, the Coulomb solution is simply given by i?point = a^^^^ exp (— r/ob) 
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I 
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I 




Figure 2.1: Top figure: Radial wave functions R of the (*HeX^) ground state for the 
choices (2.8) of the potential as labeled. Lower figures: probability densities r^R^ of 
the ^He-X~ distance for eigenstates with n < 3 and Gaussian charge distribution. All 
curves are obtained for mx oo. 



obtained with Coulomb wave functions. This is particularly true for the higher I states 
as the wave functions are pushed outwards due to the centrifugal term in (2.7). 
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Table 2.2: Complete spectrum for (''HcX") for n < 3. Binding energies are given for 
and mx oo and mx = 100 GeV (bracketed values) for the potentials as labeled and 
given in (2.8). Additionally, the expectation values (r) and rms radius {r^y^'^ for the 
"gauss" case are provided. 

("Hex-) {r^yJl = 1.67 fm 

mx ^ oo (100 GeV), Ob = 3.63 (3.76) fm 



State 


£;P°'°*[keV] 


^gauss jj^^^j 






IS 


-397 (-383) 


-348 (-338) 


1.7 


2.0 


2S 


-99 (-96) 


-93 (-90) 


6.4 


6.9 


2P 


-99 (-96) 


-99 (-96) 


5.0 


5.5 


3S 


-44 (-43) 


-42 (-41) 


14.1 


15.0 


3P 


-44 (-43) 


-44 (-43) 


12.5 


13.4 


3D 


-44 (-43) 


-44 (-43) 


10.5 


11.2 



In Table 2.1 the spectrum for the cases "point" and "gauss" [Eq. (2.8)] is given for 
the states plotted in Fig. 2.0. In addition, also the expectation value (r) as well as 
the rms radius (r^)-*^/^ are given for the "gauss" case in units of Cb = 3.63 fm. We 
solve the Schrodinger equation (2.7) for mx — > oo, i.e., for m-red = "^N, as well as for 
mx = 100 GeV (bracketed values) in order to study the influence of a finite X~ mass 
on the binding energies. The table shows that for mx = 100 GeV this leads to a shift of 
10 keV for the ground state energy but the correction quickly becomes marginal for the 
n > 1 states. The same is true when comparing the spectra for the different potentials. 
Whereas the correction to the ground state energy is substantial, 49 (45) keV, the higher 
states for the (^HeX~) system essentially coincide. It is, however, interesting to note 
that the Coulomb degeneracy is broken. We refrain from showing the energies for the 
case "h.sph" since they are the same as for the "gauss" case (except for n = 1 where a 
1 keV shift is found.) One can also see that the variational determination of the ground 
state energy in section 2.1 gave an accurate result. 

For bound states of X~ with heavier nuclei than ^He, i.e., for more compact systems, 
we expect a pronounced behaviour of the observed effects above. Analogously to the 
case (^HeX~) we can analyze (^BeX~). This is an interesting system because free ^Be 
is unstable by 92 keV and decays into two alpha particles: ^Be ^He + ^He. Indeed, 
the stable (^BeX~) system is part of a CBBN reaction chain which can open the path 
to primordial production of ^Be [63]; see Sec. 2.4. Since the lifetime of *Be is ~ 10~^^ s 
no experimental data on the charge radius of the isotope is available. In this section 
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Table 2.3: As in Table 2.1 but for (®BeX ). In addition binding energies for the 
"h.sph" case are shown; see (2.8). 

i'BeX-) {r^,)'Jl = 3.39 fm 

mx ->oo (100 GeV), ab = 0.91 (0.97) fm 



State 




keV] 


ph.sph 


[keV] 


^gauss jj^^^j 


{r) [ab] 




IS 


-3176 ( 


-2956) 


-1118 (- 


-1092) 


-1168 (-1138) 


3.8 


4.2 


2S 


-794 


(-739) 


-458 


;-437) 


-475 (-453) 


9.8 


10.6 


2P 


-794 


(-739) 


-652 


-620) 


-650 (-618) 


6.4 


6.9 


3S 


-353 


(-328) 


-243 


-230) 


-249 (-236) 


19.0 


20.2 


3P 


-353 


(-328) 


-306 


;-290) 


-307 (-290) 


14.5 


15.5 


3D 


-353 


(-328) 


-348 


;-325) 


-346 (-323) 


10.9 


11.6 



we follow [64] and adopt the value (?'c)8Be ~ '^ '^^ ^™ which is based on a microscopic 
^He + ^He model calculation [65]. Again, in Fig. 2.1 we plot the IS radial solutions of 
the Schrodinger equation (2.7) for the various potentials (2.8). The difference between 
-Rpoint and -Rgauss (-^h .sph) is now substantial. Moreover, also a slight difference between 
R 

gauss and -Rh.sph is observable for smaller radii. We therefore expect a dependence of 
the ground state energy on the adopted charge distribution. 

In Table 2.2 we provide the complete spectrum for (^BeX~) with n < 3. Expectation 
values (r) as well as the rms radius (r"^)^^"^ are also computed for the "gauss" case in 
units of Ob = 0.91 fm; msBe = 7.455 GeV. Again, we compare the energy eigenvalues 
for mx oo with the ones for mx = 100 GeV (bracketed values). As can be seen, all 
states now receive substantial corrections to the Coulomb values. Moreover, we observe 
a 50 keV (46 keV) shift in the IS energy when changing the charge distribution from 
Gaussian to square well (in r). 

Finally, we have checked all variationally determined ground state binding energies 
presented in Table 2.1 of the last section by explicit computation of the wave function. 
We find that all E^^^ given in Table 2.1 are within 1 keV of the numerically obtained 
result. Noteworthy may be the 1 keV shift for ('^HeX~). Of course, not only the 
assumed distribution of charge influences on E^, but also the error on the measured or 
theoretically predicted charge radius is a source of uncertainty. This is of pronounced 
importance for the heavier nuclei because the bound state system is more compact. For 
our purposes, however, it is not essential to pursue this issue further; see Sec. 2.4.2 
for another comment in the context of catalyzed ^Be production. In the following, we 
employ the binding energies determined from the Gaussian charge distribution. 
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2.2: IS radial wave functions for (^BeX~) from solving (2.7) with potentials 
labeled and mx oo. 



2.2.2 Continuous spectrum 

We are also in need of solutions of the Schrodinger equation (2.7) for ii^ > if we want to 
obtain charge-radius corrected bound-state formation cross sections. The normalization 
of a numerically obtained solution is more involved since the wave functions of the 
N-X~ -continuum are not bounded spatially. However, a finite charge radius leads to 
a modification of the Coulomb form of the potential only in the vicinity of the origin. 
Therefore, we can take the following approach: For r S> (^'c)?/^ solution of (2.7) has 
to be a linear combination of the regular and irregular Coulomb wave functions F^i and 
Gki, respectively. They can be expressed as 

FM = ^iYi+Yn , (2.9a) 

GM = ^^iYr-Yi) , (2.9b) 

with [66, 67] 

Yi = +^ r|!t! + i!!)' '^^'^'''^^' W,,j^y,{mr) , (2.10a) 

^* = -^WTT^""^'^'^'"'^^' W.,,J^y,{-2^kr) . (2.10b) 

Here, rj = — l/(A;ab) denotes the Sommerfeld parameter for an attractive Coulomb 
field where k is the wave vector of the relative N-X~ motion with |k| = k = {Irrij-cdE)^^'^; 
r(z) is the Gamma function [68] and Wa^i,{z) stands for Whittaker's function [69]. 

With W±a,b{^z) = eT^/2(±z)±'=' [1 + (z-i)] [69] one finds that the asymptotic 
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behavior of the wave functions (2.9) is given by 

hm Fki = +sm [kr -'q\n{2kr) -ttI/2 + Si], (2.11a) 

fcr— >oo 

lim = -cos [A;r-r/ln(2A;r) -7r//2 + <5,], (2.11b) 

fcr— >oo 

where the Coulomb phase is defined by 5i = arg r(/ + 1 + i-q). 

Now, for r ^ (^c)n^ radial solution of the Schrodinger equation to a modi- 
fied Coulomb potential can be written as^ = (2/r) [aiFki — biGki]- Retaining the 

■^This definition corresponds to normalization on the "fc/27r scale", RkuRkir^dr — 2-K5{k' — k). 
Note, however, that Gki is not regular at the origin. One has to introduce a cutoff factor if 7?^; shall be 
an entire function; see [70]. 
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asymptotic normalization 



lim R 



lOUt 



sin [A;r — r/ In (2A;r ) — 7r//2 + 5^ + o";] 



(2.12) 



it follows from comparison with (2.11) that the additional phase shift (Ji{k) is given by 
ai = tan{bi/ai). 

In the vicinity of the origin, i.e., for r < (r^))/^, the numerically obtained wave func- 
tion itl^* correctly describes the solution to Schrodinger's equation. It can be normalized 
by requiring a continuous transition at r = rfit to the outer solution 



Following the outlined approach, we solve (2.7) for the relative motion of various 
N-X~ systems. As examples, we choose ^He-X- and ^Li-X' . At the times of BBN 
the relative velocity of the ^He/^Li-X" system is Boltzmann distributed so that (k) ~ 
V nT'vedT ■ Thus, a representative value is /c = 10 MeV which corresponds to T ~ 
20 (30) keV for ^He (^Li) with mx oo. We join the inner solution with the outer one 
at rfit = 10 fm. This determines the phase shift u/. We have checked that ai is insensitive 
to the chosen value of rgt, provided r > (r^)J^^^, and that o"; — > when switching to a 
point-like nucleus. Using a Gaussian charge distribution, we find for the S-wave {I = 0) 
of the ^He(^Li)-X^ system ao = —0.21 (—0.87) whereas for the P-wave {I = 1) the phase 
shift is already significantly reduced, ai = —3.3 x 10^^ (—0.13). 

In Fig. 2.2 the corresponding wave functions for ^He-X" (top) and ^hi-X~ (bottom) 
are shown. As can be seen, the wave functions Rko for the case "gauss" receive a signifi- 
cant correction in comparison to the Coulomb case "point" which was already indicated 
by the size of the phase shifts (Tq. Of course, the curves labeled "point" coincide with 
the regular Coulomb functions {2/r)Fki [Eq. (2.9a)]. Whereas for the ^He-X~ system 
Rf^i ~ (2/r)Ffci, a deviation from the regular Coulomb P-wave is visible in the ^Li-X~ 
case. 

When considering continuum wave functions for — > the numerical evaluation of 
(2.10) is problematic. This case, however, is the most important one in the computation 
of the photo-dissociation cross section of (NX~) [Sec. 2.3.1]. Therefore, we need to 
consider the Coulomb wave functions in a different form [71], 




(2.13b) 



(2.13a) 




(2.14b) 



(2.14a) 
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where A {i/ka^, I) = 111=1 [l + (■s^flb)^] and yi,3 are related to the Whittaker functions; 
A{x,0) = 1. Defined in this way, (2.14) satisfy the asymptotic behavior (2.11). An 
expansion in powers of i.e. in energy, for yi reads [71] (see also [72]) 



3q 



q=0 p=2q 



(2.15a) 



whereas cannot be represented by a convergent expansion in powers of energy. How- 
ever, an asymptotic expansion has been obtained in [71], 

Q 3q 

y,{K, I; p) = A{k, I) i^-'" E a,,p(/)(2p)(P+i)/2y2,+,+p {V^p) + O (^-^Q-^). 

q=0 p=2q 

(2.15b) 

Here, J21+1+P and ^2/+!+^ are the respective Bessel functions of the first and the second 
kind of order 2/ + 1 +p [68] and the coefficients aq,p{l) satisfy recurrence relations (ao,o = 
1); for details see [71]. With the expansion 

A {i/ka^,l) = 1 + ^l±^M±^eal + O {k'at) (2.16) 
one obtains from (2.15) 



Fki =^/wrk I J21+1 ( \/8r/ab ) + 



12 



/(/ + 1)(2/ + 1) J2/+1 (V8r/ab 
3/2 



-3(/ + 1) (^^^ J21+3 (^8r/ab) + (^^) J21+4 (\/8r/ab) + O (k^at) 



(2.17a) 



Gki =V^ < Y21+1 ^/8r/ab) + 



12 



/(/ + 1)(2/ + I)y2m (VSr/ob 



-3(/ + 1) (^^^ Y21+3 (78^7^) + (^^^ ^ Y21+4 (a/s^T^) + O {k^at) I 

(2.17b) 

In the following section we employ those approximations in the computation of the 
photo-dissociation cross section. 



2.3 Formation of bound states 



The crucial quantity in the discussion of the catalysis of BBN reactions is the bound 
state fraction n(NX-)/^N of the light elements N. To this end we have to compute the 
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cross sections dree for radiative recombination, N + X — > (NX ) + 7, as well as cTph for 
the dissociation, (NX~) + 7bg — > N + X~ due to background photons 7bg. 

The rate (per nucleus N) for N-X~ recombination is given by Free = {o'mcv) nx- ■ 
Here, cirec has to be averaged over the distribution of relative velocities v between N and 
X~ which is Maxwellian for a sub-MeV plasma, 



Hence, 



{(TrccV) = J^^T-^I^ r dE^ CTrec^ke-^^/^, (2.19) 

V vrmred Jo 



where Ey^ = my-^^v^ /2 has been used. 

The rate of photo-dissociation Fph of (NX^ ) pairs depends on the number of photons 
whose energy E^ exceed that of the ionization potential |-Eb| of the bound state, 

1 /■°° E"^ 
n^{E^ > l^bl) = — / dE^ ^ , (2.20) 

and is given by Fph = ap\in^[E^ > \Ei,\). 

The principle of detailed balance [34] relates the cross sections via p^fTrec = 2£'^o"ph 
where |p| = my-g^v denotes the momentum of the relative motion of the N-X~ system 
and the factor of two is a statistical factor accounting for the two polarization degrees 
of freedom of the photon. From the definition of the rates Fph and Ffoc together with 
(2.19) and (2.20) it follows that"^ 

^=(^Y'e\^^\/^nx-. (2.21) 
Fph \m^cdT J 

As long as Fph(r), T^f.f.{T) ^ H(T), i.e., as long as recombination and break-up reac- 
tions happen frequently, the concentrations of N, X" , and (NX~) have time to achieve 
equilibrium values such that the reaction densities for recombination and dissociation 
are equal, Fphn(Nx-) = Frec^^N- This yields the Saha equation for the bound state 
fraction,^ 

"-(NX-) _ Free 



Fph 



2.22) 



''It is used that = E]^ + \Eh\ and that [exp (E^/T) - 1]"^ ~ exp {-E^/T) which holds well in the 
temperature regions of main interest. 

^When used in this form one may need to impose that the number of bound states cannot be larger 
than the total number recombination partners available. 
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2.3.1 Photo-dissociation and recombination cross section 

Since the early Universe is in a high-entropy state, bound states can only form efficiently 
once T <\E\^\/4Q [73]. At the relevant times, i.e., when Fph < H, only those photons in 
the high energy tail of the spectrum with > \E]~\ are capable of destroying (NX~); 
{E^) ~ 3T. In addition, the binding energy is significantly smaller than the associated 
light element mass so that a non-relativistic treatment of the photoelectric effect is 
perfectly justified. 

For the computation of the photo-dissociation cross section we can employ Fermi's 
golden rule. The probability per unit time for a nucleus bound to X~ to undergo a 
transition into the continuum is given by 

du;ph = 2Ti\Vfi\^5{-\E^,\ +E^- Ey,)dp . (2.23) 

After the transition the nucleus has kinetic energy Ey^ and momentum p.^ The density 
of final states is dp = (i^p/(27r)^ and the matrix element for absorption of a photon with 
momentum k and energy E'^ = k^ reads 

V/^ = ^e,j^^,(-k) (2.24) 

where Jjj(k) = f d^x e~'^^'^j'^'{x) is the Fourier transform of the transition current 
jj^{x) = ij)j^^ipi and denotes the photon polarization vector; see, e.g., [48]. The 
cross section is found by dividing (2.23) by the incident photon flux density. Averag- 
ing over photon polarizations [in the gauge {e^) = (0,e)], and integrating over E\^, the 
differential cross section for photo-dissociation is given by 



dQ.p AnEy 



ekXj/i|^ (2.25) 



where is a unit vector in k-direction and jjj is the spatial part of jj ■{—]<) . 

We shall consider ionization from IS as well as from 2S states. The initial state wave 
function is V'i = (47r)^"'^/^i?xg/2s- The final state has to comprise a plane wave in p 
direction together with an ingoing spherical wave [74]. In the partial wave expansion, 

^ oo 

^/ = ^ E ^'(2/ + 1) e-('^'+-')i?|p|,P,(ep . er) (2.26) 
'P' 1=0 

with unit vectors ex in x-direction; Pi are the Legendre polynomials [68]. Note the ap- 
pearance of the additional phase shift ai coming from the finite charge radius correction. 



Of course, strictly speaking, it is the energy and momentum of the relative motion. From the above 
explanations, however, it is clear that the X~ recoil is negligible. 
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Table 2.4: Listed below are the cross sections Cph for the threshold {E^ — \E\^\) 
bound-free transition from IS and 2S states. From the definition (2.30) the averaged 
cross sections (arecv) for recombination into IS for (p^~) and into 1S-I-2S for {*lleX~) 
and {^LiX~ ) are obtained in the third column. Bracketed quantities refer to the "point" 
case. In addition, a critical temperature of bound state formation Tj-oc defined by 
rph(7rcc) = HiT^-cc) is given in the last column. 



bound state 








T 

rec 




[mb] 


[mb] 


[cm^s^^moP"'^] 


[keV] 


(^HeX-) 
(6LiX-) 


1870 
118 
34 (52) 


4380 
294 (278) 
103 (123) 


3980 (IS) 

7260 (9230) 

6640 (25370) 


0.6 
8.3 
19.0 



Since the wavelength of the ionizing radiation (A = l/l-Ebl on the threshold) is much 
larger than the (NX^) dimensions, we can use the electric dipole approximation. The 
associated selection rule implies / = 1 for the continuum so that 

37>-«{5i+o"i) 

^f = ^ (ep-er)i?|p|i , (2.27) 

^IPI 



and thus (in the dipole approximation 

3i 



J/i 



yievrmredl 



- j d\ (ep • er)-R|p|iVi?is/2s • (2.28) 



Performing all angular integrations in (2.25) yields for the total photo-dissociation 
cross section 

2 

(2.29) 



fph 



•irrircdE-y 



1 f°° d 
7\^\l *^'^lp|ia;:^is/2s 



For -Rjp|i and -Ris/2S we employ our numerically obtained solutions of the previous 
section which takes into account the finite charge radius of the nucleus. Note that 
on the ionization threshold (Tph is independent of |p|. For a pure Coulomb field the 
momentum dependence cancels analytically when using the leading term in (2.17a). By 
the same token, numerically, [. . . ] in (2.29) becomes insensitive to |p|. Thus, we find a 
constant cross section for \E]j\. In this limit, using detailed balance, the averaged 

recombination cross section reads 
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from which Fph is readily obtained by using (2.21) 

T,^ = {a,,cv)(^^Y e-\^-\l^. (2.31) 

In Table 2.3 we present the results on the threshold photo-dissociation cross sections 
CTph for transistions from IS and 2S states for the elements p, ^He, and ^Li and for 
mx — *■ oo. The bracketed values are for a pure Coulomb potential whereas the other 
results are obtained by using a Gaussian charge distribution. The respective values do 
not differ very much. This is because the decrease of |i?b|(= E^) in the denominator 
of (2.29) when switching from the "point" to the "gauss" case is counterbalanced by an 
increase in the radial integral so that the net effect is small. However, the reduction 
of the total (1S+2S) recombination cross section (arecv) from the hydrogen-like case is 
drastic. This is due to the additional factor of E'^ in (2.30). In the last column we 
show the temperature for which rph(Ti.ec) = H(Tj-ec), i.e., the temperature when the 
formation of bound-states can proceed efficiently — provided that Free ^ H and that the 
bound state is not destructed by another process. 

Finally, we remark that for other (heavier) nuclei than the ones presented in Table 2.3 
the discussion of recombination can become more involved. If the light element N pos- 
sesses an excited state A^* with a level splitting smaller than the X~ binding energy, 
then recombination may also proceed into (N*X~) opening up the possibility of reso- 
nant recombination. This was pointed out in [75] where the formation of (^BeX^) was 
considered. 

2.4 Nuclear reactions with bound states and their catalysis 

After the freeze-out of weak interactions with the cease of n and p interconversion pro- 
cesses, light element fusion in SBBN proceeds via inelastic two-body nuclear reactions^ 

B + C ^ D + E I 
B + C ^ F + -f II 

with B, . . . ,F denoting the nuclei of the light elements and the arrow indicating the 
forward process, i.e., the exoergic direction with positive Q value. The reverse processes 
are typically suppressed by exp(— Q/T) such as in (2.21) (which is an atomic process.) 
Only elements with atomic mass number A < 7 are produced in relevant quantities. 

^This does not include "production" processes like that of ^Li via electron capture by ^Be or of ^He 
by beta decay of T, both of which, however, only happen at a much later time. 
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In presence of bound states of the hght elements with X during BBN the following 
additional types of inelastic reactions emerge as particularly prominent, 

(Bx-) + c ^ D + E + X- r 
iBx-) + c ^{Fx-) + -/ ir 

{BX-) + C ^ F + X- III 
{BX~) + C ^ B + {CX-) IV. 

A first observation is that, in presence of bound states, the energy gain of a nuclear 
reaction is altered. In the entrance channel, the total available internal energy is reduced 
due to the binding of B with X~ . Thus, for example, Qi* = Qi — \E^ \ whereas 
additional energy becomes available in the exit channel of II* so that Qn* = Qu — 

Since Q values of nuclear reactions are mainly in the MeV to 
multi-MeV range, usually three-body break-up reactions I* rather than {DX~) + E exit 
channels are realized. The shift in energetics can also allow for resonances which are not 
possible in SBBN. For example, type II* can be realized in resonant capture reactions 
whose intermediate excited state (FX^)* decays into the {FX~) ground state by 7 
emission. If, instead, the nucleus is in an excited state {F*X~), then also the F + X~ 
continuum acts as a concurrent channel — provided that it is kinematically accessible. 
The latter is an example of a reaction of type III which is of particular interest since it 
has no SBBN counterpart. Atomic reactions IV are called charge exchange reactions. 
They are also important to consider because they can significantly affect the relative 
concentrations of nuclei bound to X~ . They will be discussed in Sec. 2.5. 

Reactions of the form {BX^) + (CX'^) — > . . . are only of secondary importance. Their 
efficiency depends on the average relative velocity between {BX~) and {CX~) which 

— 1/2 

scales as m-^ . Thus, for weak scale relics, the suppression of the average velocity of 
^"-containing bound states relative to the velocity of light nuclei is from one to two 
orders of magnitude. 

Another observation is that the screening of the charge of B when in bound state 
with X~ will lead to a modification of the SBBN cross sections with charged "projectiles" 
C. It is customary to write the cross sections of charged-particle induced reactions in 
the form 

a{E^) = S{E^)E-'e-^-^ (2.34) 

and which defines the astrophysical S"- factor. The definition scales out the "geometrical" 
cross section vrA^ oc E^^ as well as the Coulomb penetration factor exp (— 27rr/). Note that 
during BBN (T < 0.1 MeV) the thermal energy (-Ek) of the reactants is significantly 
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smaller than the height of the Coulomb barrier — O (MeV); A = {m^eAv)~^ is the 
de Broglie wavelength of the relative motion and r] = ZsZcot/v > denotes the earlier 
encountered Sommerfeld parameter [below (2.10)]. Using (2.34) the definition of the 
thermally averaged cross section (2.19) becomes 



V vrmred Jo 



Ey 



T \E^) 



1/2 



(2.35) 



where Eq = 2mred('^CK-^B-^c)^ is called the Gamow energy. In absence of resonances 
the 5-factor is only a slowly varying function of £'k so that the integral is dominated by 
the exponential which peaks at Eq = Eq^ and which marks the energy range 
of most effective nucleosynthesis (Gamow window). 

One may then attempt to account for the bound state in the entrance channel by 
replacing Zb by Zb — 1 and correcting for the changed kinematics and energetics. In- 
deed, such a program has first been carried out in a BBN network calculation in [73]. 
When studying the effect on the charged particle induced reactions, the authors find no 
significant changes in the light element yields at the CMB inferred baryon asymmetry. 
Whereas the compactness of the bound states with the heavier of the light elements gives 
some justification to this procedure we will see in Sec. 2.5 that, e.g., the large size of the 
(pX~) system plays a crucial role in obtaining a consistent picture of BBN. 



2.4.1 Catalysis of ^Li production 

The potential influence of bound states on the BBN paradigm was already discussed 
almost twenty years ago in [76, 77, 78]. However, only recently it has been realized [79] 
that the presence of X~ at T < 10 keV can lead to a tremendous enhancement of the 
^Li output. 

In SBBN the freeze-out of ^Li from nucleosynthesis is dominated by its production 
via radiative capture and its destruction via proton burning, 

^He + D ^ "^Li + 7 Q = 1.47 MeV , (2.36a) 
'^Li + p ^ ^He + ^He Q = 4.02 MeV , (2.36b) 

respectively. The cross section for the production reaction (2.36a) is very small with 
S ~ 10~^ MeVb [80] at the lowest energies. For example, the 5-factor for ^He + 
'^He — > ''Be + 7 which is the main source for ^Li at 7/b(CMB) reads S{0) = 5.8 x 
10"^ MeVb [81]. The small SBBN output of ^Li is attributed to the inefficiency of 
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^d"ec = 4 X 10-4, TX = 105 S 




T [keV] 

Figure 2.4: Choosing X'^^^ ^ Ax IQ-^ and tx = 10^ s the yields of D, "^He, and 
^Li normalized to the proton number density rip are shown. For T > 10 keV the light 
elements follow their SBBN evolution (the output thereof is produced using the BBN 
code [83].) At lower temperatures bound states ("'HcX^) can form and catalyzed ^Li 
production proceeds via (2.37). The dashed lines show the CHeX") and ^Li abun- 
dance when the Saha approximation (2.22) is used. The resulting overestimation of ^Li 
production illustrates that a full numerical solution of (2.39) is necessary. The CBBN 
variation of D and ^He is negligible. 

the production process (2.36a).^ In this regard, also note that the destruction reaction 
(2.36b) has a large 5-factor, S'(O) = 2.97 MeVb [80]. 

Let us briefly outline the evolution of the ^Li abundance in SBBN. Once the tempera- 
ture of the primordial plasma drops below T < 0.1 MeV the deuterium bottleneck opens 
so that ^Li production can proceed via (2.36a). A sharp drop in D below T < 80 keV is 
accompanied by an associated decline in ^Li; see left part of Fig. 2.3. Net production of 
^Li via (2.36a) soon freezes out but proton burning (2.36b) continues until ~ 10 keV. 

In the previous section we have seen that once the temperature drops below 8 keV the 
photo-dissociation rate of (^He^-) freezes out. Thus, provided that tx is large enough, 
the concentration of (^HeX-) can become substantial and fusion of ^Li is then possible 

*This is usually traced back to a weak quadrupole transition (E2) in (2.36a). Theoretical calculations 
seem to suggest, however, that at the BBN relevant temperatures (T ~ 70 keV) the dipole transition 
(El) is as important [82]. Nevertheless, the dipole moment of the ''He-D cluster is almost vanishing due 
to a similar charge-to- mass ratio which yields a very small cross section; cf. [79]. 
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via the alternative path [79] 

(^HeX") +B^%i + X- Q ~ 1.13 MeV (2.37) 

and which is a reaction of type III. Whereas the size of the radiative capture cross 
section (2.36a) is governed by the selection rules of the electromagnetic transition, (2.37) 
suggests a cross section which is determined by the short distance behavior of the ^He-D 
cluster. Indeed, the original work [79] estimates 5cbbn / S'sbbn ~ 10^ which points to a 
cross section for (2.37) which is in the ballpark of photonless SBBN reaction rates and 
implies the catalysis of ^Li production. 

Meanwhile, a dedicated quantum three-body calculation of the ^He + D + X~ system 
has become available [84] which confirms the catalytic picture. The authors find S{Eq) = 
0.038 MeVb at the Gamow peak position Eq = 36.4 keV. The thermally averaged cross 
section reads [84] 

NA{cT^^t,^uv) = 2.37 X 10« (1 - 0.34 Tg) T^^^^ exp (-5.33 Tg"^/^) (2.38) 

and is given in the customary units of cm'^s~^moP^; Tg denotes the temperature in units 
of 10^ K and A^a is the Avogadro constant [85]. 

Since (^HeX~) only forms for T < 10 keV, i.e., at a time when the ^He and D 
abundances are essentially frozen out, we can incorporate the effect of catalytic ^Li 
production in the following way. At some low temperature T < 20 keV we couple the 
SBBN output into the network of Boltzmann equations 

{(^Tcc,mcv) "^HeX" - rph,4He C^HeX") 
- (^cat,6Li«> D i^ReX-) - Tx (^HeX-) , (2.39a) 

-(c^rec.-lHe^') ^HeX" + rph,4He (^HeX") 

+ {<^c.t,^uv) D (^HeX-) - Tx (^HeX") , (2.39b) 

— (<7 rec.^He''^) "b He X~ + rph,4He (^HeX-) 

+ rx(^HeX-), (2.39c) 
(o-cat,6Li^^) nh D (^HeX~ ) - (fTdes.sLi^) "b P ^Li , (2.39d) 
-(^^cat,6Li^^)nbD(4HeX-) . (2.39e) 

Light elements as well as bound state abundance are normalized to the baryon num- 
ber, N = nN/^b, (^HeX~) = ?T.(4HcX-)/^b) and X^ = nx-/n\y — better overview shall 
compensate for the slight abuse of notation. 



-HT— (^HeX-) 



-HT-!^ X- 
dT 



-HT^^Re 
dT 



-HT—^'U 
dT 

-HT-^D 
dT 
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The central input parameter for the catalytic production of ^Li (and ^Be, see below) 
is the abundance of X~ at the time of its recombination with ^He. Above 10 keV, we can 
track the resulting (^HeX~) abundance by using the Saha-equation (2.22) since photo- 
dissociation and recombination proceeds rapidly. Only at T < 10 keV, C^HeX") starts 
to build up efficiently so that we couple it into the full set of Boltzmann equations (2.39). 
We parameterize X~ by its abundance prior to decay by introducing X^^^, where the 
superscript "dec" stands for decoupling, and by the X~ lifetime tx = T^^, so that 
the (total) X~ abundance at any moment during BBN is given by X^{t) = X^^^ x 
exp(— t/rx)- The SBBN cross section {cr^^gej^iv) for residual ^Li destruction (2.36b) can 
be found in [86]. We solve (2.39) using as initial conditions the SBBN output values of 
the computer code [83]: = in^Jn^ = 0.248, D/H = 2.6xl0-^ ^Li/H = 1.14x10-^^, 
and rip/nb = 0.75; furthermore, Qes = 3.36 and /icff = 3.91. 

In Fig. 2.3 we show the evolution of the D, ^He, and ^Li number densities normalized 
to Up for the exemplary choice X~^^^ = 4 x 10~^ and tx = 10^ s. Rapid photo-dissociation 
delays bound-state formation of (^HeX"^) until T 10 keV. Once (^HeX~) forms, 
catalyzed ^Li production proceeds efficiently via (2.37) which leads to the steep rise of 
^Li at T ~ 10 keV. As can be seen, the D and ^He reservoirs are essentially unaffected 
by this. In addition, the dashed lines show the (^HeX~) and ®Li abundances when — 
instead of solving (2.39a) — the Saha approximation (2.22) for the bound state fraction 
is used. This results in an overestimation of the ^Li output. The reason is that ^He-X~ 
recombination itself is efficient only for a short period after photo-dissociation freezes 
out. Thus, in order to obtain a reasonable estimate on ^Li|cBBN) a numerical solution 
of the Boltzmann equations (2.39) is necessary. 

2.4.2 Catalysis of ^Be production 

Another dramatic catalytic enhancement is seen in the production of ^Be. The yield 
of ^Be in SBBN is tiny: Whereas the short lifetime of ^Be [see Sec. 2.2.1] renders the 
neutron capture reaction ^Be + n — > ^Be + 7 inefficient, ^Be production via fusion on 
Li isotopes yields no more than ^Be/H < 10^^^ at ryb(CMB) [87]. The catalytic path to 
^Be is shown by the following sequence [63] 



which goes through the "double bottleneck" of (^HeX ) and {^BeX ). Bound states 
(^BeX~) are formed by the radiative fusion 



X- ^ i^ReX-) ^ i^BeX-) ^ ^Be 



(2.40) 



^He + i^BeX-) (^^BeX") + 7 



(2.41a) 
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and the catalysis of Be production is triggered by the photonless recoil reaction 



(^BeX~) + n ^ ^Be + X- . (2.41b) 
The respective cross sections for (2.41a) and (2.41b) have been obtained in [63] and read 

^A{crci,t,»Bev) = 10^ T"^/^ [ 0.95 exp (-1.02/r9) + 0.66 exp (-I.32/T9) ], (2.42a) 

A^A(cTcat,9Be^) =^ 2 X 10^, (2.42b) 



in units of cm^s^^mol^^. A comment is in order here. Recently, the catalytic path (2.40) 
has been questioned in [64]. On theoretical grounds it is argued that the charge radius 
of ^Be, (?'c)8Bg = 2.5 fm, adopted in [63] is too small; see Sec. 2.2.1. Since the neu- 
tron capture reaction (2.41b) proceeds resonantly via the excited state (^Bei/2+-'^~) — ^ 
^Be3/2- + X~ , the larger charge radius {rc)l^^^ = 3.39 fm as proposed in [64] decreases 
I and shifts the resonance O (100 keV) below threshold. A final answer on the 
efficiency of (2.41b) can only be obtained by a full quantum ^He+^He+n+X~ four-body 
calculation and is announced in [64] as work in progress. 

Given the unique sensitivity to physics beyond the Standard Model a primordial 
origin of ^Be offers, we choose to incorporate (2.40) into our reaction network. The 
following Boltzmann equations describe the production of ^Be, 

-HT-^ («BeX-) = (a,,t,8Bet^> 'He (^HeX") - Tx (^BeX") , (2.43a) 
-HT-^ 9Be = (a,,t,9Be^> '^b n (^BeX") , (2.43b) 
-HT-^ n = ^ (cTfus^sHe^') nb D D + {(Ji^^i^^v) nb D T - {o^^^s^^v) nb n^He 

- {^c.t,^Bev) nb n («BeX-) - r„ n . (2.43c) 

Again, as in (2.39), the abundances are normalized to nb and written in an obvious 
notation; n = n^/n\^. For T < 15 keV, the SBBN neutron abundance can already be 
tracked well by including the processes D D ^ n ^He, T D ^ n + ^He, and 
^He-|-n — > p-|-T into the reaction network [88]. The respective cross sections (cfus^He^)) 
((Tfus 4He'f^) , and (o"(jes 3He'^) ^re taken from [86]. Note that at the time of catalyzed ^Li 
fusion the D reservoir is essentially unaffected by those residual SBBN reactions so that 
we can neglect the back-reaction on (2.39e). The yields of T and ^He are taken from the 
output from an updated version of the Kawano code [89]. For the neutron lifetime we 
use Tn = = 885.7 s [32]. 
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T [keV] 

Figure 2.5: Evolution of catalyzed ^Li and ^Be production shown together with the 
formation of the "bottle-neck" abundances of ('*HeX~) and (^BeX~) for X'^^^ = 5 x 10~^ 
and Tx ~ 5 X 10'^ s. The dashed line gives the neutron abundance while the dotted line 
shows the abundance of free X~ . 



In addition, we supplement the right hand sides of (2.39a), (2.39b), and (2.39c) with 

^{4HeX-) = -(^cat,8Be^^) "Hc (^HeX"), (2.44a) 

A^l = {cT^^t^Bev) nb n (^BeX"), (2.44b) 

^^'i = ^%ex-y (2-44C) 

respectively. In our code we can neglect the formation of {^BeX~) that proceeds via 
molecular bound states (^HeX^) [63]. This process becomes important only for a com- 
bination of large Yx- and large tx, i.e., a parameter region which is already excluded 
by ^Li overproduction. Also note that at the time when (^BeX~) form, their photo- 

(^BeX~)\ 

dissociation is not important because of the high binding energy — ;| ^ 1170 keV; 
see Table 2.2. It is important to note that we assume the SBBN value for the deuterium 
abundance. The early decays of X~ may result in an injection of nucleons into the sys- 
tem. This typically drives the deuterium abundance upward, resulting in an enhanced 
number of neutrons at later times and therefore in an increased output of ^Be, with the 
general scaling ^Be ~ const x (D/Dsbbn)^- We choose to disregard this effect, noting its 
model-dependent character. We are allowed to do so since its inclusion can only make 
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the ^Be-derived bound on the X abundance stronger. 

Figure 2.4 shows the evolution of catalyzed ^Li and ^Be production from the solution 
of the corresponding set of Boltzmann equations below T = 10 keV. The curves in 
Fig. 2.4 are based on the values X~^^^ = 5 x 10^^ and rx = 5 x 10'^ s. When the 
"bottle-neck" abundances of (^HeX~) and (^BeX~) form, the catalytic paths (2.37) and 
(2.40) to ^Li and ^Be open up, resulting in the asymptotic values ^Be/H ~ 10~^^ and 
^Li/H ~ 3x 10~^^. The dashed line shows the neutron abundance and the dotted line the 
free X~ abundance, which is dominated by its exponential decay. We remark in passing 
that residual recombinations of ^He with X~ lead to the crossing of the (^HeX~) and 
X7 lines at late time. 



2.5 Charge exchange reactions and late time catalysis 

In this Section we discuss the role of bound states of X^ with protons. The (p^^) 
system may have a large impact on the BBN predictions because (i) the proton as a 
recombination partner is the most abundant element and (ii) it is a neutral system. As 
can be seen from Table 2.3, the recombination of (pX^) bound states becomes efficient 
only after the temperature drops below 1 keV which corresponds to a cosmic time of 
t > 10^ s. Thus, the question arises [76, 73, 90] whether a revival of fusion reactions, i.e., 
a late-time catalysis, can be triggered by the (potential) high reactivity of "neutron-like" 
bound states (pX~). 

As we will see in Sec. 2.6 the presence of even a modest number density of X~ during 
the recombination with helium can lead to a production of ^Li and ^Be at levels which 
are in stark conflict with their observationally inferred primordial values. In [90] it was 
claimed that large fractions of the previously synthesized ^Li at T ~ 8 keV can indeed 
be reprocessed by (p^~) via^ 

(pX-) + ^Li^^He + 3He + X- (2.45) 

Thus, it was advocated that allowed "islands" reconcilable with observations may well 
exist in the X~ abundance/lifetime parameter space for tx- ^ 10^ s. 

In the following we shall argue that this is not the case. One reason is that any 
arising (pX~) abundance is immediately intercepted by the very efficient charge exchange 
reaction [76] 

(pX") + ^He ^ (^HeX") + p. (2.46) 
^At the time of publication of [90] the catalysis of ^Be production [63] had not yet been realized. 
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Figure 2.6: Potential energy of the proton in the field of X~ at r = and an incoming 
nucleus at r — —Re- The potential energy is plotted along the line connecting X~ 
with ''He (solid line) or ^Li (dashed line), respectively. As the distance between the 
incoming nucleus and decreases, the potential well becomes more narrow, and 
the proton ground state energy level is pushed upward. The critical deconfinement 
distance Rj^'^' is defined as the distance at which the energy of the bound state 
found variationally becomes larger than the height of the barrier Vmax to the right 
oiX-. 



This reaction may have a very large rate as its cross section is determined by the actual 
size of the (pX^) bound state that is of the order of Cb — 30 fm ((r^)-*^/^ = 50 fm). 

The charge exchange (2.46) can best be understood by employing a semi-classical 
picture. Calling R the separation between '^He and X~ (or, more generally, the sepa- 
ration between X~ and the incoming nucleus of charge Z), the one-dimensional slice of 
the proton potential energy in the field of X~ and ^He is given by, 

^W = — +1 ^' (2.47) 

r |r — K| 

and is plotted in Fig. 2.5. The limit oi R ^ oo corresponds to an unperturbed binding 
of the proton to X~ with a binding energy of = —25 keV. For Z > 1 and finite R, 
the curve has a maximum at positive values of r referred to as T^ax- As the ^He nucleus 
comes closer, R decreases. At some point, the binding energy of the proton becomes 
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Table 2.5: Deconfining distances Rc and charge exchange reaction cross sections on 
the (pX^) target for incoming nuclei with different charges Z . 





Z = 1 


Z=2 Z=3 Z=4 


Rc [fm] 


40 


95 135 160 


a = -kRI [b] 


51 


280 580 850 



positive so that the tunneling of the proton to r ^ +oo starts to become viable. For 
even smaller values of i2, one can find the distance -Rc at which the probability for the 
deconfinement of the proton approaches unity due to the fly- by of the ^He nucleus. In 
order to estimate Rc, we employ the variational calculation of the proton energy in the 
potential (2.47) by using the (unnormalized) trial wave function for the ground state, 

ilj{pi, v) = exp[-(;u - v)R/{2aax,)] (1 + uR/ba^f, (2.48) 

where fx and u are elliptic coordinates and a and b the minimization parameters. The 
coordinates are defined as fi = {ri + r2) / R and u = (ri — r2)/R, where ri and r2 are 
the proton-nucleus and proton-X~ distances, respectively. We calculate the energy of 
the ground state E'™'' as a function of the distance R. This yields the critical separation 
Rc, i.e., the distance at which E^^'^(Rc) = Vmax, and which describes the situation when 
even a metastable bound state simply cannot exist. 

The cross section for the charge exchange reaction may then be approximated by 
the geometric one with the impact parameter p = Rc, a = ttR^. The deconfining 
distances Rc together with the estimated cross sections for charge exchange are presented 
in Table 2.4. The associated thermally averaged cross sections are given by {(Tcxv) = 
(Tcx(w) = o"ex Y^8T/ ( Trmred ) and are listed in the Table 2.5 found at the end of this 
Chapter. As can be seen from Table 2.4, a ^He-X~ distance of ~ 95 fm is sufficient to 
release the proton from the bound state. Consequently, the estimate points to a very 
large cross section of almost 300 bn for the charge exchange reaction (2.46). 

For the charge exchange on '^He [Eq. (2.46)] an exact solution of the three-body Schro- 
dinger equation has recently become available [64]. The authors confirm the efficiency 
of the charge exchange and find 

NA{(Tex,mev) = 1-0 X 10^° cm^s^VoP^ (2.49) 

This compares well with the estimate ~ 4 x 10^ cm^s~^mol~^ obtained from the semi- 
classical picture when evaluated at the fiducial temperature of 1 keV.^^ Note that the 

^°The diflFerent scaling T"^/^ of the averaged cross sections obtained in the semi-classical approach 



52 Chapter 2. Bound states and catalysis of BBN 



charge exchange will mainly proceed into the n = 3 level — preferentially with highest I 
(largest number of states); cf. Table 2.1. Capture into highly excited orbits is indeed 
observed in charge exchange reactions of muons on hydrogen which gives justification to 
the employed semi-classical approach. Reference [90] finds for the charge exchange on 
^He a rate (per particle pair) of ~ 2 x 10^ cm^s~^mol~^ which is an underestimation by 
two to three orders of magnitude. 

We incorporate the charge exchange reaction on ^He in the network of Boltzmann 
equations by solving 

--H'r^ (P-'^') = (o-rcc.p-u) nb pX" - Tph.p ipX~) 

- {<^e.Me^) nb ^He (pX-) - Fx (pX") (2.50) 

for the (pX~) abundance. The back-reaction on free protons is negligible so that we 
refrain here from writing an equation for p. However, we need to supplement the right 
hand sides of (2.39a), (2.39b), and (2.39c) by 

^\i}ieX-) = (^ex,mev)nb ReipX'), (2.51a) 
A^l = -(arec,pi')nbpX- +rph,p(pX-), (2.51b) 

^Ue = ~'^(4HeX-)' (2.51c) 

respectively. Using the recent result (2.49) for the charge exchange cross section, we 
find that in the limit of infinite lifetimes, tx- oo, the abundance of {pX~) reaches 
its peak at around T = 0.7 keV. Its maximum abundance at these temperatures can be 
well approximated as 

^-1.7X10- (A,), (2.52) 

where the assumption X~ < ^He has been made [see Sec. 1.3] ensuring the linear scaling 
in (2.52).^^ That such a small (pX^) fraction only has a marginal impact on ^Li can be 
seen by comparing the destruction rate for (2.45) [64] 

^A(o-cat,des,6Li^^) = 1-6 X 10® cm^s^^moP^ (2.53) 

to the Hubble rate at the relevant temperature of T = 0.7 keV. Using X~ = 0.01 we 
find 



H 



~ 10""* (2.54) 

T=0.7 keV 



stems from the fact that their energy dependence has not been resolved; cf. [64]. 

Using instead of (2.49) the cross section {crex,4Ho'^) inferred from Table 2.4 yields the coeffi- 
cient 4 X 10"^ in Eq. (2.52); see [4]. 
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which tells us that the whole issue of {pX ) mediated destruction of ^Li (and accordingly 
of ^Be) is irrelevant so that the abundances of ^Li and ^Be fused at T ~ 8 keV remain 
unaffected at T ~ 1 keV.^^ 

We can further employ the results summarized in Table 2.4. A successive chain 
of charge exchange reactions with (pX^) can lead to molecular states that are finally 
destroyed in nuclear reactions with protons. In particular, (^LiX^) has a chance for a 
nuclear interaction with protons or helium unsuppressed by a residual Coulomb barrier 
since it is a very compact object. We can solve for those molecular bound states by 
extending our network of Boltzmann equations for T smaller than a few keV with 

-HT-^ ('LiX") = {a,,euv) (pX") ^Li - Tx (^LiX") 

- (f^cx,(6LiX-)^) (pX-) C^LiX-), (2.55a) 



,(6LiX-) 

-HT-^ (6LiX2 ) = (<Tex,(6LiX-)^') «b (p^^) (^LiX") - 2Tx eUX- 



- (^cx,{6LiX^)^) (P^-) ('LiX^ ), (2.55b) 
-HT-^ eUX,) = (^ex,(6LiX^)^) ipX-) (6LiX2 ) - 3Tx (^LiXg ) (2.55c) 
and supplement the right hand sides of Eq. (2.50) with 

^(px-) = -(^ex,6Li^) rih {pX') ^Li, -(o-ex,(6Lix-)^^) nh (pX") (^LiX") 

- (^cx,(6Lijq)^) (pX-) eUX^), (2.56a) 
= -(^Tex,6Li^) (pX-) ^Li + Tx [('LiX-) + 2(6LiX2 ) + 8(^^X3 )] (2.56b) 

A similar chain exists for ^Be where the sequence of the charge exchange reactions can 
proceed until (^BeX^). It is important to note that the efficiency of this chain reaction 
depends very sensitively on the concentration of the (pX~) bound states and on the 
mass of the X~ particle. The latter enters through the average relative velocity of two 
heavy objects, e.g., (pX~) and (^LiX"), which in turn scales as m^l/^. Therefore, in 
the limit of an infinitely heavy X~, the chain will be cut off right at the first step, 
terminating at (^LiX~). Also note that in (2.56b) we have made the assumption that 
^Li is not destroyed by the decay of X~ . Clearly some of the recoiling ^Li nuclei will be 
destroyed when released from the bound state. For the case of (^BeX~) this has been 
investigated in [75] where it was found that no significant depletion of ^Be takes place. 
On those grounds and noting that only a small fraction of ^Li is locked in bound states 
with X~ it is safe for us to disregard this effect. 



^^The same conclusion — prior to the publication of [64] — has already been reached in [4] by assigning 
the maximal possible rate for the destruction process (2.45) which is given by the unitarity bound. 
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X- ^ 10-2, mx- = 100 GeV, r^- = oo 




T [keV] 

Figure 2.7: Evolution of primordial abundances as a function of time (or temperature 
Tg) from the input X~ = 0.01, mx- = 100 GeV, and tx- ^ oo. The {pX~) abundance 
reaches its maximum of ^ 1.7 x 10^^ at T ~ 0.7 keV. Around the same temperatures, 
the abundance of unbounded CHAMPs, X^^^, starts to decline more rapidly since it is 
removed by the recombination with p followed by the charge exchange reaction on ^He. 

We run the full set of Boltzmann equations to determine the residual concentrations 
of (pX^) and of the molecular bound states of ^Li with X~ . The results are plotted in 
Fig. 2.6. As one can see, an initial concentration of per nucleon of X~ = results in 
a (pX^) abundance that never exceeds the maximum (2.52), leading to a progressively 
diminishing number of molecular states. 

We remark that for the computation of Fig. 2.6 we have not included {^LiX^) forma- 
tion via radiative recombination ^hi + X" (^LiX~) + 7. This process is not important 
for the present discussion because it does not affect (2.52); {pX^) is dominantly removed 
by charge exchange on ^He. We have checked this by including the recombination pro- 
cess into our reaction network; cf. Table 2.3. Whereas {^LiX~) would form around 
the same time as (^HeX~) it is important to note that {^LiX^) is also destroyed via 
(^LiX~) -|- p — > ^He -|- ^He -|- X~ hindering its formation; the corresponding cross section 
is given in [64] and in Table 2.5 with the name (o"cat,des2,6Li^^)- Even with -^g^. = 0.01 
we find that the final (^LiX"-) output is one order of magnitude below unbounded ^Li 
and (2.52) remains unchanged. By the same token, we find that the destruction process 
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only has a minor impact on the total ^Li abundance (below 10%). For smaller values of 
"^dec effect is accordingly weaker. 

To conclude this section, neither lithium nor beryllium synthesized in CBBN processes 
at 8 keV would be affected in any significant way by the subsequent generation of (pX^) 
bound states. Thus — as we will show later — the part of the parameter space with a 
typical freeze-out X~ abundance and a long X~ lifetime will be confidently ruled out. 



2.5.1 Relaxation after charge exchange 

Above we have stated that the charge exchange on "^He will mainly proceed into excited 
states of (^HeX^). Though we have already shown that (p^^) bound states do not reach 
an abundance level where late-time catalysis plays a role it is nevertheless amusing to 
see whether (pX^) receives further depletion by the energetic photons released from 
relaxation of ^He from its excited states n = 3 into the (^HeX~) ground state. Upon 
transition into the ground state ~ 300 keV will be injected into the plasma in form 
of photons. This happens predominantly at Tex — 0.7 keV and the photons loose their 
ability to break up (p-'^") once they are degraded below E"'^^ ~ 25 keV. 

The rate for photon-photon scattering with 7bg — scaling with E^ [91] — is rapidly 
becoming inefficient for E-y = O (100 keV) at the relevant temperature so that we can 
neglect this photon- multiplying process. Thus, the photons loose their energy mainly 
by Compton scattering on background electrons. In the low-energy limit of Thomson 
scattering the mean lifetime of a photon before scattering is 

1 Svra^ 
= ~ 10 s with (Tt = tt, (2-57) 

'^e-'^T Tex 2"^e 

where we have used rig- = 7/8r]]~,n^; see (1.45). 

Since the energy loss of 7 particles in Compton scattering is very small, we can write 
a differential equation for the systematic energy transfer to the electrons 

dE 

-^ = {aAE^)n,-. (2.58) 

Note that we have neglected a term HE^ which would account for the expansion of the 
Universe. We can do so since <C H{Tc^)~^ ~ 5 x 10^ s. The average energy loss per 
scattering is given by 



(aAE^) = I AE^^dn . (2.59) 
dil 
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At Tex the rest frame of the electrons essentially coincides with the frame of the thermal 
bath. Then da/dQ is given by the Klein-Nishina formula [92] 

2 



da 

dn 



a 



5 

2m2 \E-y 



with 



l + {E^/me){l-cose)' 
From (2.59) together with AE-y = E' — E^ and E^ <^ rric it then follows that 



(aAE^) ~ -cTT — . 



(2.60) 



(2.61) 



(2.62) 



Since the mean lifetime of the photon against Compton scattering is much shorter 
compared to the Hubble time we can integrate (2.59) from E"™ to E°^^ by neglecting 
any associated drop in temperature of the plasma. This defines a typical escape time 
Tesc, i-e. a thermalization time-scale upon which the photon looses its ability to ionize 
(pX-). We find 



Snip 



JaTVhn-yiTe, 



1 



Rout 

7 



1 



160 s. 



(2.63) 



This has to be compared with the mean lifetime of an energetic photon against 
ionization of (pX^), 



Tph 



1 



n 



(pX-)O-ph 



10^ s 



(2.64) 



and which is of the order of the Hubble time; In the last step we have used rif^px-) = 
'^(pX-) ^ given in (2.52) with X~ = 0.01. The photo-dissociation cross section has been 
obtained in Sec. 2.3.1; see Table 2.3. This tells us that the thermalization of the injected 
photon happens very rapidly so that those photons released in the relaxation process 
after charge exchange are not capable of depleting (pX~) any further. 



2.6 Constraints on the X lifetime and abundance 

In order to constrain the {tx,Yx-) parameter space from the catalytic path (2.40) 
to ^Be [63], we need to set an upper limit on its primordial abundance from existing 

^^The numerical values of the various times-scales sensitively depend on the exact value of fiducial 
temperature Tox chosen. However, the argument is not affected by that. 
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Figure 2.8: Observations of Be in Pop II halo stars. In the left panel (a), the data is 
taken from Fig. 3a of Ref. [97] and is plotted as a function of [Fe/H]. The right panel (b) 
shows the data from Fig. 6b of Ref. [99] where [O/H] provides the metallicity indicator. 
The filled dots depict the data points associated with the star G 64-12. The solid lines 
give the inferred nominal upper limits on ^Be from the weighted mean (dashed lines) 
of a sample of stars at lowest metallicity. Also shown in Fig. 2.7b is a fit of a primary 
scaling of Be; see main text. 



0.6 



observations. It is generally accepted that the galactic evolution of the abundances of Be, 
along with Li and B, are dominated by cosmic-ray nucleosynthesis. While Be is burned 
rapidly in stellar centers, it is produced in cosmic rays by the spallation reactions of 
fast protons and a particles hitting ambient CNO nuclei [93, 94]. As a consequence, the 
abundances of Be and O are linked, leading to a secondary scaling, Be oc [95]. On the 
other hand, inverse spallation reactions of CNO nuclei, both produced and accelerated 
in supernovae, will give a Be yield that is essentially independent of the metallicity of 
the interstellar medium. Such primary processes, leading to Be tx O, are expected to 
play a major role during the early galactic epochs [96] . 

The produced Be is subsequently supplemented in the outer layers of stars. Thus, 
old stars which are far from the galactic center (and thereby less affected by the galactic 
chemical evolution) bear the potential to encode any pre-galactic origin of Be. Indeed, 
Be has been observed in a number of Population II halo stars at very low metallici- 
ties [Fe/H] < —2.5. Particularly noteworthy is the detection in the star G 64-12 at 
[Fe/H] ~ -3.3 [97]. The star's high Be value of logio(Be/H) ~ -13.05 might sug- 
gest a possible flattening in the Be trend during the early evolutionary phases of our 
galaxy [97] . Whether this really points to a primordial plateau or whether this indicates 
a Be dispersion at lowest metallicities [98] is not clear at present. 
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Figure 2.9: Contour plot of CBBN abundance yields of ^Li and ^Be in the {tx, -'^doc) 
plane. The solid line shows the limit (2.65). The region above this line is excluded by 
^Be overproduction. The lower (upper) boundary of the band corresponds to ^Li/H = 
10~^^ (10^^*^). The y-axis on the right-hand side indicates the number density 
n^f normalized to the entropy density, Y^l'^. The cross shows the parameter point 
considered in Fig. 2.4. 

Figure 2.7a shows the original Be detection in the star G 64-12^^ (filled dot) along 
with a subset of data points taken from Fig. 3a of Ref. [97]. The data of Fig. 2.7b 
are taken from Fig. 6b of the recent work [99] which also uses [0/H] as a metallicity 
tracer. The latter paper discusses the implications of a new temperature scale on the 
abundances of Li, Be, and B. In principle, different assumed physical parameters which 
characterize the stellar atmosphere may result in large systematic shifts of the inferred 
abundances. In this regard, it is important to note that Be is not overly sensitive to 
the assumed surface temperature of the halo dwarfs [99]. In the following we thus shall 
take a pragmatic approach: In both Fig. 2.7a and Fig. 2.7b, we obtain the least squares 
weighted mean (dashed lines) for a representative sample of stars at lowest metallicities. 
From the variance of the fit, we can extract a nominal 3cr upper limit (solid lines) on 

^''For consistency with the rest of the data points, the ID LTE value has been plotted in Fig. 3a of 
the original reference [97]. 
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primordial ^Be. From Fig. 2.7b, we find 

logioBe/Hlhigh = -12.68 ^ ^Be/H < 2.1 x IQ-^^ . (2.65) 

Conversely, Fig. 2.7a yields ^Be/H < 10^^'^ while fitting only the last two data points 
with [0/H] < -1.3 in Fig. 2.7b would give ^Be/H < 1.3 x lO^^^. In our context, 
those values are less conservative so that we use (2.65) in the following. In Fig. 2.7b 
we have additionally fitted for a primordial component, ^Be/H|p, in combination with a 
primary scaling, ^Be/H = k (O/H)/(O/H)0. It seems, however, that a purely primary 
mechanism with k ~ 2.9 x 10~^^ fits the data best since ^Be/H|p comes out negligibly 
small. Finally, we are aware that neither of the fitted mean values in Fig. 2.7 is very 
good in terms of x^. However, a firm conjecture of a Be plateau is not the purpose of 
this work, and indeed (2.65) does provide a sufficiently conservative limit to work with. 

We can now confront the constraint (2.65) as well (1.6) with the CBBN yield of ^Be 
and ^Li obtained by solving the associated Boltzmann equations presented in Sec. 2.4 
for a wide variety of {'Tx^X'^^^ combinations. In Fig. 2.8 we obtain exclusion boundaries 
from catalyzed ^Be and ^Li production in the {tx-,X~^^^ parameter space. For conve- 
nience of the reader, the X~ number density normalized to the entropy density, 
y^ec ^ jg g^ygjj lY^Q y-axis on the right-hand side. Above the solid line, ^Be is in excess 
with respect to (2.65) and thus excluded. The shown band reflects the uncertainties in 
the observational determination of ^Li. On the lower border, ^Li/H = 10~^^ is fulfilled 
while ^Li/H = 10"^" holds on the upper border of the band. The cross indicates the 
exemplary parameter point considered in Fig. 2.4. At large lifetimes, the linear scal- 
ing of ^Li with X~ can easily be seen from the boundaries of the band. Note that we 
find ^Be/^Li in the interval between 10"'^ and 10~^, whenever CBBN is efficient, which 
confirms the observation already made in Ref. [63]. 



For a proper comparison between different assumed surface temperature scales and corresponding 
fits of primary versus secondary scaling, see Ref. [99]. 
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Table 2.6: Here we collect the key CBBN cross sections, i.e., the reaction rates per 
particle pair, iVA((Tw), used in the numerical solutions of the Boltzmann equations. 
They are given in units of cm^s~^niol~^ and Tg — T/IO^ K. The photo-dissociation 
cross sections ((Tphw) are related to the rates Fph via Fph = n-y(CTphw). 

process name rate [cm^s~^mol^"^] 

Recombination and photo-dissociation: 

p + X- ^ (pX-) + 7 {(T,ec,pv) 3980 Tg-^/^ 

{pX~) + 7bg ^ p + X- (fTph,p^^> 1.18 X 10^ T9-2 exp (-O.29/T9) 

^He + X- ^ (^HeX-) + 7 {(Trec,mev)^ 7900 Tg^^^ 

(^HeX-) + 7bg ^ ^He + X" {(Tph^Hev)^ 1-85 x lO^^ T^^ exp (-4.03/r9) 

^Li + X- ^ (^LiX-) + 7 {(Treceuv) 6640 Tg^^^ 

i^UX-) + 7bg ^ 6Li + X- (cTph,6Li^^) 2.87 x lO^o T'^ exp (-9.25/r9) 

Charge exchange: 
(pX-) + ^He ^ (^HeX-) + p (f^ex.^He^) 1-0 x 10^° [64] 

3.9 X lO^o Tg-^^ 

(pX-) + ^Li ^ (^LiX-) + p (f^ex,6Li^^) 6.45 X 10i° Tg^^ 

(pX-) + (6LiX-) ^ (^LiX^-) + p (^7e,,(6Lix-)^^) 3.37 X 109 Tg^/' (1 TeV/mx-)^/' 
(pX-) + (^LiX^-) ^ (^LiXg-) + p (a^,,(6Lix-)^) 5.25 x lO^ tJ/' (1 TeV/m^-)i/2 

^Xz destruction (from [64]): 
(pX-) + 6Li ^ + 3He + X" (acat.des.eLit-) 1-6 X 108 

(^LiX-) + p ^ 4jje + 3He + X" {a,,t,des2,'^uv) 2.6 x lO^^ Tg^'/' exp i-6.74T~'/') 

^Li and ^Be catalysis (from [84] and [63]): 

(^HeX-) + D ^ ^Li + X- ((Tcat,eLi^') 2.37 x 10^ (1 - 0.34 Tg) T~'^^^ 

X exp (-5.33 Tg"^/^) 

^He + (^HeX-) ^ (^BeX") + 7 {(Tc^t&Bev) 10^ Tg^^^ [ 0.95 exp (-1.02/r9) 

+0.66 exp (-1.32/r9) ] 

(SBeX-) + n ^ 9Be + X- {(Tc.t,^Bev) 2 x 10^ 

^ For consistency, the rates employed for (*HeX~) match the ones from [4]. From our numerical 
evaluation [Table 2.3] the respective coefficients for (urec ^hs'^) and {upi^ ijiev} read 7260 and 1.70 x 10^°. 
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Chapter 3 



Gravitinos as a probe for the 
earliest epochs 

3.1 The gravitino-stau scenario 

In the first part of tliis thesis we have discussed the implications of a generic electro- 
magnetically charged massive particle species if it is present in the early Universe 
during/after the era of BBN (t > Is). If X is a weak scale thermal relic, it freezes out 
from the primordial plasma at cosmic times t < 10^'' s so that the question of the origin 
of its longevity arises. 

Long-lived charged particles can naturally emerge in supersymmetric (SUSY) ex- 
tensions of the Standard Model. In scenarios in which the gravitino G is the lightest 
supersymmetric particle (LSP), a long-lived may be realized if the lighter stau ri is 
the next-to- lightest SUSY particle (NLSP). Assuming conserved i?-parity in this work 
the stau NLSP will be typically long-lived because it can only decay into the G LSP with 
Planck-scale suppressed couplings. Conserved i?-parity also implies that the gravitino is 
stable which makes it a promising dark matter candidate.^ 

In this part of the thesis we consider gravitino dark matter scenarios in which the 
Ti is the NLSP. In the present chapter we briefly introduce the gravitino and discuss 
immanent cosmological implications which are independent of the nature of the NLSP. 
In the next chapter we then work out the phenomenology of the gravitino-stau scenario. 

^The gravitino can also be dark matter if i?- parity is broken as long as it is ensured that the G?- lifetime 
is of the order of the age of the Universe; see, e.g., [100, 101, 102, 103] and references therein. 



62 



3.2. Supergravity and basic properties of the gravitino 63 



3.2 Supergravity and basic properties of the gravitino 

In ordinary gauge theories the generators of the of the Poincare algebra commute with the 
generators of the internal (local) symmetry such as, e.g., color SU(3)c in the Standard 
Model. Indeed, it was shown [104] that any such extension of the Poincare algebra 
in a four-dimensional quantum field theory (with non-zero scattering amplitudes) is 
necessarily trivial in the sense that both algebras decouple. 

Supersymmetry, however, is an extension of the space-time symmetry which relates 
fermionic and bosonic degrees of freedom. The associated particles form a supermultiplet. 
The spinorial supersymmetry generator Q obeys anii-commutation relations and the 
(simplest) supersymmetry algebra reads^ 

{ga,Q/3} = 27^^P^, (3.1a) 

{Qa,Q/3} = {Qa,Q/3} = 0, (3.1b) 

[^M, Qa\ = [P^.,Qp\ = 0, [P^, P.] = 0. (3.1c) 

Note that (3.1a) relates Q with the Poincare algebra; generates translations. There- 
fore, local supersymmetry implies local Poincare symmetry, i.e., invariance under general 
coordinate transformations. This is exactly what we expect from a theory of gravity so 
that local supersymmetry is also referred to as supergravity. The gauge field of super- 
gravity is the gravitino. It is a spin-3/2 Majorana particle and can be written as a 
vector-spinor if)^. 

Particles within the same supermultiplet are degenerate in mass because [P^, Qa\ = 0. 
Since we do not yet have experimental evidence for supersymmetry, we know that it has 
to be a broken symmetry if realized in nature. Local supersymmetry offers the appealing 
possibility to be broken spontaneously with a super-Higgs mechanism operating. The 
Goldstone fermion of supersymmetry breaking is absorbed by G which thereby acquires 
its longitudinal, helicity ±1/2 degrees of freedom. After supersymmetry breaking the 
gravitino has mass mg. Depending on the underlying breaking mechanism, can 
range from the eV scale up to scales beyond the TeV region [106]. 

The phenomenology of the massive gravitino is then governed by the following La- 
grangian 

C = -i^M-P-v^^^^^^SpVa - ^rngV^^[7^ ll^u + Ant. (3.2) 

The first two terms describe a free massive spin-3/2 field [107] from which it can be 
shown that the free gravitino field satisfies the Dirac equation {i$ — m^)'\pp^ = for each 

^We follow the conventions used in [105]; [ ■ , ■ ] and {■ , •} denote the commutator and anti-commutator, 
respectively. 
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component /i and is subject to the constraints "y^'ip^ = and d^ip^ = 0. The interaction 
Lagrangian reads^ 

Focusing on a minimal particle content in the observable sector, the fields (p, Xh, and A 
denote the gauge eigenstates of the scalars, chiral fermions, and gauginos of the MSSM; 
Fpa^^ is the field strength tensor of the gauge group (a) = (SU(3)c, SU(2)l, U(1)y) with 
index a of the associated adjoint representation and denotes the gauge-covariant 
derivative. All matter fields are written in terms of left-handed four-spinors Xh since 
they stem from left-chiral supermultiplets in the general supergravity Lagrangian [109]. 
For example, a right handed tau lepton is written in terms of its charge conjugate 
{t]^Y which is a left-handed spinor. Analogously the superpartner of is written in 
(3.3) as (tji)*. Details aside, most important is the fact that the interactions of the 
gravitino to the MSSM fields are fixed by the (super)symmetry and are suppressed by 
inverse powers of Mp which makes G an extremely weekly interacting particle.^ 

Being the gauge field of supergravity, the gravitino sits at the heart of any locally 
supersymmetric theory. Despite its extremely weak interaction gravitinos can be effi- 
ciently produced in the early Universe. In the next section we discuss the case of thermal 
gravitino production — a guaranteed source of potential relic gravitinos. 




3.3 Thermal gravitino production and reheating 

The observed flatness, isotropy, and homogeneity of the Universe suggest that its earliest 
moments were governed by inflation [111, 112]. The inflationary expansion is followed by 
a phase in which the Universe is reheated. The reheating process repopulates the Uni- 
verse and provides the initial conditions for the subsequent radiation-dominated epoch. 
The reheating temperature Tr can be viewed as the initial temperature of this early 
radiation-dominated epoch of our Universe. 

The value of Tr is an important prediction of inflation models. While we do not have 
evidence for temperatures of the Universe higher than 0{1 MeV) (i.e., the temperature 
required by primordial nucleosynthesis), inflation models can point to Tr well above 

^For an explicit "derivation" of dnt from the general supergravity Lagrangian see [108]. 
^Couplings of a very light gravitino can be enhanced due to its longitudinal (goldstino) modes [110] — a 
remnant of the super-Higgs mechanism. 
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Table 3.1: The gauge couplings gi and the constants Ci, ki, yi, and (3^^^ associated 



je groups U(1)y, 


SU(2)l 


, and SU(3)c. 






gauge group 


i 9i 


Ci 


ki 




A'" 


U(1)y 


1 9' 


11 


1.266 


0.653 


11 


SU(2)l 


2 9 


27 


1.312 


1.604 


1 


SU(3)c 


3 9s 


72 


1.271 


4.276 


-3 



lO^'^ GeV [112, 113]. While any initial population of gravitinos must be diluted away by 
the exponential expansion during inflation [114], gravitinos are regenerated in scattering 
processes of particles that are in thermal equilibrium with the hot primordial plasma. 
The efficiency of this thermal production of gravitinos during the radiation-dominated 
epoch is sensitive to Tr [115, 116, 117, 118, 108, 119, 120]. 

Gravitinos with > 1 GeV have decoupling temperatures of Tf' > 10^^ GeV, as will 
be shown below. We consider thermal gravitino production in the radiation-dominated 
epoch starting at Tr < Tf' assuming that inflation has diluted away any initial gravitino 
population.^ For Tr < Tp, gravitinos are not in thermal equilibrium with the post- 
inflationary plasma. Accordingly, the evolution of the gravitino number density with 
cosmic time t is described by the following Boltzmann equation [108, 119] 

dnp^ 

+ Wng = Cg (3.4) 

c.^y:SSfi+^l..?i..m (3.5) 



G 167r3M2 I 3ml i \ g, 

i — 1 \ G / 

The collision term involves the gaugino mass parameters Mj, the gauge couplings gi, 
and the constants q and ki associated with the gauge groups U(1)y, SU(2)l, and SU(3)c 
as given in Table 3.0. In expression (3.5) the temperature T provides the scale for the 
evaluation of Mj and gi. The given collision term is valid for temperatures sufficiently 
below the gravitino decoupling temperature, where gravitino disappearance processes 
can be neglected. A primordial plasma with the particle content of the MSSM in the 
high-temperature limit is used in the derivation of (3.5). 

The collision term (3.5) results from a consistent gauge-invariant finite-temperature 



^In this work, we neglect gravitino production in inflaton-decays, cf., e.g., [121, 122, 123] and references 
therein. Though this non-thermal source can give a sizable contribution it is model-dependent and 
typically small when considering high reheating temperatures with the associated gravitino yield scaling 
as Tj^^ [123]. 
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Tr [GeV] 

Figure 3.1: The thermally produced gravitino yield (3.6) as a function of Tji for 
mg = 10 MeV, 100 MeV, 1 GeV, 10 GeV, 100 GeV, and 1 TeV (from left to right) and 
A^i.2,3(Mgut) = 'Tii/2 = 500 GeV. The dashed horizontal line indicates the equilibrium 
yield of a relativistic spin 1/2 Majorana fermion. 



calculation [108, 119] following the approach used in Ref. [118]. Thus, in contrast to 
the previous estimates in [115, 116], the expression for is independent of arbitrary 
cutoffs. Note that the field-theoretical methods of [124, 125] applied in its derivation 
require weak couplings, ^ 1, and thus high temperatures T ^ 10^ GeV.^ Thus, in the 
following we focus on cosmological scenarios with Tr > 10^ GeV which is also the most 
attractive temperature range, e.g., for baryogenesis scenarios based on leptogenesis. 

Assuming conservation of entropy per comoving volume, the Boltzmann equation (3.4) 
can be solved to good approximation analytically [118, 126]. At a temperature Tiow 
Tr, the resulting gravitino yield from thermal production reads 



where the constants are given in Table 3.0. These constants are obtained with an 
effective number of relativistic degrees of freedom of ^ofrC^R) = ^efrl^R,) = 228.75 which 



For an alternative approach see [120]. 
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fohows from the fact that the entire MSSM particle content is in thermal equilibrium 
and relativistic. We evaluate gi(TR) and Mi{T^) using the one-loop evolution described 
by the renormalization group equation in the MSSM [106]: 

-1/2 



9^m 

Mi{T) 



9^{T) 



(1) 



In 



T 



GUT J 



87r2 

2 



(3.7) 
(3.8) 



with the respective gauge coupling at the Z-boson mass, gi{Mz), and the coefficients 
listed in Table 3.0. For convenience we choose to parameterize the gaugino masses in 
terms of their values at the scale of gauge coupling unification Mqut — 2 x 10^^ GeV. 
We remark in passing that when considering only the SUSY-QCD contribution in (3.6) it 
is sometimes also convenient to express the gravitino abundance in terms of the physical 
gluino mass [118]. As has been shown in [127] it is then important to employ a two- loop 
running of the gluino mass since using (3.8) for the renormalization group evolution from 
the electroweak scale to Tr, would underestimate the gravitino abundance by a factor of 
two. In this work, however, we use the running gluino mass M3 at Mqut as input where 
the effect is smaller (and working in the other direction.) 

For a standard cosmological history without release of entropy, the gravitino yield 
from thermal production at the present temperature Tq is given by yJ^(To) = Y'^^{T\o^). 
The resulting density parameter of thermally produced gravitinos reads 



nfh^ = y|P(ro) sin) h^/pc 

with /9c/[s(ro)/i2] = 3.6 X 10-9 GeV [32]. 



(3.9) 



In Fig. 3.0 the result (3.6) for the thermally produced gravitino yield Y~ {T\ 



low J 



is shown as a function of Tr for various values of (solid lines). The curves are 



obtained with m 

Mi,2,3(Mgut) 



1/2 = 500 GeV for the case of universal gaugino masses at Mqut: 
mi/2- The dashed (blue) horizontal line indicates the equilibrium yield 
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(3.10) 



which is given by the equilibrium number density of a relativistic spin- 1/2 Majorana 
fermion, = 3C(3)rV(2^^). For T > , ges{T) = Kg{T) = 230.75 since the spin- 
1/2 components of the gravitino are in thermal equilibrium. In the region where the 
yield (3.6) approaches the equilibrium value (3.10), gravitino disappearance processes 
should be taken into account. This would then lead to a smooth approach of the non- 
equilibrium yield to the equilibrium abundance. Without the back-reactions taken into 
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account, the kink position indicates a lower bound for Tf. Towards smaller mg, 
decreases due to the increasing strength of the gravitino couplings. For example, for 
= 1 GeV (10 MeV), we find Tp > 10^^ GeV (10^° GeV). 



3.3.1 Reheating phase 

In the analytical expression (3.6) we refer to as the initial temperature of the 
radiation-dominated epoch. So far we have not considered the phase in which the co- 
herent oscillations of a field (p dominates the energy budget of the Universe and where 
one usually relates Tr to the decay width of (p. In the simplest models of infiation 
the decaying field (p which reheats the Universe also drives the exponential expansion of 
the Universe. In the following we shall simply refer to (f) as the infiaton field. 

We can account for the (perturbative) reheating phase, by considering (3.4) together 
with the Boltzmann equations for the energy densities of radiation and the infiaton field, 



^^^^ + 4F/9rad = r<^p0 , (3.11a) 
^ + 3Hp^ = -r^p^ , (3.11b) 

respectively. To relate with Tr we first note that the second term on the left hand side 
of (3.11b) indicates that (p (when averaged over several oscillations) scales like matter, 
P(f) tx a~^; a denotes the scale factor. Thus, in terms of the initial infiaton energy density 
P(f)^i, the Hubble rate is given by 



as long as dominates. Assuming an instantaneous conversion of p^f^ into radiation 
when = ^H(Tji) with usually chosen to be a number between 1 and 3 then allows 
one to define a reheating temperature in terms of the decay width of the infiaton field, 

In Fig. 3.1 we show the results of a numerical integration of (3.4) and (3.11) plotted 
against the scale factor with F^ chosen such as to yield Tr = 10^ GeV for ^ = 1.8."^ 
In the lower figure we plot comoving quantities normalized to afp^pj. At a/aj ~ 10^ 



^For the actual integration we rewrite the Boltzmann equations (3.4) and (3.11) in terms of dimen- 
sionless quantities following [128]. 
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Figure 3.2: The results of a numerical integration of (3.4) and (3.11) is shown; 
is chosen such as to yield Tr = 10^ GcV for ^ ~ 1.8. In the lower figure comoving 
quantities as labeled and normalized to a^Ptj,,i are plotted. The middle figure shows 
the evolution of the temperature T of the thermal bath and in the top part of the 
figure the resulting gravitino abundance is obtained for = 100 GeV and for 
-^'^1.2,3 = ^1/2 at Mqut with TOi/2 = 500 GeV. This gives a final gravitino abundance 
of Y~^{Tiovj) = 2.9 X 10^^^; see main text for a discussion. 
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the inflaton energy density decays exponentiahy and the Universe enters the radiation 
dominated epoch which can be seen by the turn-over of the dashed curve depicting the 
radiation density: there, p^d starts to scale as a~^. This is also indicated by the evolution 
of the temperature plotted in the middle part of the figure. In the radiation dominated 
Universe T scales with a~^. Also note that Tr is not the maximum temperature of the 
Universe. When the inflaton field decays entropy is produced (dotted line in lower part) 
and the temperature scales as T oc a^^^^. The production of entropy is also the reason 
why in the top part of the figure — despite the large initial temperature — reaches 
its maximum value only in the radiation dominated regime when a^s is finally constant. 
Here, Mi^2,3 = 'm'i/2 af ^gut with mi/2 = 500 GeV and = 100 GeV has been 
chosen which gives a final gravitino abundance of ^-"^(Jiow) = 2.9 x 10^^^. 

With our result for the collision term (3.5), we find that the gravitino yield obtained 
numerically is in good agreement with the analytical expression (3.6) for ^ = 1.8. For 
an alternative Tr definition with different ^ the associated numerically obtained gravi- 
tino yield is described by the analytical expression obtained after substituting Tr with 
VU^T^ in (3.6). 

A fitting formula on the gravitino yield which also includes the effect of reheating 
and which was based on [118] has been derived earlier in [47]. However, the production 
of the helicity-1/2 component of the gravitino was neglected so that the actual yield for 
mg = 100 GeV was underestimated by about an order of magnitude. Accordingly, the Tr 
bounds given in [47, 41, 129] are underestimated in the region mg < 1 TeV. Meanwhile — 
after publication of [1] on which this section is based on — an updated treatment [130] 
has become available in which the authors now include the helicity-1/2 components as 
well as the electroweak contributions [108, 119] of the thermal gravitino production. 
For a most recent discussion on gravitino production during perturbative reheating see 
also [131]. 



3.4 Constraints on Tr 

Since for the gravitino LSP the resulting density should not exceed the dark matter 
density f^dmi 7r is bounded from above [116]. Such a bound has to be compared with 
predictions of the reheating temperature Tr, from inflation models. Moreover, Tr is im- 
portant for our understanding of the cosmic baryon asymmetry. For example, successful 
standard thermal leptogenesis [133] can typically require Tr > 10^ GeV [134, 135]. 

We update the Tr limits using the full gauge-invariant result for the relic density 
of thermally produced gravitinos, , to leading order in the Standard Model gauge 
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Figure 3.3: Upper limits on the reheating temperature Tr. On the upper (lower) gray 
band, fl^^ for Afi,2,3 = = 500 GeV (2 TeV) at Mqvt agrees with n^'^. The 

corresponding Tr limits from the requirement fi-^/i^ < 0.126 shown by the dashed 
lines for A/i/10 = M2/2 = M3 = mi/2 at A/qut and by the dotted lines for the SU(3)c 
contribution [118, 132] with M3 ~ rn-iji at Mqut- All lines are obtained with (3.6). 



couplings [108, llO].*^ In particular, this allows us to illustrate the dependence of the 
bounds on the gaugino-mass relation at the scale of grand unification Mqut- 

The reheating temperature Tr is limited from above in the case of a stable gravitino 
LSP since cannot exceed the dark matter density ridm- In the following we use [7, 32] 
the WMAP 3-year result 

^^^J? = 0.105^^21 (3.14) 

In Fig. 3.2 we show the resulting upper limits on Tr using ^ = 1.8 as a function of mg. 
On the gray band, the thermally produced gravitino density (3.9) is within the nominal 
3(7 range (3.14). The upper (lower) gray band is obtained for Mi^2,3 = ''7^1/2 at Mqut 
with mi/2 = 500 GeV (2 TeV). Prom the requirement ^i^h^ < 0.126 the dashed lines 
show the constraints for the exemplary non-universal scenario [136] Mi/10 = M2/2 = 
= mii2 at Mgut- Using the same requirement the dotted lines show the SU(3)c 

*The computation of the electroweak contributions to thermal gravitino production was subject of 
the diploma thesis [108] of the author in which similar limits on Tr were already presented — however, 
without insight on the exact sensitivity of fi^^ on the reheating process; in this section ^ — 1.8. 
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contribution [118, 132] for M3 = mx/2 -^GUT- As can be seen by comparing the 
dashed and dotted hnes, the electroweak contributions can be particularly important 
for the case of non- universal gaugino masses at Mgut- On the other hand, for gaugino 
masses which unify at Mqut, the dotted lines provide already a very good estimate. In 
this regard, note that Fig. 3.2 is updated from the one presented in [1] since the SU(3)c 
contribution to fi-^, originally obtained in [118], has meanwhile been corrected by the 
authors [132]. For completeness, we also remark that, after publication of [1], the upper 
limits on Tr shown in Figs. 5 and 6 of Ref. [137] have also been corrected. Previously, 
these figures underestimated the maximal value of by a factor of four in the region 
in which fi-^ governs the limits. 

The Tr limits shown in Fig. 3.2 are conservative bounds that do only depend on 
mg and the Mj values at Mgut- Taking into account contributions to from NLSP 
decays will make those limits stronger. In the next chapter, we will account for this 
non-thermal gravitino source in a systematic way. 



Chapter 4 

The stau as the NLSP 



In this chapter we now speciahze on gravitino dark matter scenarios in which the hghter 
stau Ti is the NLSP. Indeed, the appearance of ri as the Hghtest Standard Model su- 
perpartner is a commonplace occurrence even in models with restrictive assumptions 
on the SUSY breaking sector such as the CMSSM. The associated parameter region is 
usually not considered because of severe upper limits on the abundance of massive stable 
charged particles (see Sec. 1.3). However, in gravitino LSP scenarios ti is unstable and 
thereby a viable option. 

In this chapter, we first review the result on a frequently used range of thermal freeze- 
out abundances of ri . Employing such an estimate allows us to constrain the gravitino- 
stau scenario from BBN limits on the electromagnetic and hadronic energy release in the 
decay of ri in a rather model-independent fashion. Moreover, we employ the results of 
Chapter 2 on the catalyzed light element production of ®Li and ^Be. We shall see that the 
associated constraints pose the henceforth strongest limits on this scenario. Specializing 
to the case of the CMSSM allows us to explore concrete realizations of the gravitino- 
LSP stau-NLSP setting. In particular, we shall find that the reheating temperature Tr 
is heavily constrained by the novel CBBN bounds. We also explore the possibility of a 
non-standard cosmological history to see whether one can alleviate or even circumvent 
the strong restrictions on the parameter space. 



4.1 Generic constraints on the gravitino- stau scenario 

In order to set constraints on the outlined scenario we require knowledge on the stau 
abundance in the early Universe. In Sec. 1.3 we have found that for a standard cos- 
mological history the decoupling temperature of a weak scale charged particle satisfies 
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Tf < 7X1^^12^. Thus, with a post-inflationary reheating temperature Tr above the decou- 
phng temperature, the t\ NLSP freezes out of the primordial plasma as a cold thermal 
relic so that its yield after decoupling Y-^'^ is governed by its mass and its annihilation 
rate. Thereby, Y-^'^ becomes sensitive to the mass spectrum and the couplings of the 
SUSY model and representative values 

yi" .(0.4^ 1.5) X 10- (^5^) (4.1) 

have been used to confront the gravitino-stau scenario with cosmological constraints [31, 
138, 139, 2, 140]; m^^ denotes the mass of the lighter stau. Equation (4.1) compares well 
with our upper bound (1.15) derived from a purely dimensional analysis. The yield (4.1) 
with a coefficient 0.7 x 10^^^ is in good agreement with the curve in Fig. 1 of Ref. [31] 
that has been derived for the case of a purely 'right-handed' r ~ tr NLSP with a mass 
that is significantly below the masses of the lighter selectron and the lighter smuon, 
<C ?7iei,j:i, and with a bino-like lightest neutralino, x5 — that has a mass of = 
\.\m^. In the case of an approximate slepton mass degeneracy, ra^ < ra-^^ -p^^ < 1.1 m:j:, 
the upper value in (4.1) becomes saturated due to slepton coannihilation processes [31]. 
We shall see in Sec. 4.2.3 that approaching the Xi""^ coannihilation region, m^o ~ vtl^^^ 
even larger enhancement factors occur. On the other hand, a sizable left-right mixing 
of the stau NLSP is associated with an increase of its MSSM couplings and thus with 
a reduction of Y-^'^ . This will be discussed in Part III where a systematic investigation 
of the stau abundance and its sensitivity on the SUSY parameters will be given. In this 
section, we shall focus on the more generic Y-^'^ values described by (4.1). 

In Fig. 4.0 we collect the cosmological constraints on the gravitino-stau scenario by 
plotting versus m-^f^. Let us go through the respective limits one by one: 

Non-thermal gravitino production Each t\ NLSP eventually decays into one G LSP 
leading to a non-thermally produced (NTP) gravitino density [31, 141]: 

r!NTP;^2 ^ ydec 5(ro)/i2/pc. (4.2) 

This contributes to the relic gravitino density J7g which should not exceed the 
observationally inferred dark matter density 0dm- In Fig. 4.0 we choose as repre- 
sentative value 

yii<=c = 0.7 X 10-13 ( ""^1 \ . (4.3) 
VlOOGeV/ ^ ' 

In the light shaded region in the upper right corner agrees with J^^m^^ 

of Eq. (3.14). The shading is limited from above by a solid line which borders 
the above disfavored region in which 17 -"'"^ > Odm- Any additional contribution 
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Figure 4.1: Constraints of tlie gravitino LSP stau NLSP scenario for Y:^°'^ given 
by (4.3). In tlie liglit (gray) shaded region fi-"""^ G fi^^ holds with the region above 
being disfavored by fig > ^2^^. The thin gray straight hnes show the contours on 
which /rJdm is provided by r?-"""^. The dotted gray lines show contours of as 
labeled. Catalyzed BBN production of ^Be disfavors the the region below the thick 
solid line. In the dark shaded (blue) region ^Li/H £ 10^^^ 10^^° holds. Below, ^Li 
is overproduced due to the bound state effects. Hadronic energy release in stau decays 
disfavors the regions inside the dashed lines for different values of the adopted primordial 
D abundance (4.5). By the same token, the regions inside dash-dotted (green) curves 
are disfavored from the effect of electromagnetic energy release on D. Moreover, the 
region below the thin solid (pink) line is disfavored from overproduction of "^He/D. 



to fidnn such as a thermally produced gravitino density Jl- (Sec. 3.3), makes 
this constraint more restrictive. This is indicated by the thin gray lines labeled 
with / = 0.1, 0.01, 10"^, and 10""^, on which (4.2) obtained with (4.3) satisfies 
f O^'^P = 0.126, respectively. The timing of the ri NLSP decay into the gravitino 
LSP is governed by the two-body decay mode t\ Gt and reads 

T,,c,r-Hn^Gl) = f^il-^) . (4.4) 

n \ n / 

in the limit nrir — > 0. Contour-lines thereof are shown in Fig. 4.0 by the dotted 
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gray hnes. Starting from the upper left edge, they correspond to hfetimes of 10 ^ s 
and then — as labeled — to 0.1 s, 10 s, 10^ s, lO^s, and 10'' s. 

CBBN constraints We can now incorporate the results of Chapter 2 on the catalyzed 
fusion of ^Li and ^Be triggered by the bound state formation of (^Herj"). We use 
the CBBN constraints obtained in Fig. 2.8 with y|!.= = F-f 72, I.e., we assume 
that there exists no asymmetry between positively and negatively charged staus 
[cf. Sec. 5.6]. For Y-^'^ we use the estimate (4.3). The shaded region in Fig. 2.8 has 
corresponded to a ^Li output of ^Li/H = 10~^^ -i- 10"^*^ and which is now likewise 
associated with the dark (blue) shaded band in Fig. 4.0. On the upper border of 
the band a lithium abundance of ^Li/H = 10~^^ is attained whereas on the lower 
border ^Li/H = 10"^° holds. The region below the band is confidently ruled out 
by overproduction of ^Li.^ Moreover, we also show the CBBN constraint from 
primordial ^Be production which excludes a very similar region. 

Hadronic energy release In Sec. 1.4 we have provided an overview over the physics 
of late decaying particles during/after BBN. In our concrete scenario we can now 
implement the stringent constraint on hadronic energy release from the observa- 
tionally inferred primordial deuterium abundance.^ The limits are based on the 
severe and conservative upper bounds on the product -Evis^NLSP [here, Inlsp = ^r] 
obtained in Fig. 39 of [41] for (see references cited in [41]): 

(D/H)^^^„ = {2.78toit) X 10"^ ^ severe constraint, (4.5a) 
(D/H)high = (3.98l[];g7) x 10"^ =^ conservative constraint. (4.5b) 

Recall from our discussion in Sec. 1.2 that (4.5b) is a rather high value on D/H|p. 
Without trying to give extra credence to (4.5b), following [41], we simply take it 
as a limiting value for D/H. The average injected hadronic energy E'vis has been 
obtained in [139] from computation of the 4-body decay of the stau NLSP into the 
gravitino, the tau, and a quark-antiquark pair for a purely right-handed ti ~ tr 
NLSP. The effect of hadronic energy injection on primordial D disfavors the regions 
inside the dashed lines shown in Fig. 4.0; see also Fig 16 in [139]. 

^ Using the initial estimate [79] on the CBBN ®Li output, the associated constraint has first been 
shown in the (m^, ^ )-plane in [142]. Likewise, in this representation, the constraints from hadronic 
Ti-decays have first been obtained in [141]. However, both works [141, 142] are based on outdated 
light element yields so that we use our treatment of Chapter 2 for "Li and the update in [139] for the 
hadronic ri-decays. 

■^Additional constraints on hadronic energy release are imposed by the primordial abundances of ''He, 
^He/D, ■^Li, and '^Li/'^Li [143, 45, 40, 41, 144, 145]. However, in the region allowed by the ^Be and ''Li 
constraints from bound-state effects, i.e., Tf^ < few x 10^ s, the considered D constraint on hadronic 
energy release is the dominant one as can be seen, e.g., in Figs. 38-41 of [41] and in Figs. 6-8 of [144]. 
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Electromagnetic energy release The constraints resulting from the dissociation of 
hght elements due to interaction with the electromagnetic cascades formed in stau 
decays are obtained for a "visible" electromagnetic energy of E'vis = 0.3 E-r of the 
tau energy 

o 9 9 

m~ — m~ + mt 

released in r — > Gt. The and ^He/D constraints result from the -Evis^NLSP 
hmits given in Fig. 42 of Ref. [41] and the Dg^*^ constraint from the -Evis^LSP hmit 
given in Fig. 6 of Ref. [39] . It is the region to the right or inside of the dot-dashed 
(green) curves and the region to the right of the and thin sohd (pink) line that are 
disfavored by the primordial abundances of D and ^He/D, respectively. 

As can be seen in Fig. 4.0, the constraints from the catalytic production of ^Li and 
^Be are (essentially) the most restrictive ones. Their coinciding position is due the fact 
that the catalytic production of ^Li and ^Be both depend on (^Her^") (same timing) and 
that their output scales linearly in 1~. Indeed, the associated constraints run parallel and 
in vicinity of the T^-^ = 10'^ s contour which corresponds to the time at which (^Herj") 
formation starts to become efficient. For the adopted stau abundance (4.3) this implies 
that lifetimes in the vicinity of Tf^ ^ 6 x 10'^ s are disfavored. 

The electromagnetic Dgm and ^He/D constraints are always less restrictive than the 
CBBN constraints from ^Be and ®Li. Only the hadronic constraint D^ad competes with 
the CBBN constraints for m:^^ > 1 TeV, i.e., in the lifetime region r:^^ ~ 10'^ s in 
which is largest. Though for lifetimes shorter than about 100 s neutron-to-proton 
interconversion processes affect ^He, the associated constraint is about two orders of 
magnitude weaker (as can be seen in left panel of Fig. 1.0). We also remark that the 
elevated content of D due to hadronic and electromagnetic energy injection leads to an 
enhancement of CBBN-produced ^Li and ^Be. For example, if non-thermal processes 
boost the deuterium abundance to the level of (4.5b), it would lead to an enhancement 
of the ^Li output by a factor of ~ 2, while the corresponding enhancement factor in the 
case of ^Be is about 4. We have not included this effect since it can make our obtained 
limits only stronger. 

Here we would like to emphasize that the ^Be and ^Li constraints are the ones that 
are the least sensitive to the precise value of Y-^'^ in the region Y^^'^ > lO^^'*. This results 
from the fact that the limits are very steep in that region, as can be seen in Fig. 2.8. 
Indeed, a yield that is twice as large as (4.1) will affect the position of the ^Be and ^Li 
constraints only very mildly. In contrast, such an enhanced yield — as encountered, e.g., 
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in the case of slepton coannihilations — leads to significant changes of the dark matter 
constraint and the BBN constraints associated with hadronic/electromagnetic energy 
injection, as can be seen exphcitly in Fig. 16 of Ref. [139]. 

It should also be noted that an elevated slepton yield can lead to an additional non- 
thermal output of ®Li for r^^ > few x 10^ s. As discussed before, this is because energetic 
spallation debris of destroyed ^He nuclei from slepton decays can hit ambient '^He and 
thereby fuse ^Li [40, 41, 144]. This mechanism depends sensitively on the hadronic 
branching ratio i?h of the 4-body slepton decay into the gravitino, the associated lepton, 
and a quark-antiquark pair for which typically < 3 x 10^^ for m:^^ < 2 TeV (see 
Fig. 5 of Ref. [139]). Indeed, as discussed in Ref. [146], for those branching ratios, the 
effect of CBBN on ^Li is the dominant one in the region which is not already excluded 
by the D constraint. Thus, for m^r^ < 1.5 TeV, our obtained limits on ^Li overproduction 
are only marginally affected by the hadronic energy release of ri-decays. However, for 
larger slepton masses, i.e., for scenarios of large YJ^^^ in conjunction with B\i > 10"^, 
the hadronic production of ^Li becomes efficient so that only a simultaneous treatment 
of both effects can decide on the accurate ^Li BBN output.^ Note that this can make 
our presented limits on ^Li only stronger. Thus, we are on the conservative side when 
neglecting such additional contributions. We note in passing that with a highly fine- 
tuned rrir^-mQ degeneracy leading to E'vis — > 0, any bound on energy release can be 
evaded. However, the CBBN bounds remain. 

Let us also comment on the reliability of the novel CBBN constraints from ^Li and 
^Be overproduction and address the implications of the associated restrictions on the 
{m^,mr^) parameter space: 

Reliability of CBBN constraints As already emphasized in Sec. 1.1, observations of 
^Li are extremely difficult. Whereas in the cold interstellar medium the lines of ^Li 
and ^Li are well resolved, measurements of the isotopic ratio ^Li/^Li in the outer 
layer of stars is complicated because the absorption lines of ^Li are not resolved 
spectroscopically with respect to the lines of '^Li due to thermal blending [18]. 
Indeed, the claim of a "^Li-plateau" is being challenged in the recent papers [147, 
148]. The presence of ^Li is inferred from an enhancement of the "red wing" of 
the ^Li absorption line. It is argued that such a line asymmetry could also be 
mimicked by Doppler-shifts due to atmospheric convective motions. Accordingly, 
some of the observations may eventually turn out to provide only upper limits. 

^Using the catalysis of BBN reactions in order to seek for a simultaneous solution of both, the ^Li 
and '^Li problem [cf. Sec. 1.1] has, e.g., been made in [145, 90, 75]; for a most recent discussion which 
also includes the lithium output due to X decays see [42]. 
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Moreover, ^Li is more fragile than ^Li and would burn more efficiently at lower 
temperatures. Therefore, if there is a (yet unconfirmed) stellar mechanism (see, 
e.g., [149]) that resolves the lithium problem, i.e., that depletes ''Li by a factor of 
two or three, ^Li would have been depleted by a larger factor. 

Given those issues, we have adopted here a generous range on the observationally 
inferred upper limit on primordial ^Li with ^Li/H = lO"^'^ being a very high value. 
In this regard it is also important to note that the catalytic effect on ®Li is very 
strong so that — for the purpose of setting constraints — we are not overly sensitive 
to the precise value of the upper bound. This can be seen by the fact that the dark 
shaded (blue) band in Fig. 4.0 is rather thin while spanning one order of magnitude 
in fused ^Li. In this respect. Sec. 2.5 becomes important in which we show that 
the destruction of large fractions of the previously synthesized ^Li by (p Tj" ) is not 
feasible. This rules out the possibility that allowed islands in the parameter region 
with large /large r:j;^ — which was advocated to remain viable in Ref. [90] — exist. 

Unlike ®Li, ^Be is firmly detected in a significant number of stars at low metallicity, 
and its observational status is not in doubt. (For the latest data on the ^Be 
abundance in metal-poor stars, see, e.g., [97, 150, 151, 98].) Also note that stellar 
depletion would affect ^Be less than either ^Li or ^Li since both ^Li and ^Li are more 
fragile than ^Be. Moreover, the nuclear physics rates that enter in the calculation 
of ^Be catalysis are dominated by resonances. Given the wealth of experimental 
information on the ^Be resonances [152, 153], this may eventually allow for very 
reliable calculations of the catalytic rates. Though it has recently been argued 
that the resonance in the final step (2.41b) in the fusion of ^Be is shifted below 
threshold [cf. Sec. 2.4.2] it is a neutron induced reaction so that such a shift will 
affect the efficiency of the reaction but may not be fatal. Moreover, also note that 
we have adopted a very conservative upper limit on primordial ^Be in (2.65) which 
can already be seen by mere optical inspection of Fig. 2.7b. Taking as a grain of 
salt that the final efficiency of the ^Be reaction is not fully established but noting 
the powerful physics potential a primordial origin of ^Be offers we have chosen to 
incorporate this constraint. Moreover, given the fact that both, ^Li and ^Be, show 
the same sensitivity on the gravitino-stau parameter space any conclusions drawn 
from ®Li are only corroborated and not altered by ^Be. We eagerly await further 
investigation of the critical catalyzed nuclear rates which shall give a final answer 
on the CBBN output of '^Be. 

Implications of CBBN constraints From Fig. 4.0 it is immediately clear that the 
new CBBN constraints imply a lower limit on given m^. Clearly, a minimum 
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value of the hghtest Standard Model superpartner directly affects the testability 
of such SUSY scenarios at future colliders. 

For the gravitino dark matter scenario in which i^dm is exclusively provided by 
j^NTP^ i.e., which are situated in the light gray shaded band, formerly (marginal) 
allowed islands with 200 GeV ^ rn^ < 400 GeV and m^-^ < 1.5 TeV are now 
confidently ruled out. Moreover, considering other values than / = 1 the gravitino 
mass mg is constrained to values well below 10% of the slepton NLSP mass m^f^j for 
m:^j < 0(1 TeV). Thereby the kinematical determination of proposed in [154] 
remains cosmologically disfavored at the next generation of particle accelerators. 

Of course, the CBBN constraints only emerge if r:^^ is large enough to allow for 
(^He-rf) bound state formation. Thereby, scenarios with a gravitino mass of 
< 200 MeV and m^-^ > 80 GeV — the latter of which is supported by the 
non-observation of long-lived charged sleptons at the Large Electron Positron Col- 
lider (LEP) [32] — are unconstrained from CBBN. Accordingly, for gauge-mediated 
SUSY breaking which typically predicts small values of mg, the CBBN constraint 
can be irrelevant. However, in gravity-mediated SUSY breaking the gravitino 
mass sets the scale for the soft breaking parameters so that > 10 GeV are 
the most natural values. Then, the CBBN constraints impose a lower limit of 
m^^ > 400 GeV. 

To extract further implications from the gravitino-stau scenario we resort in the 
next section to concrete supersymmetric realizations by full specification of the SUSY 
parameters. Thereby we obtain further insight on the superparticle mass spectrum. 
Moreover, this allows us obtain a stringent upper bound on the reheating temperature 
of the Universe. 

4.2 The gravitino-stau scenario in the CMSSM 

We now consider gravitino dark matter scenarios in the framework of the CMSSM where 
one assumes universal soft SUSY breaking parameters at Mgut- The CMSSM yields 
phenomenologically acceptable spectra with only four parameters and a sign: the gaugino 
mass parameter mi/2j the scalar mass parameter mo, the trilinear coupling the 
mixing angle tan (3 in the Higgs sector, and the sign of the higgsino mass parameter /i. 

In the CMSSM with the gravitino LSP, the next-to-lightest SUSY particle is either the 
lightest neutralino x5 or the lighter stau ri.^ The BBN constraints on electromagnetic 

"A stop ti NLSP is not feasible in the CMSSM [155]. 
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and hadronic energy injection disfavor the Xi NLSP for > 100 MeV [141, 156, 
137]. For the slepton NLSP case, the BBN constraints associated with hadronic/electro- 
magnetic energy injection have also been estimated and found to be much weaker but 
still significant in much of the parameter space [141, 156, 137, 139]. In the following, 
however, we shall see that the novel CBBN constraints drastically change this picture.^ 

Let us give an overview of what is presented in the remainder of this chapter. In 
Sec. 4.2.1 we consider the CBBN restrictions on the ri lifetime. Employing the results 
on m^-^ from a renormalization group analysis this allows us to relate the stau mass with 
the high-scale parameter mi/2- The obtained lower limit on the latter parameter can 
be translated into an upper bound on the reheating temperature. This will be done 
in Sec. 4.2.2. Both limits will be based on the estimate (4.3) of the stau decoupling 
yield. We contrast the obtained semi-analytical limits on with exemplary 

CMSSM parameter scans in Sec. 4.2.3. For those examples we will also explicitly consider 
in Sec. 4.2.4 the possibility of a non-standard cosmological history. In this context, we 
also check on the viability of thermal leptogenesis in Sec. 4.2.5. 

4.2.1 Lower limit on mi/2 

In the previous section we have realized that the new bounds emerging from the thermal 
catalysis of nuclear reactions yield the most dominant restrictions on the gravitino-stau 
parameter space. We have also noted that the discussion of those bounds is facilitated by 
the fact that both, the constraint on ^Li as well as the one on ^Be production, essentially 
disfavor the same region in the (mg, m;jr^ ) parameter space. In the following we focus on 
^Li and adopt as an upper limit on its primordial abundance 

^Li/Hlp < 6 X 10-1^ (4.7) 

We thereby resort in this work to a more conservative point of view than in our main 
discussions of the published works [1, 2, 3] which were based on ^Li/H|p < 2 x lO^^^ [39]. 
In this way, we give account to the concerns cast in the previous section (without, 
however, going to the very extreme using 10"^".) 

The limit (4.7) is shown by the dash-dotted (red) line in Fig. 4.1 as a function of 
the yield of negatively charged staus Y-^'^ and r^:^. As in Fig. 2.8 the curve is obtained 
by numerical integration of the Boltzmann equations presented in Chapter 2. In addi- 
tion, the solid lines show contours of constant ^Li output for other values than (4.7) as 

^Seeking a solution to the lithium problems, bound state effects within the framework of the CMSSM 
have first been considered in [145]. We set a different focus: We constrain the CMSSM parameter space 
from which we mainly derive an upper bound on Tr. 
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Figure 4.2: The solid lines are contours of constant ^Li/H as labeled and produced in 
CBBN. They are obtained by solving the full Boltzmann equations presented in Chap- 
ter 2. In addition, the dot-dashed (red) line shows the adopted limiting primordial 
abundance (4.7) with the region above being disfavored from ^Li overproduction. The 
dashed lines show the ^Li output of the Boltzmann network when the Saha approxima- 
tion (2.22) for the CHe-rf) bound state fraction is used. 

labeled. Moreover, the dashed gray lines show the ^Li contours when the Saha approxi- 
mation (2.22) is used for the bound-state abundance. The associated overestimation of 
^Li once more demonstrates the importance of a full numerical solution of the Boltz- 
mann equations. We remark that Fig. 4.1 contains an improvement with respect to 
the corresponding figure in [2] in the sense that proton-burning (2.36b) of ^Li has now 
been included and which leads to some reduction of the final ^Li output. Moreover, we 
now use a (^Herj") recombination cross section which includes the finite charge radius 
correction and accounts for recombination into IS and 2S states; see Chapter 2. 

Using the estimate (4.3) with Y-^'^ = Y~'^^/2 we find from Fig. 4.1 that the amount 
of ^Li produced in CBBN can be in agreement with (4.7) only for stau lifetimes of 

T^,<6x 103 s. (4.8) 

As can be seen from the supergravity prediction (4.4) of r^^ , the requirement (4.8) implies 
a lower limit on the splitting between and provided m^-^ < 0(1 TeV). Because of 
this hierarchy, the factor (1 — m~/m~^)-^ can be neglected in Eq. (4.4) in the following. 

Let us now turn to the CMSSM. We employ the computer program SPheno 2.2.3 [157] 
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to obtain the low-energy supersymmetric particle spectrum from the high-scale input at 
Mgut- In the region in which ti is the NLSP, we find 

< 0.21m^/2 (4.9) 
by scanning over the following parameter range: 

mi/2 = 0.1 ^6 TeV, 
tan ^ = 2 ^ 60, 
sgn^ = ±1, 
—AniQ < Aq < 4mo, 

and with mo as large as viable for a ri NLSP.^ 

For small left-right mixing, ri ~ tr, (4.9) can be understood qualitatively from the 
estimate for the mass of the right-handed stau near the electroweak scale [158] 

~ 0.15m^/2 + "^0 - sin^ Gw^z cos 2/3 . (4.10) 

since uiq <C rny2 ^ large part of the ti NLSP region. In fact, (4.9) tends to be 
saturated for larger rriQ, i.e., in the stau-neutralino-coannihilation region where the mass 
of the lightest neutralino m~o ~ m^, . This can be understood since the neutralino is 
bino-like in this region and typically m~o — 0.19m^^2-'' the remaining part of the stau 
NLSP region, smaller values of m^^^ satisfying, e.g., m~ = 0.15m^^2 can easily be found. 

To be on the conservative side, we set the stau NLSP mass rrir^ to its maximum value 
at which (4.9) is saturated: m~ = O.llrriy^. Using (4.4) this allows us to extract a 
lower limit on the universal gaugino mass parameter from the constraint (4.8) 

'"'/2>0-87TeV(^)''\ (4.11) 

Since for a ri NLSP typically tuq <C my 2^ it is the gaugino mass parameter mi/2 
which sets the scale for the low energy superparticle spectrum. Thus, depending on mg, 
the bound (4.11) implies rather high values of the superparticle masses. This is particu- 
larly true for the masses of the squarks and the gluino since their renormalization group 
running from Mqut to Q ~ 0(1 TeV) is dominated by M3{Q) ~ m;^/2Q^s(<3)/as(^GUT)- 
Therefore, for > 10 GeV, the cosmologically favored region is associated with a mass 
range that will be very difficult to probe at the Large Hadron Collider. 

^We choose mt = 172.5 GeV for the top quark mass. In addition, we use the Standard Model 
parameters mb(mb)^ = 4.2 GeV, af^{mz) = 0.1172, and Q-^^(mz) = 127.932. 

^This estimate is relatively independent of tan/3 and valid in the mi/2 region in which also the LEP 
bound on the Higgs mass [32], rrih > 114.4 GeV, is respected. 
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We stress that the scan over the entire natural CMSSM parameter space has enabled 
us to set a bound on mi /2 which depends on the gravitino mass but is independent of the 
CMSSM parameters.^ We, however, also remark that we have used the estimate (4.3) 
without accounting for the dependence of Y-^'^ on the SUSY parameters (other than 
rrir^). It will be thus important to reflect the bound (4.11) on exemplary scenarios where 
Y-^'^ is computed in each point of the parameter space. This will be done in Sec. 4.2.3. 
Moreover, in Part III we will carry out a detailed study of the dependence of the stau 
decoupling yield on the SUSY parameters. There, we will indeed find that it is possible 
to evade (4.11) in exceptional cases. 

Finally, we also note that in the derivation of (4.11) we have only made use of the 
catalytic BBN effects. In Sec. 4.1 we have seen that only the D constraint on hadronic en- 
ergy release can compete with the CBBN constraints. Accordingly, the D constraint can 
only tighten the bounds on (and Tr in the following). Thus, taking a conservative 
point of view, we are allowed to neglect this complication. 



4.2.2 Upper bound on Tr 

The amount of gravitinos produced in thermal scattering is sensitive to the reheating 
temperature Tr, and to the masses of the gauginos and hence to [119]. The associ- 
ated gravitino density can be approximated by^ 

This follows from (3.6). Here we use that the running gaugino masses Mj associated 
with the gauge groups SU(3)c, SU(2)l, and U(1)y satisfy M3 : M2 : Mi ~ 3 : 1.6 : 1 
at a representative scale of 10^ GeV at which we also evaluate the respective gauge 
couplings. Furthermore, we only need to take into account the production of the spin- 
1/2 components of the gravitino since (4.11) implies Mf /"im?- ^ 1 for mp; > 1 GeV. 

For a given the reheating temperature Tr is limited from above because 

cannot exceed the dark matter density (3.14). Using the derived lower bound (4.11) 
allows us to extract the upper limit: 

Ta<5xl0'GeV(^)"\ (4.13) 

This constraint is a slowly varying function of mg: (mg/10 GeV)^/^ = 0.6 2.5 for 

* Similar limits have also been discussed in models in which the ratio is bounded from 

below [159]. 

^For a discussion on the definition of Tr see the discussion in Sec. 3.3.1; here, ^ = 1.8. 
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mg = 1 GeV -7- ITeV. Therefore, (4.13) poses a strong bound on Tr for the natural 
gravitino LSP mass range in gravity-mediated supersymmetry breaking scenarios.^'' 

Note that the constraint (4.13) rehes on thermal gravitino production only. In ad- 
dition, gravitinos are produced in stau NLSP decays with the respective density (4.2). 
While the precise value of Y-^'^ depends on the concrete choice of the CMSSM parame- 
ters, the upper limit (4.13) can only become more stringent by taking 0-"'"^ into account 
(provided that (4.3) is not substantially depleted.) 



4.2.3 Exemplary parameter scans in the CMSSM 

Taking into account gravitinos from thermal production and from late decays of the 
lightest Standard Model superpartner we can confront our numerical findings with the 
above derived semi-analytical limits on and Tr for various values of mg. Con- 

sidering concrete CMSSM scenarios allows us to compute the thermal and non-thermal 
gravitino production in each point of the parameter space without relying on typical 
values of the decoupling yield Ynlsp of the NLSP such as (4.1). 

Earlier studies of Tr constraints within the CMSSM used the result of [118] to explore 
the viability of Tr > 10^ GeV [156, 137]. Our study presents also scans for Tr as low 
as 10^ GeV based on (3.6) which includes electroweak contributions to thermal gravitino 
production [108].^^ 

In Fig. 4.2 the solid (black) and dotted (blue) lines show respectively contours of 
1nlsp(2o) and rriNLSP in the (771^/2) ""t-o) plane for Aq = 0, fi > 0, tan/3 = 10 (left 
panel) and tan/? = 30 (right panel). Above (below) the dashed line, m~o < 
(m^r < rn~o). The medium gray and the light gray regions at small mi/2 are excluded 
respectively by the mass bounds m~± > 94 GeV and m/jO > 114.4 GeV from chargino 
and Higgs searches at LEP [32] . The leftmost dotted (blue) line indicates the LEP bound 
rrir > 81.9 GeV [32]. For tan/3 = 30, tachyonic sfermions occur in the low-energy spec- 
trum at points in the white corner labeled as "tachyonic." For those scans we employ 
the computer program SuSpect 2.34 [161] to calculate the low-energy spectrum of the 
superparticles and the Higgs bosons. ^'^ Assuming a standard cosmological history, the 
yield iNLSp(2o) is obtained from the I^nlsp^^ values provided by the computer program 
micrOMEGAs 1.3.7 [162, 163]. 

^"Similar, but less restrictive limits have been obtained in [140] by relaxing the CMSSM-specific 
splitting (4.9). 

^^Meanwhile, after pubhcation of [1], related works appeared [42, 160]; cf also Sec. 3.4. 
^■^In this section, we have used the following values: mt = 172.5 GeV, mb(mb)^^ = 4.23 GeV, 
af^{mz) = 0.1172, and Q-^^(mz) = 127.90896. 
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The contours shown in Fig. 4.2 are independent of and Tr. Therefore, they can 
be used to interpret the results shown in the figures below. Note the sensitivity of both 
l?;(To) and on tan/3. By going from tan/3 = 10 to tan/? = 30, Y^{Tq) decreases by 
about a factor of two at points that are not in the vicinity of the dashed line, i.e., that 
are outside of the r-x? coannihilation region. While m^f becomes somewhat smaller by 
increasing tan (5 to 30, the tan (3 dependence of mr;o is negligible. 

Let us now explore the parameter space in which the relic gravitino density matches 
the observed dark matter density (3.14), 

0.075 < ^ij^h^ + < 0.126 . (4.14) 

Now, Tr and appear in addition to the traditional CMSSM parameters. We focus 
on > 1 GeV since the soft SUSY breaking parameters of the CMSSM are usually 
assumed to result from gravity-mediated SUSY breaking. However, we do not restrict 
our study to fixed relations between and the soft SUSY breaking parameters such 
as the ones suggested, for example, by the Polonyi model. 

In Fig. 4.4 the light, medium, and dark shaded (green) bands show the ("ii/2)^o) 
regions that satisfy the upper limit (4.14) for Tr = 10^, 10^ and 10^ GeV, respectively, 
where tan (3 = 10, = 0, /i > 0. The four panels are obtained for the choices (a) = 
10 GeV, (b) = 100 GeV, (c) = 0.2 mo, and (d) = tuq. In the dark-gray 
region, the gravitino is not the LSP. The regions excluded by the chargino and Higgs 
mass bounds and the line indicating m~o = m~ are identical to the ones shown in the 
upper panel of Fig. 4.2. The dotted lines show contours of the NLSP lifetime (4.4). For 
the Xi NLSP, we calculate r~o from the expressions given in Sec. IIC of Ref. [141]. 

The tnlsp contours in Fig. 4.4 illustrate that the NLSP decays during/after BBN. 
Successful BBN predictions therefore imply cosmological constraints on mg, m-NLSP, 
and Inlsp [141, 156, 137, 139]. Indeed, as stressed before, it has been found that the 
considered x? NLSP region is completely disfavored for > 100 MeV by constraints 
from late electromagnetic and hadronic energy injection [141, 156, 137, 145]. In the 
Ti NLSP region, the constraints from electromagnetic and hadronic energy release are 
important but far less severe than in the x? NLSP case. 

Including the constraints from the bound-state effects this picture changes. As we 
have already seen in the previous section, in most of the r NLSP parameter space, the 
bounds from the catalysis of ^Li and ^Be can be much more severe than the ones from 
late energy injection. We incorporate the {t^^ , Y~- )-dependent CBBN constraint on ^Li 
from Fig. 4.1 which is shown in Fig. 4.4 by the long dash-dotted (red) line. The ti-NLSP 
parameter space to the left of this line is excluded. 
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Figure 4.3: Contours of 1nlsp(7o) (solid black lines) and mNLSP (dotted blue lines) 
in the (toi/2,'71o) plane for Aq = 0, ^ > 0, tan/3 = 10 (upper panel) and tan/3 ~ 30 
(lower panel). Above (below) the dashed line, m^o < rrir {rrir < ^^l^- The medium 
gray and the light gray regions at small mi/2 show the mass bounds m-± > 94 GeV 
and TTiho > 114: A GeV from chargino and Higgs searches at LEP [32]. 



88 Chapter 4. The stau as the NLSP 



= lOGeV, tan/3 = 10, = 0, /x > = 100 GeV, tan/3 = 10, = 0, > 




(c) mi/2[GeV] (d) mj/a [GeV] 



Figure 4.4: The (toi/25 "^o) planes for tan/3 = 10, Aq — 0, fj, > 0, and the choices 
(a) rriQ = 10 GeV, (b) — 100 GeV, (c) mg = 0.2 mo, and (d) = mo. In each 
panel, the light, medium, and dark shaded (green) bands indicate the regions in which 
n^h'^ e ildm'*^ for = 10^7 10*7 and 10*^ GeV, respectively. The medium gray and 
the light gray regions at small mi/2 are excluded respectively by chargino and Higgs 
searches at LEP. In the dark gray region, the gravitino is not the LSP. The dotted 

lines show contours of the NLSP lifetime. Below the dashed line, my < m^o. With 

^1 

the T NLSP, the region to the left of the long-dash-dotted (red) line is cosmologically 
disfavored by bound-state effects on the primordial ''Li abundance [79]. The effects of 
late hadronic energy injection on the primordial D abundance [139] disfavor the r NLSP 
region between the short-dash-dotted (blue) lines in panel (b) and the one above the 
corresponding lines in panels (c) and (d). In addition, in (d) the electromagnetic ^He/D 
constraint and the gluino mass contour rrig — 2.5 TeV are shown as labeled. 
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We thereby update our figures presented in the pubhshed work [1]. In those figures, 
the ^Li output was taken from the work [79] which provided the initial estimate on the 
efficiency of the catalyzed production. Meanwhile, the dedicated quantum-three-body 
calculation [84] became available which lead to a reduction of the S'-factor of (2.37) 
by roughly one order in magnitude. As discussed in Chapter 2, we have incorporated 
this state-of-the-art result together with other improvements in our Boltzmann network 
equation. 

In the lifetime region t^^ < few x 10^ s which is unconstrained by CBBN only the 
constraint from deuterium on the hadronic energy release becomes important. Following 
the approach explained in Sec. 4.1, we incorporate the constraint on hadronic energy 
release for D. In Fig. 4.4 these are shown by the short-dash-dotted (blue) lines. The 
D constraint disfavors the region between the corresponding lines in panel (b) and the 
region above the corresponding lines in panels (c) and (d). In panel (a) the D constraint 
does not appear. In addition, for orientation, in panel (d) we also include the elec- 
tromagnetic ^He/D constraint as a thin (pink) line and show the gluino mass contour 
TUg = 2.5 TeV which is a near to vertical thick (violet) line. 

Indeed, one finds in each panel of Fig. 4.4 that the highest Tr value allowed by the 
considered BEN constraints is about 10'' GeV. The bands obtained for Tr > 10® GeV are 
located completely within the region disfavored by the ^Li bound. In previous gravitino 
dark matter studies within the CMSSM that did not take into account bound-state effects 
on the primordial ^Li abundance, much higher temperatures of up to about 10^ GeV 
were believed to be allowed [156, 137, 119]. 

The constraint Tr < lO"^ GeV remains if we consider larger values of tan /?. This 
is demonstrated in Fig. 4.4 for tan/3 = 30, Aq = 0, fi > 0, and (a) = 10 GeV 
and (b) = tjiq. The shadings (colors) and line styles are identical to the ones in 
Fig. 4.4. 

Let us comment on the dependence of the considered BBN constraints on the assumed 
primordial abundances of D and ^Li. As can be seen in Figs. 4.4 and 4.4, the constraint 
from late hadronic energy release is quite sensitive on the assumed primordial D abun- 
dance. In contrast, even if we relax the restrictive ®Li bound on 5^lsp/2 by two orders 
of magnitude, we still find Tr < 10^ GeV. For example, the ^Li constraint relaxed in this 
way would appear in Fig. 4.4 (b) as an almost vertical line slightly above mi/2 = 2.5 TeV. 

Limit on mi/2 The limit (4.11) emerges since 171^^ scales with mi/2 [see Fig. 4.2] and 
since r:^^ is fixed once mg and m^^^ are specified. The choice = 10 GeV in 
Fig. 4.4 (a) and Fig. 4.4 (a) allows for an immediate comparison of the exemplary 
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rrig = lOGeV, tan/3 = 30, = 0, > = mo, tan = 30, = 0, > 




(a) mi/2[GeV] (b) mi/2 [GeV] 

Figure 4.5: CMSSM planes as in Fig. 4.4, but for tan/3 = 30, = 0, ^ > 0. In (a) 

= 10 GeV and in (b) = mo has been chosen. 



CMSSM scenarios with (4.11). Only in the vicinity of the dashed line, i.e., in 
the Ti-Xi coannihilation region, the position of the ^Li constraint approaches its 
conservative lower limit on mii2- This is because ri becomes heavier for larger mg 
which shortens r^^^ for fixed mg. Contrariwise, the splitting between the actual 
position of the ^Li constraint and (4.11) is larger for smaller mo- This is slightly 
more pronounced in Fig. 4.4 (a) than in Fig. 4.4 (a) and results from the fact 
that the increase in tan [i leads to a decrease in m^^ so that r^^ becomes larger for 
fixed mg. 

That the lower limit (4.11) can imply high values of the superparticle masses is 
illustrated by the vertical (violet) lines in Figs. 4.4 (d) and 4.4 (a) which show the 
gluino mass contour m-g = 2.5 TeV. In this regard, also note that the mass of the 
lighter stop is mj^ ~ 0.7mg in those ti NLSP regions with mh > 114.4 GeV. Since 
there the gaugino mass parameter sets the scale for the low energy superparticle 
spectrum, depending on m^, the bound (4.11) implies high values of the super- 
particle masses which can be associated with a mass range that will be difficult to 
probe at the LHC. 

Limit on Tr The limit (4.13) relies on thermal gravitino production only, 17-^ ~ Tr. 
Thus the upper limit on Tr becomes more stringent by taking $7-"'"^ into account. 

(jr 

In Figs. 4.4 (d) and 4.4 (b) we have fixed mg = niQ. Thereby, the non-thermal 
production (4.2) becomes more important for larger values of niQ. In addition, 
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Y-^'^ takes on its maximum at a given mx/2 in the t~Xi coannihilation region. 
This leads to the bending of the bands (4.14) towards lower mi/2- This indicates 
that (4.13) indeed seems to provide a good estimate. 

While the constraint < 10^ GeV is found for each of the considered rela- 
tions, one cannot use the ^Li bound to set bounds on m^^ without insights into 
mg. The ^Li (and similarly the ^Be) bound disappears for t:^ < 10^ s [79] which is 
possible even for m^-^ = 0(100 GeV) provided mg is sufficiently small; see (4.4). 
However, the constraints on Tr become more severe towards small as is shown 
in Fig. 3.2. Thus, the constraint Tr < 10'' GeV cannot be evaded by lowering 
TTig provided < Tp. An upper limit on of 10^ GeV can be problematic for 
inflation models and baryogenesis scenarios. This finding can thus be important 
for our understanding of the thermal history of the Universe. 

4.2.4 Late-time entropy production 

The constraints shown above are applicable for a standard thermal history during the 
radiation-dominated epoch. Such a standard cosmological evolution may, e.g., be ac- 
complished by considering only a minimal framework such as the MSSM.^^ 

On the other hand, focusing on a gravitino LSP we also explicitly consider the gravity 
sector. For example, it is well known that supergravity and string theories generically 
suffer from the appearance of scalar fields which can give rise to "cosmological moduli 
problems" [166, 167, 168]. Typically, the interaction of such an exotic (eventually) 
massive field (p to the MSSM sector is suppressed by a high-energy scale such as Mqut 
or Mp. The problem arises because (p easily drops out (or never has reached) thermal 
equilibrium. 

For a massive particle species which is in equilibrium, the energy density becomes 
exponentially suppressed once the temperature drops below its mass. However, if <j) is 
frozen-out, p<^//Orad starts to scale as T^^ so that, eventually, the energy density in cj) 
dominates over the one in radiation. Thus, if (p lives sufficiently long, it is possible that 
a substantial amount of entropy is released in its out-of-equilibrium decay. Since an 
entropy release tends to erase any pre-existing quantity which itself has originated from 

^^Even inflation seems to be feasible within the MSSM with the inflaton being a gauge invariant 
combination of squark and slepton fields [164]. Strictly speaking, however, neutrinos do not obtain 
masses within this framework which are needed to explain the observed neutrino oscillations [165]. 

^^Other entropy production events after inflation are, e.g., Q-ball decays; cf. [169] and references 
therein. Gravitino dark matter scenarios with late-time entropy production have been considered previ- 
ously for gauge-mediated SUSY breaking where Tr > [170, 171, 172, 173, 174]. 
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a non-equihbrium process — such as r/b — this can be problematic. Moreover, additional 
unwanted relics can be produced in the decay of cp leading, e.g., to too much dark matter. 
However, as we shall see in the following, entropy production can also have a positive 
effect [175]. 

The change d{a^s) = dS = dQ/T in entropy per comoving volume is found from the 
"heat added" in the decay, dQ = a^T^p^dt [176]. The associated evolution of the entropy 
per comoving volume is thus described by 

dS _ T^p^a^ ^'>-2 X 1/3 



dt 



^5cffJ T^p^a^S-'/^ (4.15) 



We have already indicated by the chosen notation that the situation is somewhat similar 
to that at the end of inflation. Thus, we can solve (4.15) by simultaneously considering 
the Boltzmann equation for p^ in the form (3.11b) together with the Priedmann equation 
which governs the evolution of the scale factor a of the Universe; denotes the decay 
width of (j). 

From the discussion of Sec. 3.3.1 we know that the temperature after the decay can 
be expressed in terms of F^, 



T, 



after 



10 



5efr(T'after)7r2 



1/4 



Vr>Mp , (4.16) 



which satisfies F,^ = 3-?/(Taftor)- Note that primordial nucleosynthesis imposes a lower 
limit on this temperature. It mainly emerges from the fact that a re-thermalization of all 
the three neutrino species after the reheating process does not happen instantaneously. 
Insufficient thermalization affects the Hubble rate, thereby the n/p freeze-out value, and 
thus the ''He output in BBN [see Sec. 1.3]. The bounds derived in [177, 178, 179, 180] 
are in the range 

Tafter > 0.7-4 MeV . (4.17) 

In the upper panel of Fig. 4.5 we show the evolution of 5, a^p^, and a^Prad for 
two exemplary scenarios respecting (4.17). The scale factor a is normalized by^^ ai = 
a(10 GeV) = 1 GeV~^ and the temperature dependence of /leflf is taken into account as 
determined in [28]. For ^^(10 GeV) = 0.1 /Orad(10 GeV) and Tafter = 6 MeV, S increases 
by a factor of A = 100 as shown by the corresponding solid line. For P(/,(10 GeV) = 
8prad(10 GeV) and Tafter = 4.9 MeV, S increases by a factor of A = 10^ as shown by 
the corresponding dotted (blue) line. 



Giving a the dimension of a length makes the radial coordinate in the Friedmann- Robertson- Walker 
metric dimensionless. 
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Figure 4.6: Top: Evolution of S, a p^, and a piad as a function of T for the 
normalization ai = a(10 GeV) = 1 GeV^^. The solid lines are obtained for 
P4,{10 GeV) = 0.1prad(10 GeV) and Taftor = 6 MeV, the dotted (blue) hues for 
P0(1O GeV) = 8pi.ad(10 GeV) and Tafter = 4.9 MeV. Bottom: Evolution of the temper- 
ature for the scenario which is depicted by the solid lines in the left panel. 
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The lower panel of Fig. 4.5 shows the associated evolution of the photon temperature 
for the case A = 100. As long as the comoving entropy is conserved, T oc h^^^^a~^. This 
causes T to decrease slightly less slowly than during the quark-hadron transition at 
T ~ 200 MeV which we have assumed to be adiabatic [181]. As already encountered 
in Sec. 3.3.1, the temperature scales with a~^/^ during entropy production, i.e., during 
reheating, the temperature does not increase but rather drops less slowly [176]. Below 
Tafter the Universe again expands as a radiation-dominated one. 

In the following, we restrict our study of entropy production at late times, Tbefore — 
Tiow Tpj,, so that the thermal production of gravitinos is not affected. We assume 
that the production of gravitinos and NLSPs in the entropy producing event, such as 
the direct production in decays of (/>, is negligible. Moreover, in this section, we focus 
on scenarios in which the decoupling of the NLSP is not or at most marginally affected 
by entropy production, i.e., either Tr S> Tkftcr S> T^^^^ or p^ad ^ Pep for T > j^nlsp^ 
Note that the latter condition excludes the event shown in Fig. 4.5 with A = 10^; see, 
however. Sec. 4.2.5. Thus, the thermally produced gravitino yield and — in the case of 
entropy production after NLSP decoupling — also the non-thermally produced gravitino 
yield are diluted: 

SiTio 

•S* (Tafter) 



YaiT^iter) = 9^ Ig(Tiow) . (4.18) 



In the case of late-time entropy production before the decoupling of the NLSP, we 
parameterize this by writing 

lf^(ro) = ^r|P(TioJ. (4.19) 

In this case, Inlsp(To) ^-nd thereby fi-^^ and the BBN constraints remain unaffected. 

Conversely, in the case of late-time entropy production after the decoupling of the 
NLSP (and before BBN) both, ^^^(ro) and yNLSp(7b), are reduced: 



>NLSp(To) = ^yNLSp(Tiow). (4.20) 

Accordingly, 17-^ and 17-"'"^ become smaller and the BBN constraints can be relaxed. 

In Fig. 4.7 we show how late-time entropy production before (left) and after (right) 
NLSP decoupling affects the ^Li constraint and the region in which 0.075 < n^h^ < 0.126 



^^Note however, if kinematically allowed, a (f> field will also typically decay into SUSY particles [182, 
183, 122]. Thus, the constraints discussed below shall therefore be considered as conservative/optimistic 
bounds. 
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for Tr = 10^ GeV. The (mi/2) ""t-o) planes are considered for tan/3 = 10, j4o = 0, ^ > 0, 
= 100 GeV (upper panels) and = mo (lower panels). The dark shaded (dark 
green) region is obtained without late time entropy production, 6 = A = 1. The medium 
and light shaded (medium and light green) bands are obtained with a dilution of 
(qTP _^ qNTP) by (5 = 10 (A = 10) and 6 = 100 (A = 100), respectively. The dot-dashed 

(j (j 

(red) line illustrates that the ^Li bound is independent of 6, as shown in the panels on the 
left-hand side, and becomes weaker (i.e., moves to the left) with increasing A, as shown 
in the panels on the right-hand side. Other curves and regions are identical to the ones in 
the corresponding panels of Fig. 4.4. Note that we do not show the D constraint on late 
hadronic energy injection since it is not sensitive to 6 and vanishes already for A = 10; 
an exception is the severe D constraint which still appears for A = 10 in panel (a). BBN 
constraints on Xi NLSP scenarios with entropy production after NLSP decoupling will 
be studied elsewhere. 

Comparing panels (b) and (d) of Fig. 4.4 with panels (a) and (c) in Fig. 4.7, we find 
that a dilution factor of 6 = 10 (100) relaxes the Tr bound by a factor of 10 (100). Since 
the BBN constraints are unaffected by 6, the cosmologically disfavored range of NLSP 
masses cannot be relaxed. With the dilution after NLSP decoupling, the relaxation of 
the Tr constraints is more pronounced. Here also the cosmologically disfavored range 
of NLSP masses can be relaxed [175]. However, as can be seen in panels (b) and (d) of 
Fig. 4.7, the ^Li bound is persistent. With a dilution factor of A = 100, large regions 
of the {1711/2,^0) plane remain cosmologically disfavored. For even larger factors of A, 
however, the ^Li bound can be evaded as will be shown explicitly below. 

Figure 4.7 shows that inflation models predicting, for example, Tr = 10^ GeV become 
allowed in the CMSSM with gravitino dark matter for 5 = A ~ 100. Here it is not 
necessary to have late-time entropy production in the somewhat narrow window between 
NLSP decoupling and BBN. This is different for the viability of thermal leptogenesis in 
the considered scenarios where Tf' > Tr. 

4.2.5 Viability of thermal leptogenesis 

The constraint Tr, < 10^ GeV obtained in the considered CMSSM scenarios for a stan- 
dard cosmological history strongly disfavors thermal leptogenesis. However, if entropy 
is released after NLSP decoupling, a dilution factor of A > 10^ can render thermal 
leptogenesis viable for Tr > 10^^ GeV. 

Thermal leptogenesis — in its simplest form, with hierarchical heavy right-handed Ma- 
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Figure 4.7: The effect of late-time entropy production before (left) and after (right) 
NLSP decouphng on regions in which 0.075 < n^K^ < 0.126 for Tr, = 10^ GeV. The 
(mi/2, ^o) plane is shown for tan/3 = 10, Aq = 0, jj, > 0, ttIq = 100 GeV (upper panels) 
and ttIq = mo (lower panels). The dark shaded (dark green) region is obtained without 
late-time entropy production S = A ~ 1. The medium and light shaded (medium and 
light green) bands are obtained with a dilution of nj^^ (f]TP+f]^TP) by (5 = 10 (A = 10) 



and 6 = 100 (A = 100), respectively. The r NLSP region to the right of the dot-dashed 
(red) line is cosmologically disfavored by the primordial ®Li abundance. Other curves 
and regions are identical to the ones in the corresponding panels of Fig. 4.4. The severe 
D constraint for A = 10 appears only in panel (d). 
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jorana neutrinos — usuahy requires Tr > 10^ GeV [134, 135, 184, 185]. However, late- 
time entropy production dilutes the baryon asymmetry which is generated well before 
NLSP decoupling, 

rjhiTa.fter) = ^ ??b(7before) • (4.21) 

Therefore, the baryon asymmetry before entropy production must be larger by a factor 
of A in order to compensate for the dilution. For A ~ 10^, this can be achieved in the 
case of hierarchical neutrinos for Mri ~ Tr > 10^^ GeV, as can be seen in Fig. 7 (a) of 
Ref. [186] and in Fig. 2 of Ref. [187]. Here, Mri is the mass of the lightest among the 
heavy right-handed Ma jorana neutrinos. 

That it is indeed possible to produce such a large amount of entropy in the narrow 
time window between NLSP decoupling an BBN, we show in Fig. 4.5 a scenario in which 
a dilution factor of even A = 10^ is generated in the out-of-equilibrium decay of a heavy 
particle cp [dotted (blue) lines]. However, because of p</,(10 GeV) = 8prad(10 GeV), the 
Hubble rate is already enhanced already during the decoupling phase of the NLSP. This 
leads to an increase of r/^LSP j^j^j yNLSp(7f ^^^)- In the results shown below, we indeed 
account for this by using a modified version of the micrOMEGAs code. Thereby, the Inlsp 
contours shown in Fig. 4.2 do not apply in this section. After entropy production, the 
net effect is still a significant reduction of Inlsp(^o) — provided that (p does not decay 
into SUSY particles. 

In Fig. 4.7 we consider now two concrete scenarios with P(/,(10 GeV) = 8/9rad(10 GeV) 
as initial condition and with A = 10'^ and 2 x 10'^ corresponding to the respective 
reheating temperatures Tafter = 48 MeV and 24 MeV. We choose the (mi/2) "io)-pla'ne 
with tan/3 = 10, Aq = 0, fj, > 0, and = ttt-q. Here the shaded (green) bands indicate 
the region in which 0.075 < Q^h'^ < 0.126 for Tr = 2 x 10^^ Q^y ^nd A = 10^ (dark) and 
2 X 10^ (medium). In addition, the corresponding evolution of the ^Li bound is shown by 
the dot-dashed (red) lines. The regions in the ri-NLSP region above the curves remain 
cosmologically disfavored. The gray regions are identical to the ones in Fig. 4.4. 

We find that the ^Li bound can indeed be evaded for A = 2 x 10'^. Gravitino dark 
matter scenarios with successful thermal leptogenesis in the ri NLSP region are located 
on the light-shaded (light green) band. As can be seen in Fig. 4.7, the ri NLSP region 
with mi/2 ^ 700 GeV where m^^ < 250 GeV (cf. Fig. 4.2), is no longer disfavored by 
the ^Li bound provided A > 10'^. Such scenarios are particularly promising since the 
long-lived ri NLSP could provide striking signatures of gravitino dark matter at future 
colliders [154, 188, 189, 190, 191]. 

^^Right-handed (heavy) neutrinos (and their superpartners) are not part of the (C)MSSM. Again, this 
section contains an updated discussion of the pubhshed work [1]; see Sec. 4.2.3. 
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Figure 4.8: The effect of entropy production after NLSP decoupling for Tr = 2 x 
10^^ GeV and A > 10^ in the (mi/2,mo) plane for tan/3 = 10, Aq ^ 0, ^ > 0, and 
mg = mo. The shaded (green) bands show the region in which 0.075 < il^h^ < 0.126 
for A = 10^ (dark) and 2 x 10^ (medium). The dot-dashed (red) lines illustrate the 
corresponding evolution of the ^Li bound with the regions below being cosmologically 
allowed. 

Finally, let us remark that the exact value of the actual amount of entropy required 
is model dependent. Though we have accounted for the fact that the presence of the 
energy density during NLSP decoupling can affect ^nlsp ' have neglected a possible 
branching ratio of the (j) decays into ri and/or gravitinos. Moreover, such a scenario is — 
of course — fine-tuned since entropy has to be released in a very narrow time window. 
However, with our improved treatment of the CBBN yield of ^Li it becomes slightly easier 
(when compared to [1]) to circumvent the stringent ^Li bound. This is indicated by 
the fact that we were able to allow for an increased NLSP decoupling yield by choosing 
P(j,{W GeV) = 8prad(10 GeV) which lead to the generous choice of Tafter > 20 MeV. 
Since we have not provided explicitly a model for cf) we understand the results presented 
in this section as a proof-of-concept study that a non-standard cosmological evolution 
can evade even the most stringent of all BBN constraints while still allowing for successful 
leptogenesis. 



^^The value is now in the region provided by the more recent Ref. [84] in which also our employed 
CBBN rate for ^Li production was obtained. The authors, however, only use the Saha approximation 
for the (''Her]") bound state abundance; see Sec. 2.3. 
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Chapter 5 



Thermal relic stau abundances 



In Part I of this thesis we worked out the effects which a long-hved electricahy charged 
massive particle species has on the predictions of primordial nucleosynthesis. In 
Part II we then considered concrete realizations of in gravitino dark matter scenarios 
in which ri is the lightest Standard Model superpartner. There, a central parameter in 
the investigation is the abundance of ti during/after BBN. For a standard thermal 
history it is determined by t^^-^ and by its thermal freeze-out value which we shall call 
1~ = Y-^'^ in the following; recall that Y~^^ = [n-^'^ + n~T)/ s is the total stau number 
density normalized to the entropy density prior to its decay . 

In particular, we saw that (i) governs the non-thermally produced relic den- 
sity (4.2) of gravitino dark matter that originates from ti decays, (ii) controls the 
hadronic and electromagnetic BBN constraints by quantifying the total energy density 
m^-^sY:^ which is eventually be released in the decay and (iii) parameterizes the abun- 
dance 1^/2 of recombination partners with the light elements leading to the catalysis of 
BBN. 

Importantly, the points (i)-(iii) have in common that they give rise to upper limits 
on 1~: (i) For example, for = 50 GeV and a thermally produced gravitino density 
J^TP ^ 0.99 J^dm (0.9 0dm), one finds from (4.14) Y^ < 10'^^ (10"^^); see also Fig. 13 
of [139]. (ii) The BBN constraints on hadronic and electromagnetic energy release can 
be as restrictive as Y^ < 10"^^ (10^^^); cf. Fig. 12 of [139] and Figs. 14 and 15 of [130]. 
(iii) Catalyzed production of ^Be (and ^Li) imposes restrictive upper limits of IV ^ 
2 X 10-^^ (2 X 10-^^-2 X 10"^^) for t^, > 10^ s; see Fig. 2.8. 

In Chapter 4 we have either made use of representative values (4.1) of 1^ or performed 
exemplary CMSSM parameter scans which gave us 1^, e.g., for fixed values of tan /3 and 
At- In this chapter we now calculate Y^ by taking into account the complete set of 
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stau annihilation channels in the MSSM with real parameters for SUSY spectra for 
which sparticle coannihilation is negligible. Using our own code for the computation 
of the resulting thermal relic stau abundance Y^, we examine explicitly (i) the effect of 
left-right mixing of the lighter stau, (ii) the effect of large stau-Higgs couplings, and 
(iii) stau annihilation at the resonance of the heavy CP-even Higgs boson H^. We 
consider both the "phenomenological MSSM" (pMSSM) (see, e.g., [161]) in which the 
soft SUSY breaking parameters can be set at the weak scale, and the CMSSM. 

Within the framework of the pMSSM, we show examples in which 1^ can be well 
below 10~^^. Even within the CMSSM, we encounter regions with exceptionally small 
values of 5^ < 2 x 10^^^. We stress that the results in the following are independent 
on the nature of the LSP. However, we discuss the implications of these findings for the 
gravitino-stau scenario.^ We also address the viability of a ri-rf asymmetry. The key 
quantities for the significant reduction could be probed at both the LHC and the 
ILC. 

The work presented in this chapter is based on the publication [5]. A calculation of 
the thermal relic abundance of long-lived staus has also been part of a detailed study [36] 
which focuses on gauge interactions and on the effect of Sommerfeld enhancement. In 
contrast, the most striking findings of our study — in which Sommerfeld enhancement is 
not taken into account — are related to the Higgs sector of the MSSM. 

We also remark that our work has some overlap with [192] which appeared as [5] was 
being finalized. Whereas [192] focuses on the potential suppression in 1~ due to enhanced 
annihilation into the lighter Higgs final state, our work provides an investigation of stau 
decoupling based on a complete set of annihilation channels. Thereby, we also study 
enhanced stau annihilation into heavier Higgs final states and consider annihilation at 
the heavy Higgs resonance. 

The outline of this chapter is as follows. In the next section we review basic proper- 
ties of the staus to introduce our notations and conventions for the stau mixing angle. 
Section 5.2 explains the way in which we calculate Y^ and provides the complete list 
of stau annihilation channels. In Sec. 5.3 we analyze the dependence of the most rel- 
evant stau annihilation channels on the stau mixing angle. Effects of large stau-Higgs 
couplings and stau annihilation at the resonance are studied in Sects. 5.4 and 5.5, 
respectively. The viability of a ri-rj* asymmetry is addressed in Sec. 5.6. In Sec. 5.7 we 
present exemplary parameter scans within the CMSSM that exhibit exceptionally small 
Y^ values. Potential collider phenomenology of the parameter regions associated with 
those exceptional relic abundances and potential implications for gravitino dark matter 

^For simplicity, we call the parameter region in which m~o < m^^^ the x?-NLSP region in the following. 
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scenarios are discussed in Sects. 5.8 and 5.9, respectively. 



5.1 Stau mixing and mass eigenstates 

In this section we review some basic properties of the stau to set the notation. In 
absence of inter-generational mixing, the stau mass-squared matrix in the basis of the 
gauge eigenstates (tl,tr,) reads 

9 fm'i + injr rrirX* \ „ , , , 

-^1 = ^/"^ , , = (RrV { , ] Rr (5.1) 



rrirXr mi. + m 




with 



m^L = "^II + (-^ + si^^ co*^ 2/? (5.2) 

TTiRR = - sin2 6'H/ M| cos 2/3 (5.3) 

Xr = Ar- n* tan /3 . (5.4) 

Here, m^^ and m^f:^ are the soft SUSY breaking masses, Aj- is the trilinear coupling, ^ 
is the Higgs-higgsino mass parameter, and tan/3 = V2/V1 denotes the ratio of the two 
Higgs vacuum expectation values. In this work we restrict ourselves to the MSSM with 
real parameters. Then X* = X^ so that the mass eigenstates ri and T2 are related to tl 
and tr by means of an orthogonal transformation 

with i?^ = . ; M (5.5) 

with % denoting the stau mixing angle. Imposing the mass ordering m^^ < 
choosing < < vr, the mixing angle can be inferred from the elements of 

2mT-XT- 2mT-XT 2mT-Xr 

ian.2dr = — , 7y— = , sm.2d:f = — 7^ ^ , (5.6) 

mf T — m?>r, m~ — m~ 

LL ritx Tl T2 

where the sign of the second relation determines the quadrant of O^f. In the first relation, 
we have introduced 6 = m^L — nT-jiR- particular, % = 7r/2 corresponds to a purely 
right-handed stau, ti = tr, whereas maximal mixing occurs for 6^ = tt/A and 37r/4. The 
physical stau masses are then given by 




1 



2 2,2, 



5^ y/W+^W^ (5.7) 

from which we see that an increase of \Xt-\ leads to a reduction of m^^. 



5.2. Calculation of the thermal relic stau abundance 103 



5.2 Calculation of the thermal relic stau abundance 

We have undertaken the effort to set up our own full-fledged relic abundance calculation. 
Let us in the following give a description of our approach to compute the stau yield 1?^. 
Throughout this work we assume a standard cosmological history with a temperature T 
of the primordial plasma above the stau decoupling temperature Tf so that the lighter 
stau Ti was once in thermal equilibrium. Then, the total stau yield 1~ = 1-^+1-* is found 
by solving the Boltzmann equation (1.9) with Yx = Yt Using the Maxwell-Boltzmann 
approximation, the stau equilibrium yield Y-'^ is given by 

and the thermally averaged annihilation cross section by [28] 

(<Tv){T) = , ^ \ — - r ds V^Ki ( P'is<s) , (5.9) 

where Ki is the modified Bessel function of order i and Pcff = s — Am-^ j 2. 

Note that (crt>) contains all the information from the particle physics side. It is 
obtained by computing the total stau-annihilation cross section, 

^ = \ ^tot with (Ttot = C^n n^rr + X] '^n r*^x , (5-10) 

X 

where the sum for the annihilation of t\ pairs^ has to be taken over all final states X. 
The factor 1/2 is convention but gives (1.9) its familiar form. The complete list of annihi- 
lation processes in the MSSM with real parameters — save for coannihilation processes — is 
given in Table 5.0.^ In addition, this table shows all possible particle exchanges, where 
s, t, and u are the Mandelstam variables which denote the respective channel. A number 
of annihilation processes proceeds also via a four-point vertex. Those are marked in the 
column named "contact." Already by mere optical inspection, we immediately see that 
the Higgs sector plays potentially an important role in the determination of the stau 
yield Y^. 

For all channels in Table 5.0, we generate Fortran code for the squared matrix el- 
ements by using the computer algebra packages FeynArts 5.4 [194, 195] and 

^Counting wise we distinguish between ri ti — ^ X and the conjugate process ri ri ^ X. In absence 
of CP violation in the SUSY sector, their cross sections agree so that we can solve a single Boltzmann 
equation (1.9) for obtaining VV. 

•^For a purely right-handed stau ri — tr, the stau annihilation channels and associated cross sections 
have already been presented in Ref. [193] in the context of Xi-ri coannihilation. 
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Table 5.1: The complete set of stau annihilation channels in the MSSM with real 
parameters for scenarios in which sparticle coannihilations are negligible. The mass 
eigenstates of the Higgs fields are denoted by /i", H'^, v4°, and and the ones of the 
neutralinos, the charginos, and the tau sneutrino by Xi 4, xt2^ ^^'^ ^r, respectively. 
Because of the absence of a tiTiA^ coupling (cf. Sec. 5.4), s-channel exchange of the 
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FormCalc 5.3 [196, 197]. For a chosen point in the SUSY parameter space, we then 
compute the radiatively corrected superparticle spectrum by running the spectrum gen- 
erator SuSpect 2.40 [161]. Its output allows us to set all SUSY parameters so that we 
can compute the total cross section a"tot(s) given by (5.10) and subsequently the ther- 
mally averaged cross section (5.9). Numerically, the computation of (5.9) is the most 
demanding part in the relic abundance calculation. In particular, we take special care 
about the following cases: 

• iJ^-resonance: Resonant stau annihilation via exchange is one of the central 
points in this part. In the generation of the matrix elements, we have therefore 
included the total //"-width Tfjo in the respective s-channel propagators. 

• Propagator poles: A diverging t(u)-channel propagator can be encountered when a 
production threshold is met. We overcome this problem by including a "sparticle- 
width" of O.Olm^:^ in the respective propagators in the vicinity of dangerous thresh- 
olds. A particularly interesting example with a diverging t(u)-channel propagator 
is given by the process ri r* — > jH^ if ^/s = rrif^o is fulfilled since then the H^- 
exchange in the s-channels of other processes is resonant simultaneously. 

• Bessel functions: The Bessel functions in (5.8) and (5.9) exhibit an exponential 
behavior for large arguments x ^ 1 [68] 



For small temperatures T, the arguments of Ki and K2 in (5.9) become large 
simultaneously. Therefore, in order to ensure numerical stability, we expand the 
Bessel functions in (5.9) for rrir^/T > 35 as in (5.11) and cancel the exponents 
analytically. 

We find the starting point for the numerical integration of (1.9) by solving [163] 



where g*{T) is given by (1.10). Tfi marks the point at which the stau starts to de- 
couple chemically from the background plasma, ^^(Tfi) — Y~^{Tii) ~ Ay~'^(rfi) with 
A = 0.1 [163] chosen in our code. Since we use a globally adaptive Gaussian integration 
routine to calculate (5.9), the computation of {crv){T) is time-demanding. Therefore, 
we evaluate (5.9) on a grid of different temperatures and use cubic spline interpolation 
to obtain values in between. We then solve the Boltzmann equation (1.9) by numerical 




(5.11) 




(5.12) 
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Figure 5.1: Schematic overview of our approach to calculate Y^: (i) We generate ma- 
trix elements from which we obtain the differential cross section da /d cos shown here 
is an example of the tt channel, (ii) Integration yields the invariant cross section a. 
The effect of "thermal weighting" of a is shown in (iii) by plotting (5.9) before inte- 
gration for the single tt channel (as a function of Pcff-) Summing up all annihilation 
channels and integration of (5.9) yields (crv) in (iv) at a grid of temperatures (crosses), 
(v) Taking into account the temperature dependence of and hcs and — upon cubic 
spline interpolation of (crv) in (iv) — we arrive at iu (vi) from integration of (1.9) . 
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integration from Tfi to zero. There, we fully take into the account the temperature 
dependence of and h^s by interpolating the respective tabulated values provided as 
part of the relic density code DarkSUSY 4.00 [33]. The freeze out temperature can then 
be defined by Tf = (Tfi + Tf2)/2 where Tfa is given by Yp{Ti2) = Yr{Ti2)/W [163]. For 
T < T{2, residual annihilations will further reduce Y:^ so that we refer to the decoupling 
yield Y-'^'^ as the quantity at the endpoint of integration. As already pointed out in the 
introduction, for simplicity, we call this yield Y^. Moreover, we will quantify T in terms 
of x = m^f-^/T and in particular Tf in terms of Xf = m^f^/Ti. We have also schematically 
depicted the approach to Y^f in Fig. 5.0 for an exemplary scenario (which gives rise to 
resonant stau annihilation. See the figure caption for details.) 

Note that we have additionally modified the FeynArts MSSM model file for the 
generation of the matrix elements in two ways: The first version, which we use throughout 
Sects. 5.3-5.5, allows us to set all g^^fc^Higgs and all trilinear Higgs couplings by using 
the computer tool FeynHiggs 2.6.3 [198]; see also Sects. 5.4 and 5.5. The second 
version allows for a direct comparison with the existing computer code micrOMEGAs 
2.0.6 [162, 163, 199]. We have transcribed their routine [200] for the computation of 
the running quark masses to Fortran, adopted all ^feQ^ -Higgs couplings, and modified 
all Higgs-self couplings of our matrix elements to match with their implemented version 
of the MSSM [201]. Using this second version, we find perfect agreement between our 
codes. ^ 



5.3 Dependence of stau annihilation on the stau mixing 
angle 

In order to isolate the distinct features of the different annihilation processes we need 
to have full control over the superparticle mass spectrum. Therefore, in the following, 
we will not rely on any constrained model (such as the CMSSM) where the soft-SUSY 
breaking parameters are subject to stringent boundary conditions at some high scale 
(such as Mgut)- In those models, the mass spectrum is found only after renormaliza- 
tion group (RG) evolution from the high scale down to the electroweak scale. Instead, 
we choose to work in the framework of the "phenomenological MSSM" (pMSSM), see, 
e.g., [161]. There, all soft-SUSY breaking parameters can be set at the scale of elec- 

*For our computation we use the Standard Model parameters rut = 172.5 GeV, mb(mb)'^^ — 
4.25 GeV, aF^(Mz) = 0.1172, a-i^(Mz) = 127.932, and Mz = 91.187 GeV. Since micrOMEGAs 
has hard-coded sm9w = 0.481 from which it computes Afw using the on-shell relation with Mz, we 
follow their convention to allow for a better comparison of our results with micrDMEGAs. 
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troweak symmetry breaking — a low scale — which we fix to ~ 2mr^. In particular, one 
can also trade the Higgs mass-squared parameters m'jj^ and rn'jj^ against /u and the pseu- 
doscalar Higgs boson mass m^o.^ Choosing fi as an input parameter is very convenient 
for two reasons: First, together with the specification of the gaugino masses Mi^2 we 
have control over the gaugino/higgsino mixture of the neutralinos Xi- Second, fi enters 
directly into the stau-Higgs couplings, whose importance will become clear in the next 
section. Furthermore, in the following, we choose to set all soft-SUSY breaking scalar 
masses (apart from m^and m^f:^) to a common value Ms = 1 TeV. Thereby, we essen- 
tially decouple all sfermions which are not of interest for us. This ensures also that we 
never enter accidentally any coannihilation regime. Finally, for simplicity, we set also 
all trilinear parameters to a common value A. Given /j,, Aj- = A, and tan /?, and thereby 
Xr, we can then fix m^^^ and 9:^ to arbitrary values by adjusting m'^^ and 6 in Eqs. (5.6) 
and (5.7). 

In the following, we will focus on two distinct regions of the SUSY parameter space. 
In the beginning, we will choose m^o to be very large m^o = 1 TeV ^ Mz- This cor- 
responds to the decoupling limit of the MSSM where the following (tree-level) relations 
hold [202] 

mlo ~ M| cos^ 2/3, mjjo ^ M| sin^ 2/?, (5.13) 

9 9 9 9 , X i sin^ 4/3 

=mjio +M4, cos^ (/? - a) - ^ ^ ' . (5.14) 

AO 

Therefore, m^o ~ mji^o — "TJ-jji up to corrections O (^M^/rrij^o^ so that any of the stau 
annihilation channels into heavy Higgs bosons is kinematically blocked. Furthermore, 
cos {13 — a) = up to corrections O (M|/m^o) implies that the H^VV coupling {V = 
W) becomes very small so that we loose the ff'^-exchanges in the stau annihilation 
channels with a VV final state. At the same time, the light Higgs boson takes on its 
Standard Model value for the hPVV coupling. Complementary to that we will consider 
also regions of the SUSY parameter space with smaller m^o, e.g., in the next section, 
where we will put a stronger focus on the Higgs sector and its connection to Y^. 

In Fig. 5.1a we show the ^-dependence of the masses of the heavier stau, m-^^j, 
(curved solid line) and the tau-sneutrino, rriy^, (dashed line) for fixed m^^ = 130 GeV 
and the input parameters tan (3 = 10, m^o = ^ = —A = 1 TeV, and 6Mi = M2,3 = 
1 TeV. Because of SU(2) gauge invariance, sets also the soft-breaking mass for 
the tau-sneutrino hence approximately m~ ~ m-^ + 5 so that becomes lighter than 
Tl for 9^ < 18° {5 is negative in that region). In addition, we plot the masses of the 

^Though the advocated procedure may require fine-tuning in the electroweak symmetry breaking 
conditions, it conveniently provides us with running parameters at the scale of stau annihilation. 
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Figure 5.2: (a) The dependence of my^ (curved solid line) and m^;^ (dashed line) on 

the stau mixing angle % for the input parameters m^:^ = 130 GeV (horizontal solid 

line), tan/3 = 10, m^o = = Mg = = 1 TeV, and 6Mi = hh = /i = 1 TeV (i.e., 

Xi ~ B) for which m^o = 169 GeV (dash-dotted line) and to/jO = 116 GeV (dotted 

line), (b) Dominant stau annihilation cross sections times the relative velocity of the 

incoming staus as a function of for Pos = 10 GeV and the same input parameters 

as in (a). The curves show the channels with the following final states: h^h'^, 77, tt, 

WW, TT, ZZ, jZ (at 9r = 40°, from top to bottom). In addition, we plot (TttVi- 

for the case of a wino-like neutralino, Xi — with m^o = 175 GeV as obtained with 

Ml = 6M2 = 1 TeV (thin gray line). No lines are shown for 9^ < 18° where mp^ < rrir^ . 
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lightest neutralino, m~o = 169 GeV (dash-dotted line), the lighter stau, = 130 GeV 
(horizontal solid line), and the lightest Higgs, m^^o = 116 GeV (dotted line). We note in 
passing that m^^^ may deviate slightly from its anticipated input value due to radiative 
corrections. We then correct for this by an adjustment of so that we indeed ensure 
m^^ to be constant. 

In Fig. 5.1b we plot the dominant stau annihilation cross sections times the rela- 
tive (non-relativistic) velocity in the center-of-mass frame of the incoming staus, ?;r = 
2Pes/m^^, for the same parameters as in Fig. 5.1a. Owing to an (approximate) Maxwell- 
Boltzmann distribution of the stau velocity, (-Pes) |Tf ~ y^m^^^Tf, we choose Peff = 10 GeV 
as a representative value. ^ The curves show the annihilation channels with the following 
final states: hPhP, 77, rr, WW, tt, ZZ, 7Z (at 6*^ = 40°, from top to bottom). All 

^This value is actually at the somewhat lower end, given m^-^ > 100 GeV and Tf ~ mpi/25. However, 
aVr depends only weakly on PcH, and the thermally averaged (aiv) will be shown in the upcoming figures. 
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channels except 77 show a strong dependence on 6^. The h^h^ (rr) channel peaks at 
0~ = 7r/4 — a feature which we will discuss in detail in Sec. 5.4. For the rr channel, 
the overall size of the cross section is governed by m^o since this channel proceeds only 
via t(tt)-channel exchanges of neutralinos. Our chosen input values lead to a bino-like 
neutralino, Xi — and cTt-t drops for an increasingly 'left-handed' stau. (For com- 
parison, the thin gray line shows fTrr^r for the case of a wino-like lightest neutralino, 
Xi = W , of similar mass, m~o = 175 GeV, as obtained by changing the gaugino mass 
input parameters to Mi = 6M2 = 1 TeV.) The annihilation into a WW pair becomes 
important for an increasing tl component in ri, i.e., towards smaller 9^, since the t{u)- 
channel exchange with the tau-sneutrino opens up; the tiVt-W (tiTiWW) coupling is 
proportional to cos 9^ {cos^O^). The modulation of the 7Z channel can be understood 
by considering the structure of the nriZ coupling oc (1 — Asiv? 9w + cos 20:^). Note 
that the first two terms practically cancel out. For stau annihilation into a ZZ pair 
there is an additional contribution from r2-exchange with the respective rir2-^ coupling 
oc sin20:f. Having discussed the dominant ri annihilation channels in a simple manner, 
we also warn the reader that interferences between the different Feynman diagrams of a 
given channel may well lead to a counter-intuitive behavior. In this regard, see Ref. [36] 
for a thorough discussion of rirj* annihilation into vector bosons. For the limiting case of 
a purely 'right-handed' stau, ri ~ -Fr, — > 7r/2), we recover the relative importance of 
the annihilation cross sections into 77, 7Z, ZZ, and rr with bino t(ti)-channel exchange 
found in Ref. [31]. 

Figure 5.2 shows the ^^^-dependence of (upper panel) and of the relative importance 
of the dominant thermally averaged cross sections, {criv) / {atotv) , at 2; = 25 (lower panel) 
for the same input parameters as in Fig. 5.1. The lines in the lower panel are associated 
with the same dominant annihilation channels as in Fig. 5.1b. In addition, the relative 
importance of the sum of the displayed cross sections, (o'disp^')/(ctot^'), (thin line, as 
labeled) is shown to demonstrate that the displayed channels constitute indeed (up to at 
most about 10%) the dominant part of {atotv) for the chosen set of input parameters. In 
the upper panel, the total stau decoupling yield obtained with our own relic abundance 
calculation is shown by the thick line and the one computed with micrOMEGAs, Y~^^, by 
the thin gray line. For 9^ < 25°, both curves start to deviate from each other since one 
enters the Vt--ti coannihilation region in which the stau decoupling yield increases. This 
coannihilation effect leads also to the rise of the thin gray line that shows (dtot^) Y~^^ 
in arbitrary units (a.u.) in the lower panel. Note that the same line illustrates 1^ oc 
l/(c"tot'f) for > 25°, where the result of our relic abundance calculation agrees with 
1^^. Interestingly, for the given input parameters, Y^ is not overly affected by the 
variation in 0;jr in this region, which reflects the fact that (crtotv) and thereby {crv) vary 
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Figure 5.3: Dependence of the stau yield Yf (upper panel) and of the relative impor- 
tance of the dominant thermally averaged cross sections, {oiv) / {a totv), at x = 25 (lower 
panel) on the stau-mixing angle 9^ for the same input parameters as in Fig. 5.1. In the 
upper panel, the thick line shows the stau yield obtained with our relic abundance 
calculation and the thin gray line the one obtained with micrOMEGAs to which we refer 
as Fi"^. In the lower panel, the line styles are associated with the same dominant 
annihilation channels as in Fig. 5.1b. In addition, we show (as labeled) the relative im- 
portance of the sum of the displayed cross sections, (crdispt^)/(o'toti'), and {atotv) 1-"^ 
in arbitrary units (a.u.). No lines are shown for 9^ < 18° where mp^ < rriri- 
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by less than a factor of about 1.5 at the relevant time of decoupling. In the next sections, 
we will demonstrate that this picture changes significantly for certain other choices of 
the input parameters. 



5.4 Effects of large stau-Higgs couplings 

Owing to the scalar nature of the stau, there exists a remarkable difference between the 
standard neutralino decoupling and the scenario in which the long-lived stau freezes out 
from the primordial plasma. For the neutralino LSP, the ^ parameter enters into the 
annihilation cross sections only indirectly by influencing the gaugino/higgsino mixture 
of Xi- This stands in strong contrast to the case in which a scalar particle is the lightest 
Standard Model superpartner: the sfermions couple directly to dimensionful parameters 
of the theory, namely, the trilinear couplings A and the Higgs-higgsino mass parameter 
/I. The corresponding operators in the MSSM Lagrangian always contain a Higgs field. 
In particular, the stau-Higgs couplings are given by 

^MSSM 9 ^ T*C[T*,Tp,n]Tpn (5.15) 

^ a,l3=L,R 

with Tl = h^, H^, A^. We have pulled out the factor g/M-w so that the 'reduced' 
couplings C[T*,Tf^,7i] among the gauge eigenstates tl and tr are given by [203] 



2 J cp 2 \ cp cpj 

— [Ar—+fl—] -SwMzSa+f3 + m^— ^ 

2 \ cp cpj Cfil 



(5.16) 



(vn \ 

+i— (^^tan/? + Li)\ 
. . ^ 2 ^ (5.17) 
(A^tan/3 + ^) j 

where C^* can be obtained from (5.16) upon the replacement a ^ a — 7r/2. 
Whenever convenient, we use the shorthand notation = sin^^vi/, c-y = cos 7, and 
= sin 7. The parameters Ar and fi only appear off-diagonal and they are multiplied 
with the associated fermion mass, the tau mass rrir- 

Using C = RrCR~, one obtains the couplings of the mass eigenstates ri and T2. In 
this regard, it is important to note that the coupling of the CP-odd Higgs boson to 
the lighter stau vanishes, C[t*,ti, A^] = 0. Therefore, we have not listed the process 
TiT* 7^" in Table 5.0. By the same token, there is also no s-channel exchange of A'^ 
in any of the annihilation channels. Note that this statement remains valid even after 
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the inclusion of radiative corrections: There is no induced mixing between hP[H^) and 
in absence of CP-violating effects in the SUSY sector. 

Let us now turn to the probably most interesting couplings in the context of rirf 
annihilation, namely, the ones of the lighter stau to and . The 'reduced' tiTi 
coupling reads 

(5.18) 

This is a complicated expression. However, if we choose m^o to be large, rrij^o » Mz, 
we can simplify (5.18) by using cos (/3 — a) = [cf. (5.14)], 

C°L[r*, h^] ^ Qc^^^ - sl,C2e}j Mlc^p -ml- "^XrS^e, ■ (5.19) 

Thereby, we make an interesting observation: In the decoupling limit (DL), the tiTi 
coupling becomes proportional to the left-right entry rrirXr of the stau mass-squared 
matrix (5.1) and to S2e^. Therefore, it comes as no surprise that the tit^ annihilation 
cross section into h?h? peaks at 9:^ = 7r/4 — the point of maximal tl-tr mixing — as can be 
seen, e.g., in Fig. 5.1b. Analogously, one finds that the tiTiH^ coupling is proportional 
to {Ar tan /3 + /i) S2e^ in the decoupling limit. Complementary, the TiT2h^ /H^ couplings 
exhibit in this limit the same combination of A, fi, and tan/3 as their riri counterparts 
but those terms are now multiplied by C2e^ instead of S29f 

After the above discussion, it is clear that there exists the possibility to enhance the 
total stau annihilation cross section atot — and thereby to decrease IV l/(c7tot^^) — by 
choosing a proper combination of large A, fj,, and tan p. In the remainder of this section, 
we will explore this possibility for two exemplary pMSSM scenarios. 

Before proceeding let us make some technical comments. Large values of the previ- 
ously mentioned parameters may well lead to large radiative corrections.^ In order to 
arrive at a proper loop-improved tree-level result, we re-evaluate the entire Higgs sec- 
tor using FeynHiggs. In particular, we have modified our generated matrix elements 
in a way that allows us to set all trilinear Higgs couplings to their loop-corrected val- 
ues.® Note that this goes well beyond a simple a Oeff prescription. Only then, we 
mostly find better agreement of our cross sections for stau annihilation into two Higgses 
with the ones computed by micrOMEGAs. The latter program uses CalcHEP [200] for 

^In this context, note that we introduce a large mt-m^^ ^ sphtting when choosing Ms ~ 1 TeV. 

*The author grateful to T. Plehn and M. Ranch for providing us, for cross-checking, with their 
implementation of a Fortran routine which calculates the Higgs self-couplings using the effective potential 
approach [204]. 
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the generation of the matrix elements. There, the trilinear Higgs self-couplings have 
been expressed in terms of mf^o , mjjo, and m^o which effectively reabsorbs a bulk of the 
radiative corrections [201]. We therefore think that we do slightly better whenever we 
encounter some disagreement between the mentioned cross sections. Though the overall 
effect on is typically small, it can be at the level of 20% (see below). Finally, it is 
well known that a large A parameter may lead to charge/color breaking minima (CCB) 
in the scalar MSSM potential; see, e.g., Ref. [205]. SuSpect performs some basic checks 
which we take into account to make sure that we do not violate the constraints associ- 
ated with CCB. We remark that our pMSSM scenarios are chosen such as to allow us to 
extract the important features of primordial stau annihilation in the most transparent 
way — without emphasis on naturalness considerations. 

In Fig. 5.3 we demonstrate the effect associated with a large TiTih? coupling by 
presenting the ^^^-dependence of (upper panel) and of the relative importance of the 
dominant thermally averaged cross sections, {(Tiv) / {atotv) , at a; = 30 (lower panel) for 
the pMSSM scenario associated with = 130 GeV, tan (3 = 50, mj^o = Ms = M3 = 
—A = 1 TeV, and 6M1 = M2 = /U = 1 TeV. In this scenario, m/jO stays in the range 
117 — 119 GeV and the lightest neutralino is bino-like with a mass of m~o = 169 GeV. 
Stau annihilation into heavy Higgses remains kinematically forbidden. The curves in 
the lower panel are associated with stau annihilation into h^h^, WW, tt, ZZ, 77, and 
7Z (at 9:^ = 80°, from top to bottom). As is evident, the annihilation into h^h^ is 
enhanced already well before 6^ = it/4:. At the peak position, a^ohov^ ~ 8.8 x 10'^ pb 
for Peflf = 10 GeV (no thermal average), which is still three orders of magnitude below 
the unitarity bound for inelastic s-wave annihilation, CTuVr = 8-7r/(mfjPefr) [206, 36]. 
Also the cross sections for stau annihilation into WW and ZZ are strongly enhanced 
towards 6^ = 7r/4 since the s-channel contribution of tiTi h^* VV becomes very 
important. At their respective peak positions, aww^r — 250 pb and azzi^r — 130 pb 
for Peff = 10 GeV. (Because of the dominance of the h^h^ channel, the corresponding 
maxima do not show up in Fig. 5.3 where {(Jiv) / {a totv) is shown.) By the same token, 
the cross sections of all (kinematically allowed) channels with a fermion-antifermion final 
state (e.g. rr) — which are subdominant in the scenario considered in Fig. 5.3 — experience 
an enhancement for 9^ ^ tt/4:. In total, there is an enhancement of {(Jtotv) that delays 
the thermal freeze out of the staus significantly, i.e., Xf ~ 33 for 9:^ ^ /A. As can be 
seen in the upper panel of Fig. 5.3, the decoupling yield is thereby reduced dramatically 
down to a minimum value of 1^ = 7.4 x 10~^^ for maximal left -right mixing of the staus. 

In the previous pMSSM examples, annihilation into final states containing heavy 
Higgs bosons is kinematically forbidden. We can allow for those channels by reducing the 
input value my^o . Indeed, scenarios in which all Higgs bosons are very light in conjunction 
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Figure 5.4: Analogous to Fig. 5.2 but for the pMSSM scenario associated with rriri ~ 
130 GeV, tan/3 = 50, and m^o = ^ = Ms ^ 6Mi = Afa.s ^ ~A = 1 TeV and for 
X = 30. The stau decoupling yield takes on its minimum value of IV ~ 7.4 x IQ-i^ at 
Or = 45°. The displayed stau annihilation channels are associated with the following 
final states: h°h°, WW, tt, ZZ, 77, and 7Z (at 6*? = 80°, from top to bottom). No 
lines are shown for 6^ < 4° where < my, . 
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Figure 5.5: Analogous to Fig. 5.3 but for the pMSSM scenario associated with m?^ = 
150 GcV, tan/3 = 50, m^o = 130 GeV, Ms = M3 = -A = 1 TeV, 3Mi = M2 = 
/i = 1 TeV and for x — 30. The stau decoupling yield reaches its minimum value of 
IV = 4.1 X 10^^^ at 9if — 7r/4. The displayed stau annihilation channels are associated 
with the following final states: h°h^, H"H°, bb, WW, ZZ, tt, and 77 (at 

9^ = 75°, from top to bottom). For an optimized presentation of those channels, the 
line indicating the relative importance of the sum of the displayed cross sections is 
scaled down by a factor of 1/2: (adispv) /2{atotv) ■ No lines are shown for 9r < 4° where 
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with large tan/3 have been studied in the literature, see, e.g., [207, 208] and references 
therein. We thus consider now the following pMSSM scenario: m^o = 130 GeV, m^^ = 
150 GeV, tan/? = 50, Ms = Ms = -A = 1 TeV, and 3Mi = Ms = /x = 1 TeV. In 
Fig. 5.4, the associated 6'^;-dependence of and of {(Jiv) / (atotv) at x = 30 for the now 
dominant channels is shown in a similar way as in Fig. 5.3; only the relative importance of 
the sum of the displayed cross sections is scaled down by a factor of 1/2, ((Tdispw)/2((Ttot^^)) 
to allow for an optimized presentation of the single dominant channels. Throughout the 
considered 0^ range, the masses of both CP-even Higgs bosons are relatively light and 
remain rather constant: mf^o = (118 it 1.5) GeV and mjjo = (128.5 it 1) GeV. Here 
the dominant annihilation channels are associated with the following final states: h^H^, 
h^h^, H^H^, bb, WW, ZZ, tt, and 77 (at 6^ = 75°, from top to bottom). As can 
be seen, stau annihilation into h^H^ is now more dominant than the one into h^h^ 
and also the H^H^ channel becomes important, where each of those channels is indeed 
associated with an (absolute) annihilation cross section {(Tiv) that peaks at 0:^ = 7r/4. 
Also the annihilation into bb is significant — a process which we will discuss in detail in 
the following section. In this respect, one should stress that all processes with s-channel 
exchange are here less suppressed by m^o in the respective propagator than in 
the previously considered scenarios. Note that the asymmetry of {(Jiv) / (atotv) of those 
dominant channels {hPH^, h?h^, H^H^, bb) with respect to a reflection at 9^ = 7r/4 
is dominantly caused by the ^^^-dependent modulation of the WW channel. As in the 
pMSSM scenario considered in Fig. 5.3, there is again an significant enhancement of 
(o'tot^^) that delays the stau freeze out such that Xf ~ 33 at 6^ ~ 7r/4. Thereby, the 
efficient annihilation into final state Higgses is accompanied by a significant drop in 
down to Y:f = 4.1 x IQ-^'' at 6*^ = 7r/4 as can be seen in Fig. 5.4. At this minimum, there 
is a 20% disagreement between 1~ from our calculation of stau decoupling (solid line) 
and the micrOMEGAs result Y~^^ (thin gray line) which is a consequence of the different 
treatments of the Higgs sector described above. 

Let us finally remark that the Higgs couplings to fermions and vector bosons as 
well as the Higgs self-couplings develop a strong dependence on m^o once we leave the 
decoupling regime (rriy^o < 200 GeV); for a comprehensive review see, e.g., Ref. [209].^ 
Changes in m^o can therefore be accompanied by shifts in the relative importance of 
the corresponding annihilation cross sections. This underlines the fact that the details 
in the Higgs sector may very well be crucial for the determination of the relic abundance 
of a long-lived ri. 

^The Higgs sector is also particularly sensitive to the mixing in the stop sector. In the considered 
pMSSM scenarios, \Xt \ = \At — ^cot[3\ ^ Ms which corresponds to the "typical-mixing scenario" [210]. 
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5.5 Resonant stau annihilation 

By inspection of Table 5.0 it becomes clear that primordial stau annihilation can also 
proceed resonantly via s-channel exchange of the heavy CP-even Higgs boson for 
mf{0 ~ 2m;p^. While the LEP bound on the stau mass m^f:-^ > 82 GeV [32] forbids 
to become on-shell {m^^^ ~ 140 GeV, e.g., [209]), the s-channel exchange of is 
absent^*^ because of C[ti, rj*, A'^] = (see Sec. 5.4). Again, our choice to work in the 
framework of the pMSSM proves to be very helpful. Since the resonance occurs for 
2mfj ~ rrifjo, one runs quickly into the decoupling limit in which mjjo is governed by 
the input parameter m^o according to the simple relation (5.13). This allows us to scan 
through the resonance easily. 

Let us explore resonant stau annihilation by considering the exemplary pMSSM sce- 
nario associated with m^^^ = 200 GeV, = 83° (i.e., a mostly 'right-handed' ri), 
tan/3 = 40, and —A = ^ = 4Mi = M2,3 = Ms = 1 TeV, for which we vary m^o 
(and thereby rrijjo) to scan through the resonance. Figure 5.5 shows the resulting rrijjo- 
dependence of 1^^ (upper panel) and of {civ) / (atotv) at x = 25 for the dominant anni- 
hilation channels (lower panel). Those channels are now associated with the following 
final states: bb, tt, tt, 77, h^h^ , WW, 'yZ, ZZ, and ti (at mjjo = 350 GeV, from top 
to bottom). In Table 5.0 all resonant channels can be identified. Close to the resonance 
condition Im^^ ~ ^hOj the most important processes are stau annihilation into bb and 
TT. This is because the couplings of those final state fermions to are tan /? enhanced: 
for tan/3 » 1, the ffH^ coupling ?71jS/3_q, tan/3 with f = b,T [203]. The (broad) 
peak associated with the resonance^^ already builds up for rrifjo > 2'm:f^ = 400 GeV. At 
zero relative velocity, this would be a region in which the resonance cannot occur. 
However, since ri is in kinetic equilibrium at the time of freeze out, resonant annihi- 
lation takes place already for 2m:^j < rrifjo [29]. For m^/o < 2mfj = 400 GeV, the 
processes containing s-channel exchange proceed with a slightly faster rate (if kine- 
matically allowed). The impact of the resonance on the thermal ri freeze out and 
the resulting 1~ is substantial. Since the total width of is Fj^o = (6 — 10) GeV for 
niHO = (300 — 500) GeV in the considered pMSSM scenario, the reduction of IV extends 
over a relatively large mjjo range. In this regard, note that Tjjo could be substantially 
larger had we not essentially decoupled all sfermions — except ti , T2, and Vr — by choosing 
Ms = 1 TeV. For rrij^o ~ 404 GeV, i.e., at the dip of the resonance, we find Xf ~ 33 
and a minimum stau decoupling yield of 1^ = 9.7 x 10^^^ (dark line). Thus, despite 

^"Even in absence of SUSY-induced CP violation, resonant annihilation via ^''-exchange may still 
proceed through ri-rg coannihilation. However, this scenario requires considerable fine-tuning in the 
stau mass-squared matrix since ri and T2 have to be nearly degenerate. 

^^Notice that we plot (atv) /{atotv) so that the actual shape of the resonance looks somewhat different. 
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Figure 5.6: Dependence of IV (upper panel) and of {(Jiv) / (atotv) at x = 25 (lower 
panel) on m^o for the pMSSM scenario associated with — 200 GeV, 9^ — 83°, 
tan/3 = 40, and -A = ^J. = 4Afi = il/2,3 = Ms = 1 TeV. In the upper panel, the 
dark line shows the stau yield obtained with our relic abundance calculation and the 
thin gray line the one obtained with micrOMEGAs. The stau decoupling yield takes on 
its minimum value of 1? = 9.7 x IQ-i^ at m^/o = 404 GeV. In the lower panel, the 
displayed dominant stau annihilation channels are associated with the following final 
states: bb, tt, tt, 77, WW , 7Z, ZZ, and ti (at m/^o = 350 GeV, from top to 

bottom). 
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the (stih) moderate value of tan/3 = 40, a significant reduction of 1^ is encountered. 
Indeed, 1^ can be even further suppressed for a larger value of tan (3. Let us remark 
that an accurate determination of 1^ in the resonance region requires to take special 
care of the bbH^ vertex. This coupling is well known to receive substantial radiative 
corrections for sizable values of tan/3. Therefore, we rely again on the computer tool 
FeynHiggs to compute all quark-antiquark-Higgs couplings and the total width Tf^o. 
Also the micrOMEGAs code takes special care of the bbH^ vertex. We therefore think 
that the difference between the yields shown in the upper panel of Fig. 5.5 reflects the 
theoretical uncertainty involved in the determination of T^jo as well as the bbH^ vertex. 

5.6 On the viability of a Ti-t^ asymmetry 

Given the strong bounds on the abundance of negatively charged ri from bound-state 
effects during BBN, i.e., from CBBN of ^Li and ^Be, it is natural to ask whether it 
is possible to have an excess of positively charged -rj^'s over negatively charged ti's. 
The generation of a particle-antiparticle asymmetry requires a departure from thermal 
equilibrium. Therefore, one might think that a Ti-f^ asymmetry can be produced at the 
time of the stau freeze out if the (slepton number violating) process tiTi tt occurs 
at a different rate than its conjugate counterpart. Such a situation might indeed occur 
if we allow for (CP-violating) complex values of the parameters A^, fi, and Mi^2 in the 
SUSY sector. However, the staus are still tightly coupled to Standard Model particles 
so that they remain in kinetic equilibrium with the primordial plasma. Therefore, any 
excess of rj* over ri arising will be washed out quickly by the inelastic scattering process 
T^T ^ tit}'^ Indeed, it is well-known [29] that processes of the latter type occur at 
much larger rates than the rates for the mutual annihilation of the decoupling particle 
species. The same argument given in [29] can be adopted to our case. At the time of 
freeze out, the reaction rates of interest can be estimated as 

fin ->■ TT : n:i=^n:i=^ar,7i~^rT ~ r^m|^e"^"'^i/'^fj?.^?:j^^^ , (5.20) 
r^T ^ ?ir : n^*n,a^.,_,^^- ~ T^/^m'J^e-^^i/^a^*r-.r,r , (5.21) 

—(*) 

since t^ is approximately Boltzmann distributed. For simplicity, we have treated the 
tau lepton r as a (still) relativistic species. By taking the ratio of (5.21) with respect 

Additional equilibrating processes are, e.g., tiW~ ^ tiW^ or tiH^ ^ tiH^ , which are however 
Boltzmann-suppressed. Also note that a lepton asymmetry of the order of the baryon asymmetry is 
expected because of charge neutrality of the Universe; cf. [211] and references therein. 
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to (5.20), 

e™-i/^ ~ 10^ for m^jT ~ 25 , (5.22) 

we find that the equilibrating process is by far more dominant. Here, we have used that 
cr;f-^;f^^^T- and a^'T-^Trr are not too different. In fact, both processes proceed at tree level 
exclusively via exchange so that one cannot decouple (5.21) from (5.20) by a simple 
adjustment of the neutralino mass spectrum. 

5.7 Exceptionally small abundances within the CMSSM 

We have shown above that the total stau annihilation cross section can be significantly 
enhanced. The thermal freeze out of ti's is thereby delayed such that their abundance 
prior to decay, 1%, is suppressed. In the following we focus on the CMSSM to see whether 
the effects discussed in Sects. 5.4 and 5.5 do appear also in models in which the pattern 
of soft-SUSY breaking parameters fulfills certain boundary conditions at a high scale. 
Note that we compute 1^ with micrOMEGAs in this section since coannihilation processes 
are not included in our relic density code. In addition, we employ SPheno 2.2.3 [157] 
for the computation of the mass spectrum and the low energy constraints associated with 
-B(6 — > 57) and the anomalous magnetic moment of the muon a^. Let us now proceed 
by discussing two exemplary CMSSM parameter scans. 

Figure 5.6 shows contours of constant in the (W'i/2)?7io) plane for tan/3 = 43, 
Aq = 0, and a negative sign of the /i parameter. The contour lines represent the values 

= 10~^^, 4 X 10~^^, 10~^^, and 4 x 10^^^, where darker shadings imply smaller values 
of 1%. The dashed lines are contours of m^f^ = 100, 300 and 600 GeV (from left to 
right). The light-shaded region at mi/2 ^ 450 GeV is excluded by the mass bound 
m/jO > 114.4 GeV from Higgs searches at LEP [32]. The white area indicates the region 
in which either correct electroweak symmetry breaking is not established (in the very 
upper left corner) or in which m~o < ra:^^. Since /x < 0, the plane is actually in tension 
because of (negative) SUSY contributions a^^^"^ to the anomalous magnetic moment of 
the muon, a^ = {g- 2)^/2. 

Figure 5.7 presents a scan over the ("^-1/2)^0) plane for tan/3 = 55, Aq = 2mo, and 
/X > with contours of = 4 x 10"^^ lO^^"^, 4 x lO"^'', 10"^^, and 4 x 10"^^ (darker 
shadings indicate smaller 1~ values) and = 100, 300, and 600 GeV (dashed lines, 
from left to right). The large light-shaded region in the lower left corner is excluded by 
the robust bound m^^ > 82 GeV [32] from collider searches of charged sleptons (or by 
the appearance of a tachyonic spectrum). The LEP Higgs bound m^o < 114.4 GeV [32] 
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Figure 5.7: Contours of (as labeled) in the (toi/2,™o) plane for tan/? = 43, 
Aq = 0, and fi < 0, where darker shadings imply smaller values. The dashed lines 
are contours of nir^ — 100, 300, and 600 GeV (from left to right). The light-shaded 
region at rni/2 < 450 GeV is excluded by the LEP bound mf^o < 114.4 GeV [32]. 
In the white area either m^o < rrir-^ or correct electroweak symmetry breaking is not 
established (in the very upper left corner), where the thin contours indicate the Higgs 
funnel in the NLSP region. Table 5.1 provides detailed information for the SUSY 
model represented by the point "A" that is indicated by the star. 



is situated within this region in close vicinity to its boundary for m-o ^ 400 GeV and 
is indicated by the solid line for rriQ > 400 GeV. In the region to the left of the dotted 
line, B{b 57) > 4.84 x 10"^ [212], which is in tension with the bounds from inclusive 
6 — > S7 decays. 

Let us now discuss some generic features of the stau yield within the CMSSM on the 
basis of Figs. 5.6 and 5.7. We note beforehand that our more general statements on the 
Ti NLSP region in the CMSSM are corroborated by a parameter scan over the following 
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Figure 5.8: Contours of Yf (as labeled) in the (toi/2,™o) plane for tan/3 = 55, 
Af) = 2mo, and /i > 0, where darker shadings imply smaller values. The dashed 
lines are contours of niri ~ 100, 300, and 600 GeV (from left to right). The large light- 
shaded region in the lower left corner is excluded by bounds from direct Higgs and SUSY 
searches (or by the appearance of a tachyonic spectrum) . In the region to the left of the 
vertical solid and dotted hues, rriho < 114.4 GeV [32] and B{b s-f) > 4.84x10-'' [212], 
respectively. In the white area, m^o < m^^. Table 5.1 provides detailed information for 
the SUSY models represented by the stars "B" and "C" (as labeled). 



range 

mi/2 = (0.1 - 6) TeV, tan/3 = 2 - 60, 

- Amo < Aq < 4mo, sgn/i = ±l. (5.23) 

In both figures an almost horizontal, narrow band of low appears in which 2mri ~ 
rriffo holds so that stau annihilation proceeds via resonant production of the heavy CP- 
even Higgs boson H^. We have marked points the centers of the respective regions with 
"A" and "B" for which we provide detailed information in Table 5.1. Given a present 
uncertainty of ~ 3 GeV in the determination of mi-fi [213], we note that the LEP Higgs 

^''Here, we disregard CMSSM parameter points in which SPheno flags an error in the spectrum calcu- 
lation. 
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bound has to be treated with some care. For example, a (vertical) m^jO = 112 GeV 
contour would be situated at mi 12 — 400 GeV in the resonance region of Fig. 5.7. 
Accordingly, one could consider the entire resonance region shown to be compatible 
with direct Higgs searches. However, due to the large value of tan/3 = 55, the bound on 
h ^ 57 is very strong so that a large part of the resonance region remains excluded by 
this constraint. In this regard, it is interesting to see (Fig. 5.6) that Im^^ ~ mj^^o also 
appears in the t\ NLSP region for lower values of tan /?. In the center of both resonance 
regions, the yield becomes as low as 1^ = 4.2 x 10^^^ (point A) and = 2.5 x 10^^^ 
(point B). Despite the heavier mass of the lighter stau (see Table 5.1), the suppression 
of is still more pronounced in Fig. 5.7 than in Fig. 5.6. This is because the bottom 
Yukawa coupling becomes larger with increasing tan /3, as discussed already in Sec. 5.5. 
In fact, annihilation into hh final states is in both cases by far the dominant process 
with relative importance of 76% (point A) and 87% (point B). The extension of both 
resonance regions is due to the total width of of respectively V no ~ 9.6 GeV (point 
A) and Fj^^o ~ 22 GeV (point B); note the logarithmic scales in Figs 5.6 and 5.7. We 
note in passing that the appearance of the H'^ resonance does not imply the absence of 
the neutralino funnel region which is indicated by the (unshaded) contour lines in the 
white area of Fig. 5.6 

Of course, the question arises whether the appearance of the resonance region is 
encountered more generically within the framework of the CMSSM. In principle, it is 
not easy to provide a simple quantitative connection between m^^ and m//o for arbitrary 
values of the CMSSM parameters. However, without emphasis on an overall applicability, 
a qualitative picture can be drawn. Let us start with the mass of the CP-odd Higgs boson 
m^o which can be written as [214, 215] 

m^o ~ 1 / sin^ /? {ml + 0.52m^/2 + - . . . ) . (5.24) 

Here, the ellipsis stand for contributions from the bottom and tau Yukawa couplings. 
For tan /? > 20, m^o ~ '^o ^ 2.5m^^2 ~ • • • ) s-^id the corrections from the bottom and 
tau Yukawa couplings become important so that m^o is driven towards lower values 
note that sin^ /? ~ 1 for tan /? > 20. Indeed, this property can be used to constrain tan /? 
from above by confronting m^o with the lower bound from LEP, m^o > 93.4 GeV [32]. 
On the other hand, for large mi/2) one also enters the decoupling limit of the MSSM so 
that m^o and mjyo will be nearly degenerate in mass; cf. (5.13). This can be also seen 
from the exemplary points presented in Table 5.1. Therefore, also m^Q will be driven 

^''The latter relation ignores contributions from A-terms which can be important but complicate the 
envisaged illustrative picture; for the derivation, we have used mt(mt) = 163 GeV in Eq. (2.25a) of 
Ref. [215]. 
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towards lower values for growing tan (5. Now, left-right mixing of the lighter stau for not 
too large values of tan (j is small within the CMSSM, ti ~ tr, so that approximately 
m~ ~ m-Q + 0.15m^^2 [215]. Therefore, 2mr^ < m^jo is the relation that holds usually 
in the region in which ti is the lightest Standard Model superpartner. However, for 
large tan /?, the contributions from the bottom Yukawa coupling in (5.24) can become 
strong enough (growing with rriQ [215]) to overcome any additional decrease of m^^ 
due to left-right mixing so that the resonance condition 2171^^ ~ mjjo can indeed be met. 
Nevertheless, from scanning over the CMSSM parameter range (5.23) it seems to us that 
the resonance condition 2m^j ~ nijjo is not easily realized in the part of the ti NLSP 
region in which ti-Xi coannihilations are negligible. Conversely, it is clear that relaxing 
the universality conditions for the soft-SUSY breaking masses at Mqut will make it easier 
to find parameter regions in which the resonance condition 2m:^j ~ niffo is satisfied. Of 
particular interest in this respect is the model with non-universal Higgs masses (NUHM) 
with niHi / / w-o at -^gut- There, one can adjust the input parameters in order 
to realize resonant stau annihilation. Indeed, this model is qualitatively the same as the 
class of pMSSM scenarios considered in the previous sections, where niHi and mjj^ are 
traded (at the low-scale) against my^o and /u by using the electroweak symmetry breaking 
conditions. 

Low values are also realized in the narrow vertical region around ~ 1-1 TeV 
in Fig. 5.7. At the representative point "C" of that region, 1~ = 2.2 x 10^^^ and the main 
stau annihilation channels are the ones into h^h^ (90%) and WW (6%); see Table 5.1. 
For larger values of mi/2, Yt exhibits its well known behavior and grows with m:^^. To 
the left of the = 4 x 10~^^ contour, the yield increases quickly since the annihilation 
into h^h^ becomes kinematically forbidden. Indeed, regions of low 1~ which are due to 
the aforementioned annihilation channels are a commonplace appearance in the CMSSM 
parameter space. They are found slightly above the lowest feasible values of mi/2, i-e., 
close to the boundary of the region which is excluded by direct Higgs and SUSY searches 
and where m^-^ > rrij^o still holds. This is because ri is light in that region since the SUSY 
particle spectrum scales with mi/2 (typically, mo <C mi/2 for ri NLSP). Moreover, we 
find that the LEP Higgs bound drops hardly below — 450 GeV for tan/3 > 40 

and mg < 100 GeV.^^ Due to a strong correlation between the gaugino mass parameter 
mii2 and the size of the /i parameter, /u^ ~ (1 — 3) m^y2 [216], the value of ^ in the 
experimentally allowed region is large. Recall from Sec. 5.4 that the TiTih? coupling 
is ~ sui29:fXr (m^o ^> Ml) so that \X^\ = \Aj. — ;utan/3| will become sizeable by 
increasing tan/3. This leads then to efficient stau annihilation into h?h? final states. 

^^The position of the LEP Higgs bound (which appears as a near to vertical line for low mo) is very 
sensitive to the value of mt. Lowering mt shifts the bound towards larger values of 
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Indeed, in those CMSSM regions, also | sin 2%| is maximized so that IV already starts to 
drop below the estimate (4.1) for tan/3 > 40. Note, however, that the left-right mixing 
of Ti within the CMSSM is somewhat constrained. Neglecting r- Yukawa contributions, 
the RG-evolution induced splitting reads — m-^ ~ 0.37my^ [215] and indeed ti 
remains mainly right-handed: By scanning over the parameter space, we typically find 
65° Ot ^ 115° and thus | sm.29:f\ < 0.75 in the ti NLSP region in which m^^^ > mj^o 
and m/io > 114.4 GeV holds. 

5.8 Prospects for collider phenomenology 

If a SUSY model with a long-lived ti of m^f^ < 0.7 TeV is realized in nature, the ti 
discovery potential will be promising at the LHC with a luminosity of 100 fb"^ [141]. For 
m^^ < 0.25 TeV (0.5 TeV), ri's can also be examined in precision studies at the ILC with 
a cm. range up to ^/s = 0.5 TeV (1 TeV). Once long-lived ri's are produced, one should 
be able to distinguish them from muons by considering the associated highly ionizing 
tracks and with time-of-flight measurements. One should then also be able to infer m^-^ 
from measurements of the ri velocity and its momentum [217] and complementary from 
(threshold) studies of the process e"'"e~ tit^ at the ILC. 

Both mechanisms leading to exceptionally small values come with testable pre- 
dictions: certain ranges of the stau-mixing angle 9^ together with large values of tan /3, 
|//|, and/or \At-\ and, in the case of resonant stau annihilation, also mjjo — 2m^^. In 
particular, the large stau-Higgs couplings lead to an enhanced production of light Higgs 
bosons in association with staus via e^e~ — > Trflh^ and 77 — > Trflh^ . The associated 
cross sections can then be relatively large at the ILC with a sufficiently high cm. en- 
ergy [218]. In addition, the above reactions with instead of hP in the final state 
can have also relatively large cross sections if and ti are sufficiently light. These 
reactions will then allow for an experimental determination of the stau-Higgs couplings 
and clarify whether its values are compatible with an extremely small value of [218]. 
Moreover, a measurement of rriijo pointing to m^o ~ 2m:^^ could be an experimental 
hint for resonant stau annihilation in the early Universe. 

Indeed, the scenarios considered could allow for a determination of both m^jO and 
rriijo already at the LHC. Because of the large values of tan /?, the dominant production 
mechanism for h? /H^ will be the associated production of the neutral Higgs bosons with 
bottom quark pairs, pp — > hhhP /H^; see, e.g., [219, 220, 221] and references therein. In 
fact, associated hhhP /H^ production with h^/H^ —>■ ^~^fi~ is considered as one of the 
most promising processes for measurements of rrij^o at the LHC despite the relatively 
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Table 5.2: Exemplary CMSSM points A, B, and C shown in Figs. 5.6 and 5.7. In 
addition to the quantities explained in the main text, values of the gluino mass nig and 
of the mass of the lighter stop m^_^ are given together with the relative importance of 
the dominant stau annihilation channels, Xf — rrir-^/Tf, and the decoupling yield Yp. 
For each point, we list gravitino dark matter scenarios with = 100 (50) GeV and 
associated values of the stau lifetime t^-^ , the non-thermally produced gravitino density 
h^, and the maximum reheating temperature T^^^. 
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9.6 
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-63 
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740 



76 
10 
9 
2 

30 
4.2 



10" 



GeV] 



100 
(50) 

5.7 X 10^ 
(7.5 X 10^) 

1.2 

(0.58) 

1.9 X 10^ 
(9.5 X 10^) 



B 



600 
748 
55 



197 
673 
80 



115 
390 
22 
666 
473 



1375 
1091 



87 

< 1 
11 

< 1 
30 
2.5 



100 

(50) 

6.5 X 10^ 
(6.4 X 10^) 

0.7 

(0.35) 
1.1 X 10^ 
(5.5 X 10^) 



C 



1138 
30 

55 



127 
739 

75 



117.9 
799 
41 

1262 
-164 



2446 
1757 



< 1 
90 

< 1 
6 

32 
2.2 



100 

(50) 

8.5 X 10^ 
(8.7 X 10"^) 

0.64 
(0.32) 

3.1 X 10^ 
(1.5 X 10^) 
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small h^/H^ /f^"*"/^" branching ratio [222]. In SUSY scenarios with a sufficiently 
light long-lived ri NLSP, these processes will be complemented by associated hhh^/H^ 
production with /H^ — > rirj*, where measurements of the invariant mass of the rirj* 
pair could potentially provide a unique way to infer m/jO and rrijjo at the LHC. In fact, 
hP /H^ — > Tirf will occur most prominently exactly in the regions associated with the 
exceptional 1^ values due to the enhanced stau-Higgs couplings. Having outlined these 
proposals, we leave a dedicated study for future work. 

Table 5.1 illustrates that the kinematical reach of both the LHC and the ILC could 
be sufficiently large to allow for the studies mentioned above. In none of the given points 
does m:f^ exceed 200 GeV so that rirj* pair production would already be possible at the 
ILC with y/s < 0.5 TeV. There, one could also produce TiT^h? final states in scenarios 
A and C. Even the condition rrijjo — 2m^^ could be probed in both scenarios A and B 
that allow for resonant stau annihilation. 



5.9 Implications for gravitino dark matter scenarios 

We have seen in this thesis that IV is subject to stringent cosmological constraints. 
Indeed, to decide on the cosmological viability of a SUSY model, one has to confront the 
associated 1^ values with those constraints. For gravitino LSP scenarios with unbroken 
R-parity, we have obtained restrictive cosmological constraints in Part II. In particular, 
in Sects. 4.1, 4.2.1, and 4.2.2 we have derived constraints and implications thereof under 
the assumption that >~ can be described by (4.3). However, while (4.3) is quite reliable 
for Ti ~ tr [31, 138, 1, 36], we have shown in the previous sections that 1~ (for a given 
m^^) can be about two orders of magnitude smaller than (4.3). 

Generally speaking, in this chapter we have shown that islands exist in which can 
be significantly below (4.3) even within the CMSSM and for a standard cosmological 
history. Thus, in gravitino dark matter scenarios with such exceptionally small IV values, 
our understanding of the cosmological constraints and the associated implications could 
change significantly. To demonstrate this point, let us indicate for which 1~ values the 
existing cosmological constraints are respected: 

• For 1~ < 10^^^, the upper limit on 1~ imposed by the non-thermal production of 
gravitinos in ri decays, 0-"'"^ < / ildm — given explicitly in (22) of Ref. [139] — is 
respected for < 500 GeV even if only a small fraction / = 0.01 of dark matter 
is assumed to originate from ri decays; cf. Fig. 13 of Ref. [139]. This applies equally 
to other scenarios with an extremely weakly interacting LSP — such as the axino 
LSP [223, 188, 224]— originating from n decays. 
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• For IV ^ 10~^^, the BBN constraints associated with effects of hadronic energy 
release on the primordial D abundance can be respected for ti ~ tr and m^^ up 
to 10 TeV independent of the ri lifetime; cf. Fig. 11 of Ref. [139]. For a sizable 
admixture of tl in ri , this IV constraint can become more restrictive in particular 
with the enhanced stau-Higgs couplings allowing for exceptionally small values. 
Nevertheless, these exceptional values are typically associated with m^^ < 300 GeV 
where the IV limit is significantly more relaxed: IV ^ 10^^^ for ri ~ tr. A 
tightening to IV ^ 10^^^ (10^^^) will then require an increase of (hadronic) i?vis 
by a factor of 10^ (10^). On the other hand, sufficiently degenerate and will 
always be associated with small values of -Evis and thereby with relaxed IV limits 
from energy release, even in the case of strongly enhanced stau-Higgs couplings. 

• For IV ^ 10~^^ (10~^^), the BBN constraints associated with effects of electromag- 
netic energy release on the primordial D (•^He) abundance can be respected inde- 
pendent of the Tl lifetime; cf. upper panels of Fig. 12 (100 GeV < m^f:-^ < 10 TeV) 
of Ref. [139] and Figs. 14 (m^^ = 100 GeV) and 15 (m^^ = 300 GeV) of Ref. [130]. 

• For IV < 2 X 10"^^ (2 x lO^^^^ ^ 2 x 10"^^), the BBN constraints associated with 
bound state effects allowing for CBBN of ^Li and ^Be can be respected even for 

> 10^ s; see Fig. 2.8. Recall, that these values correspond to upper limits on 
the primordial fractions of ^Be/H and ^Li/H of 2.1 x 10^^^ and to the generous 
range 10^^^ 10^^'^, respectively. 

Thus, the SUSY models which come with thermal relic stau abundances of IV ^ 2 x 10~^^ 
can respect each of those cosmological constraints independently of the stau lifetime if 
a primordial ^Li/H abundance of about 10"^'^ is viable. In particular, the limit (4.8) of 
"^Ti ^ 6 X 10^ s and its implications discussed in Chapter 4 are then no longer valid even 
for a standard cosmological history with primordial temperatures of T > Tf. Thereby, 
the regions with IV < 2 x 10^^^ are associated with particularly attractive gravitino dark 
matter scenarios: 

• The gravitino mass can be within the range 0.1 < mg < m;^^ for which its kinemat- 
ical determination could be viable [154, 190, 191]. Together with measurements 
of rriri "^rn a kinematically determined mg would allow one to measure the 
Planck scale Mp at colliders [154, 190, 191]. Indeed, an agreement of the Mp value 
determined in collider experiments with the one inferred from Newton's constant 
Gn would support the existence of supergravity in nature [154]. 

• For sufficiently close to m^i, the spin-3/2 character of the gravitino becomes 
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relevant so that it could be probed in principle by analyzing the decays ri 
Gtj [154]. 

• With 1- < 2 X 10^^^, 0-'^^ is neghgible so that basically all of ^l^m can be 
provided by gravitinos from other sources such as thermal production. Indeed, if 
also gravitino production in decays of scalar fields such as the inflaton [121, 123] 
is negligible, reheating temperatures of > 10^ GeV could become viable for 
TJiQ ~ 100 GeV and not too heavy gaugino masses; see, in particular. Sec. 3.4. This 
would mean that thermally produced gravitinos could provide the right amount of 
dark matter and that thermal leptogenesis (with Tr > 10^ GeV as a benchmark 
value [134, 135]) would be a viable explanation of the cosmic baryon asymmetry, 
i.e., there would be no gravitino problem. 

• With a kinematically determined m^, one would be able to probe the reheating 
temperature Tr at colliders and thereby the viability of thermal leptogenesis [119]. 

• For T:^^ > 10^ s, the small If values could still allow for the primordial catalysis of 
^Li and ^Be in agreement with existing astrophysical observations; see Sec. 2.6. 

Table 5.1 illustrates that gravitino dark matter scenarios of the type discussed above 
can even be accommodated within the CMSSM. For gravitino masses of 50 GeV and 
100 GeV, we list the associated values of r:^, , of ^^^h"^, and of the maximum reheating 
temperature T^'^^ under the assumption that other gravitino sources can be neglected. 
The stau lifetime t^^ is given in (4.4) and the T^^^ values imposed by 17-^ /i^ < 0.126 can 
be inferred from (4.12). At each CMSSM point and for both mg values, ri is very long 
lived, T:^: > 10® s, and gravitino production from ti decays is negligible, fi-^^/i^ < 10~^. 
In all cases, the gravitino mass = 100 GeV is sufficiently close to m:^^ so that 
the spin-3/2 character of the gravitino can in principle be probed [154]. A reheating 
temperature of Tr > 10^ GeV is viable only for the points A and B with significantly 
below 1 TeV, i.e., at the points at which resonant stau annihilation leads to the reduction 
of Y^. Because of r^^ > 10® s, the hmit from CBBN of ^Be is at Fj^^^ ~ 2 x lO'^^ 
for each point as can be inferred from Fig. 2.8. This bound disfavors point A while the 
points B and C are associated with 1^ values very close to this limit and thereby with 
^Be/H (®Li/H) values of about 2.1 x IQ-^^ (lO^iO), 
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In this thesis we have worked out the cosmological impUcations of a long-hved electro- 
magnetically charged massive particle species also called CHAMP. Our working hy- 
pothesis has been that X possesses a weak-scale mass mx > O (100 GeV) and typically 
a lifetime rx ^ 1 s. We have assumed that the temperature T of the early Universe was 
high enough so that X has achieved chemical equilibrium with the primordial plasma. 
Then, following a standard cosmological evolution, X experienced a thermal freeze-out 
once T < mx/25. This makes X to what is called a thermal relic (prior to its decay). 

BBN with a long-lived CHAMP 

We have started our investigation with a brief introduction into the framework of stan- 
dard Big Bang Nucleosynthesis (SBBN) where we also have given account to some of 
the latest measurements from which primordial light element abundances are inferred. 
In a simplified discussion of X-decoupling we have argued that its expected cosmological 
abundance prior to decay reads lO'^^ < Yx < lO-^^^mx/lOOGeV). We have recalled 
that the long-lived CHAMP scenario is strongly constrained by BBN limits on electro- 
magnetic and hadronic energy release in the X-decay. For a reliability check on hadronic 
BBN constraints we have worked out the Coulomb stopping power of charged hadrons 
in the plasma. In particular, we have developed on a refined approach taking into ac- 
count peculiarities in the plasma-screening of the Coulomb interaction and paying close 
attention to the velocity dependencies of the cross sections. We find reasonably good 
agreement with the treatment used in Ref. [41] from which we have incorporated the 
associated constraints. 

Subsequently, the effects of X~ on BBN due to its binding onto the light nuclei N 
have been considered. Given that the X~ -catalysis of thermal nucleosynthesis reactions 
had only been discovered recently [79], we have laid out in detail the central points 
of CBBN. Using the variational approach, we have obtained ground state energies for 
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bound states (NX~) by taking into account the finite nuclear charge radius of N. This 
leads to a reduction (in magnitude) from the naive point-like Coulomb values, e.g., for 
(^HeX~) by 13% and for (^BeX~) by 60% which has also been confirmed upon numerical 
solution of the Schrodinger equation. For the examples of ("^HeX") and (^BeX~) also 
the complete spectrum for n < 3 has been computed. We have further obtained the wave 
functions for the N-X~ continuum. This has allowed us to calculate the cross sections 
for (NX~) photo-dissociation, (uphv), and radiative recombination, (arecv), including the 
finite charge radius correction and taking into account recombinations into IS as well as 
2S states. Those rates (per particle pair) are important since they control the fractional 
bound state abundance and thus the timing and efficiency of CBBN. For example, for 
(^HeX-*) and (^LiX") we find a reduction of ((Trcc'w) from the hydrogen-like case by 17% 
and by 74%, respectively. 

From observations of beryllium in Population II halo stars at very low metallicities, 
we have extracted a nominal upper limit on primordial beryllium of ^Be/H < 2.1 x 10^^^. 
This limit allows one to set constraints on models in which the primordial A = 8 divide 
is bridged by catalytic effects. Considering the primordial catalysis of ^Be [63], we have 
derived rx-dependent upper limits on the X~-yield prior to decay, F-^f . For a typical 
relic abundance > 3 x 10"!^ (IQ-^''), we find that this ^Be limit translates into 
an upper limit on the X~ lifetime of rx ^ 6 x 10^ s (10^ s). Furthermore, we have also 
worked out the catalytic production of ^Li which, depending on the adopted upper limit 
on primordial ^Li, gives rise to similar bounds. 

We have clarified that the presence of (pX^) bound states cannot relax the Y^l'^ 
limits at long lifetimes tx in any substantial way. Indeed, we have shown explicitly by 
solving the associated full set of Boltzmann equations that late-time effects of (pX^) 
bound states can affect the lithium and beryllium abundances synthesized at T ~ 8 keV 
by not more than a few percent. Any substantial formation of {pX~) at T ~ 0.7 keV 
is immediately intercepted by the very efficient charge exchange reaction of {pX~) with 
■^He. This comes as no surprise given the large size of the {pX~) system ~ 30 fm and 
the fact that the proton deconfinement probability approaches unity already for a ^He- 
X~ distance of ~ 95 fm. In particular, we find that the fractional density of protons 
in bound states does not exceed the level of ~ 10"^ for Yx- < ^4He- By the same 
argument, the ^Be yield also remains unaffected by late-time catalysis. Thus, we find 
that the possibility of allowed islands in the parameter region with typical Y^'^S and large 
Tx — which was advocated in Ref. [90] — does not exist. 
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The gravitino-stau scenario 

In the second part of this thesis we have considered cosmological constraints and their 
impUcations for models in which the gravitino is the LSP and the stau is the NLSP. We 
have first focused on G as a dark matter candidate. Using the full gauge-invariant result 
for the thermally produced gravitino abundance to leading order in the Standard 
Model gauge couplings [108, 119], we have studied bounds on the reheating temperature 
Tr from the constraint ^ ^dm- In particular, taking into account the dependence of 
on the masses of the gauginos has allowed us to explore the dependence of the Tr 
bounds on the gaugino-mass relation at the scale of grand unification A/gut- We have 
explicitly studied the effect of G regeneration during a post-inflationary perturbative 
reheating phase. Thereby, we have made contact between the notion of Tr as the initial 
temperature of the radiation-dominated epoch in the analytical expression (3.6) for the 
G abundance and the definition of Tr in terms of the decay width T^p of the inflaton. 

Applying the rx-dependent upper limits on Y^^J: derived from the primordial catal- 
ysis of ^Be and ^Li in Part I, we have analyzed the emerging constraints in the gravitino- 
stau scenario, i.e., for = X~ . For typical values (4.3) of the stau NLSP yield after 
decoupling, the ^Be and ^Li constraints have been found in close vicinity to each other 
so that they lead to the same implications. For example, for = 10 GeV, the CBBN 
constraints impose the lower limit m^^ > 400 GeV with rising tendency for growing 
rriQ. For ti being the lightest Standard Model superpartner such a limit directly af- 
fects the testability of those SUSY scenarios at future colliders. Furthermore, for a 
primordial limit of ^Li/H ^ 6 x 10^^^ the calculated ^Li abundance drops below this 
observational bound only for r^j ^ 6 x 10^ s (likewise for ^Be). Taken at face value, 
we find that this constraint translates into a lower limit on the gaugino mass parameter 
"^1/2 ^ 0-87 TeV (mg/10 GeV)^/^ in the entire natural region of the CMSSM parameter 
space. This limit implies a restrictive upper bound 7r < 5 x 10^ GeV (mg/10 GeV)^/^. 

Using exemplary (7711/25^-0) CMSSM planes where we explicitly compute Y-^'^ in 
every point, we have further explored gravitino dark matter scenarios and the associ- 
ated Tr bounds for niQ > 10 GeV and for temperatures as low as 10'' GeV. Taking 
into account the ^Li CBBN constraint as well as the constraints on electromagnetic and 
hadronic energy injection from ri-decays, we have illustrated that in the considered re- 
gions of the CMSSM parameter space Tr < lO'^ GeV indeed is the highest cosmologically 
viable temperature of the radiation-dominated epoch in case of a standard thermal his- 
tory of the Universe. Moreover, in the ri NLSP region the lower bound on mi/2 typically 
implies a very heavy superparticle mass spectrum where, e.g., rrig < 2.5 TeV can be well 
excluded and which makes such scenarios hard to probe at the LHC. The bound on Tr 
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imposes a serious constraint for inflation models. Moreover, thermal leptogenesis seems 
to be strongly disfavored in the considered regions of the CMSSM parameter space. 

With late-time entropy release, the obtained limit Tr < 10'' GeV can be relaxed. 
For example, the dilution of the thermally produced gravitino yield by a factor of 10 
relaxes the Tr bound by about one order of magnitude in regions where i7 dominates 
fig. In the case of entropy production after NLSP decoupling, the yield of the NLSP 
prior to its decay, InlsPi is reduced so that the BBN constraints can be weakened. 
Although the ^Li bound is persistent, we find that it disappears provided Inlsp is 
diluted by a factor of A > 10^. We have discussed the viability of thermal leptogenesis 
in a cosmological scenario with entropy production after NLSP decoupling. We find 
that successful thermal leptogenesis can be revived in generic regions of the CMSSM 
parameters space for A/ri ~ Tr > 10^^ GeV and A > 10^, where Mri is the mass of the 
lightest among the heavy right-handed Majorana neutrinos. There, the collider-friendly 
r NLSP region with m^f < 250 GeV reopens as a cosmologically allowed region in the 
CMSSM with the gravitino LSP. 

Thermal relic stau abundances 

In the final part of this thesis we have carried out a thorough study of primordial stau 
annihilation and the associated thermal freeze-out. Taking into account the complete 
set of stau annihilation channels within the MSSM with real parameters for cases with 
negligible sparticle coannihilation, the resulting thermal relic ri yield Y-^'^ has been 
examined systematically. While we have often (implicitly) focused on the ri ~ tr case 
in Part II by employing (4.3), we have investigated cases in Part III in which ri contains 
a significant admixture of tl including the maximal mixing case as well as ri — tl. 

We find that the variation of the stau mixing angle 0^ affects the relative importance 
of the different annihilation channels significantly but not necessarily the resulting 1^'^'^ 
value for relatively small values of tan (3. By increasing tan however, we encounter 
a dramatic change of this picture for large absolute values of the Higgs-higgsino mass 
parameter /x and/or of the trilinear coupling Ar, which are the dimensionful SUSY pa- 
rameters that govern simultaneously stau left-right mixing and the stau-Higgs couplings: 
Stau annihilation into and H^H^ can become very efficient (if kinemati- 

cally allowed) so that Y~^^ can decrease to values well below 10~^^. The scalar nature 
of Tl allows those parameters to enter directly into the annihilation cross sections. This 
mechanism has no analogue in calculations of the thermal relic density of the lightest 
neutralino x5- 
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The stau-Higgs couplings are crucial also for the second Y~^^ reduction mechanism 
identified in this work: Even for moderate values of tan /3, we find that staus can anni- 
hilate very efficiently into a bb pair via s-channel exchange of the heavy CP-even Higgs 
boson provided the MSSM spectrum exhibits the resonance condition 2m:^-^ ~ m^i^o. 
We have shown explicitly that the associated Y-^'^ values can be below 10^"*^^ as well. 
This mechanism is similar to the one that leads to the reduction of the Xi density in the 
Higgs funnel region in which neutralino annihilation proceeds at the resonance of the 
CP-odd Higgs boson A'^. 

We have worked with an effective low energy version of the MSSM to investigate 
the 6';p-dependence of Y-^'^ and the two 1^^^-reduction mechanisms in a controlled way. 
In addition, we have shown that the considered effects can be accommodated also with 
restrictive assumptions on the soft-SUSY breaking sector at a high scale. Within the 
CMSSM, we encounter both mechanisms each of which leading to Y-^'^ ~ 2 x 10~^^ in 
two distinct regions of a single {'mi/2, mo) plane. 

We have discussed possibilities to probe the viability of the presented l^^'^-reduction 
mechanisms at colliders. While a m^o measurement pointing to mjjo 2± 2m^^ would 
support resonant primordial stau annihilation, studies of Higgs boson production in 
association with staus, e+e^ (77) TiT^h^,TiT^H^ could allow for an experimental 
determination of the relevant stau-Higgs couplings, for example, at the ILC. Moreover, 
we have outlined that associated hhh^/H^ production with / tit^ has the 
potential to allow for a determination of both nii^o and rrijjo at the LHC if a SUSY 
scenario with large tan /? and large stau-Higgs couplings is realized. 

With the obtained small Y-'^'^ values, even the restrictive constraints associated with 
CBBN could be respected so that attractive gravitino dark matter scenarios could be 
revived to be cosmologically viable even for a standard cosmological history. Within this 
class of models, collider evidence for supergravity, for the gravitino being the LSP, and 
for high values of the reheating temperatures of up to 10^ GeV is conceivable, which could 
thereby accommodate simultaneously the explanation of the cosmic baryon asymmetry 
provided by thermal leptogenesis and the hypothesis of thermally produced gravitinos 
being the dark matter in our Universe. 
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